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Abstract 

The  Non-Intrusive  Load  Monitor  (NILM)  measures  equipment  performance  by  measuring 
and  analyzing  the  source  power  to  the  equipment  at  a  single  point  in  the  electrical  system. 
Previous  studies  have  proven  the  usefulness  of  the  NILM  system  in  characterizing  the  state 
of  mechanical  systems  onboard  U.S.  Coast  Guard  vessels  and  at  the  U.S.  Navys  Land  Based 
Engineering  Site  (LBES)  in  Philadelphia,  Pennsylvania. 

This  thesis  seeks  to  augment  the  NILM  system  by  exploring  a  more  user  friendly  Graph¬ 
ical  User  Interface  (GUI)  to  allow  shipboard  crews  to  utilize  the  NILM  while  in  operation. 
Previous  applications  of  NILM  required  post-event  data  analysis  in  the  laboratory.  An 
additional  monitor  was  installed  on  the  Low  Pressure  Air  Compressor  (LPAC)  #1  at  the 
LBES  facility  to  investigate  abnormalities  detected  in  the  operation  of  LPAC  # 2  by  previous 
research.  The  ability  of  the  NILM  to  function  at  the  highest  levels  of  the  electrical  distri¬ 
bution  system  was  also  explored  at  the  LBES  facility  with  the  installation  of  two  additional 
NILM  systems  on  the  main  switchboards  supplying  power  to  the  auxiliary  system  loads.  Fi¬ 
nally,  a  brief  overview  of  the  analysis  software  of  the  Multi-Function  Monitor  (MFM),  a  key 
component  in  modern  ships  Zonal  Electrical  Distribution  Systems  (ZEDS),  is  presented  to 
explore  the  possibility  of  the  NILM  and  MFM  systems  operating  in  conjunction  to  improve 
the  operation  of  future  ZEDS. 

Thesis  Supervisor:  Robert  W.  Cox 

Title:  Assistant  Professor  of  Electrical  and  Computer  Engineering 
University  of  North  Carolina  at  Charlotte 

Thesis  Supervisor:  Steven  B.  Leeb 

Title:  Professor  of  Electrical  Engineering  and  Computer  Science  &  Mechanical  Engineering 
Thesis  Supervisor:  Patrick  Hale 

Title:  Director,  Systems  Design  and  Management  Fellows  Program 
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Chapter  1 


Introduction 


1.1  NILM  System 


The  Non-Intrusive  Load  Monitoring  (NILM)  system  consists  of  two  major  components,  an 
electrical  enclosure  that  processes  sensor  inputs  and  a  computer  that  analyizes  the  output  of 
the  electrical  enclosure.  On  a  3-phase  system,  the  inputs  to  the  electrical  enclosure  are  one 
current  transducer  (CT)  and  two  voltages  taps  use  to  determine  the  voltage  between  two 
phases.  The  CT  is  normally  installed  on  the  A-phase  and  the  voltage  taps  are  installed  on 
the  B-  and  C-phases.  The  electronics  in  the  electrical  enclosure  collects  these  measurements 
and  sends  them  to  the  computer  where  the  real  and  reactive  powers  are  calculated.  These 
values  are  sent  to  the  computer  which  stores  these  values  into  hourly  log  files.  These  files 
can  be  analyzed  to  determine  the  operating  characteristics  of  the  shipboard  systems  and  to 
perform  diagnostics  on  these  systems. 

NILM  research  has  be  conducted  at  the  Massachusetts  Institute  of  Technologys  (MIT) 
Laboratory  for  Electromagnetic  and  Electronic  Systems  (LEES)  for  over  two  decades.  The 
main  focus  of  the  research  for  shipboard  applications  to  this  point  has  been  to  understand 
how  the  dynamics  of  various  shipboard  systems  can  be  monitored  using  NILM  and  their 
operating  modes  and  failure  modes  detected.  The  logging  of  the  observed  data  occurred  in 
real  time  as  the  events  happened  but  the  analysis  was  conducted  at  a  later  time. 
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1.2  Motivation  for  Research 


1.2.1  NILM  Graphical  User  Interface 

The  nature  of  the  past  research  was  focused  primarily  on  understanding  the  systems  being 
monitored.  There  is  now  a  significant  level  of  knowledge  on  a  few  specific  systems  such 
that  the  focus  can  be  shifted  to  that  of  making  the  system  useable  and  meaningful  to  the 
customer,  in  this  case  the  ship's  crew.  In  order  to  do  this  effectively,  the  analysis  of  the 
events  that  occurred  and  their  classification  had  to  be  automated. 

A  second  aspect  that  was  explored  was  the  interface  through  which  the  operator  interacts 
with  the  NILM  computer.  Traditionally,  the  interface  with  the  NILM  system  was  a  simple 
menu  window  with  limited  functionality.  A  significant  number  of  operations  required  exiting 
to  the  Linux  command  line  interface  where  standard  Linux  commands  could  be  used.  The 
processing  of  data  involved  utilizing  scripts  written  in  PERL  and  MATLABrMprograms, 
both  of  which  were  normally  written  individually.  It  is  unreasonable  to  assume  that  the 
ship’s  personnel  for  the  U.S.  Navy  or  the  U.S.  Coast  Guard  will  have  sufficient  time  or 
training  to  process  the  data  in  this  fashion.  Also,  this  wrould  not  lend  itself  to  real  time 
processing  of  the  data  and  indication  of  the  operation  of  the  system  being  monitored. 


1.2,2  NILM  at  the  Land  Based  Engineering  Site 

The  non-intrusive  nature  of  the  NILM  system  has  had  limited  examination  in  a  shipboard 
environment  up  to  this  time.  In  order  to  understand  the  operation  of  a  particular  system, 
the  NILM  equipment  was  normally  installed  electrically  close  to  the  system  being  monitored 
to  minimize  interference  and  noise  in  the  signals.  When  installed  in  this  manner,  multiple 
system  would  be  required  to  monitor  many  systems  onboard  the  ship.  This  would  not  be 
truly  non- intrusive  as  these  systems  would  require  many  interfaces  and  NILM  equipment 
located  throughout  the  engineering  spaces. 

The  Low  Pressure  Air  Compressor  (LPAC)  monitored  at  LBES  was  found  in  previous 
work  to  have  an  unexplained  abnormality  in  its  operation  [3].  In  order  to  determine  the 
cause  of  this  abnormal  signal  in  the  operation  of  the  LPAC,  further  investigations  are 
required  and  the  installation  of  a  NILM  on  the  second  LPAC  for  comparison. 
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1.3  Objectives  and  Outline  of  Thesis 


The  research  presented  in  thesis  is  a  continuation  of  research  conducted  by  LCDR  William 
Greene,  USN  [18],  LCDR  Thomas  D.  McKay,  Jr.,  USN  [23],  LCDR  Patrick  Bennett,  USN 
[3]  and  LT  Mark  Fiber,  USCG  [32].  LCDR  Greene  was  on  the  team  that  installed  the  first 
set  of  two  NILM  sensors  at  the  LBES  Facility  on  the  LPAC  and  the  Fuel  Oil  (FO)  Pumps 
supplying  the  LM2500  gas  turbine  main  engines.  He  also  presented  the  first  look  at  using 
NILM  in  conjuction  with  the  Multi-Function  Monitors  (MFMs)  installed  on  ships  equipped 
with  a  Zonal  Electrical  Distribution  System  (ZEDS).  LCDR  McKay  installed  a  third  NILM 
at  the  LBES  facility  on  the  Universal  Engine  Controller  (UEC).  Analysis  of  the  data  from 
the  three  NILM  systems  allowed  major  events  which  occur  in  the  engineering  plant.  LCDR 
Bennett  continued  the  study  of  the  LBES  facility,  installing  two  additional  NILM  systems 
on  the  Lube  Oil  (LO)  and  Gas  Turbine  Generators  (GTGs)  to  provide  additional  data  for 
event  change  detection.  LT  Piber  researched  using  NILM  to  diagnose  and  solve  problems 
associated  withthe  Coast  Guard  Cutters. 

The  objective  of  this  thesis  is  to  further  research  into  areas  associated  with  the  two  main 
areas  of  the  MIT  LEES  NILM  group,  the  U.S.  Coast  Guard  (USCG)  Cutters  and  the  LBES 
Facility.  A  Graphical  User  Interface  (GUI)  was  investigated  for  use  with  the  NILM  system 
and  was  specifically  designed  for  use  with  the  CHT  system  onboard  the  USCG  Cutters. 
At  the  LBES  Faeilty,  research  into  the  incorporation  of  the  NILM  system  into  the  MFM 
systems.  The  MFM  algorithm  if  examined  in  detail  and  specific  functions  related  to  NILM 
analysis.  A  second  NILM  was  installed  on  the  remaining  LPAC  at  the  LBES  facility  and 
the  systems  monitored  to  determine  the  cause  of  the  abnormal  operation  noted  be  Bennett 
[3].  Two  additional  NILM  systems  were  installed  on  the  main  electrical  distribution  panels 
for  the  LBES  facility. 
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Chapter  2 


Graphical  User  Interface  (GUI) 
For  the  NILM 


Previous  research  with  the  NILM  system  focused  on  characterizing  the  operation  of  specific 
systems  such  that  their  operating  states  can  be  detected  and  used  in  diagnostics  to  determine 
normal  and  abnormal  operating  conditions.  For  example,  detecting  whether  a  pump  is 
operating  in  slow  speed  or  fast  speed.  As  the  systems  were  monitored  by  the  NILM,  the 
operating  conditions  that  were  logged  by  the  operators  were  compared  to  the  recorded  data. 
This  provides  a  means  correlating  the  electrical  power  signals  to  known  events.  Once  the 
events  and  their  associated  electrical  signals  have  been  correlated,  the  power  signals  can  be 
monitored  and  the  operating  state  of  the  equipment  detected. 

The  installation  of  the  NILM  systems  for  the  U.S.C.G.  Cutters  occurred  while  the  ships 
were  in  port.  Data  collection  and  storage  occurred  while  the  ships  operated  at  sea.  Upon 
return  to  port,  the  data  was  downloaded  and  analyzed  in  a  laboratory  setting.  The  NILM 
system  is  effective  in  detecting  abnormal  modes  of  operation  [32,  24]  however  the  notification 
to  the  crew  is  always  delayed  by  the  need  to  analyze  the  data  after  returning  to  port.  An 
automated  classification  algorithm  was  developed  by  LCDR  Proper  [34].  This  thesis  shows 
the  development  of  a  GUI  for  the  Collection  Hold  and  Transfer  (CHT)  system. 

2.1  CHT  System  Description 

The  CHT  system  onboard  the  U.S.C.G.  270  foot  Famous  Class  Cutters  is  a  vacuum  assisted 
collection  system.  The  system  consists  of  a  3G0-  gallon  collection  tank,  2  sewage  discharge 
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pumps,  2  vacuum  pumps  and  the  associated  piping  and  valves  necessary  for  controlling  the 
flow.  The  system  collects  waste  from  22  locations  within  the  ship.  The  vacuum  pumps 
maintain  a  vacuum  within  a  normal  operating  range  of  12  to  19  inches  Hg.  The  sewage 
discharge  pumps  maintain  the  level  within  the  collection  tank  as  indicated  by  3  level  probes 
[17].  A  simplified  system  diagram  is  shown  in  Figure  2-1. 


To  seal 
tank 


To  holding 
tank 


2.2  CHT  GUI  Description 

The  CHT  GUI  was  programmed  in  MATLAB™to  take  advantage  of  the  advanced  math- 
ematics  and  GUI  building  features.  The  CHT  GUI  program  was  designed  in  a  modular 
fashion  with  a  systems  perspective.  The  GUI  would  not  be  a  stand  alone  program,  it  would 
take  the  output  of  the  classification  program  developed  by  LCDR  Proper  [34]  and  present 
the  information  to  the  user. 

2,2.1  Design 
Requirements  Identification 

The  requirements  for  the  CHT  GUI  were  to  provide  an  interface  for  the  ships  crew  to  utilize 
the  NILM  software  without  requiring  extensive  training.  The  design  for  the  GUI  needed  to 
be  intuitive  such  that  crewmembers  could  easily  understand  the  information  quickly.  The 


22 


NILM  system  would  be  intalled  in  an  auxilary  engineering  space  onboard  the  ship,  not  in 
an  office.  The  crewmembers  needed  to  be  able  to  operate  the  GUI  in  this  environment  in 
the  course  of  normal  duties. 

Design  Philosophy 

A  modular  design  was  developed  with  a  clear  definition  of  the  interfaces  between  modules. 
Each  module  was  designed  to  perform  a  specific  function  and  provide  the  resultant  data  or 
indications  to  the  other  modules  through  the  interface  variables.  The  GUI  was  designed 
specifically  for  the  CHT  system  onboard  the  IJSCG  cutters,  but  with  a  modular  design, 
modules  could  be  changed  to  make  it  applicable  to  other  systems.  The  modular  system 
would  also  facilitate  upgrades  to  the  system  as  better  processing  algorithms  are  developed. 
A  symplified  design  spiral  approach  was  used  in  the  development  of  the  GUI. 

2.2.2  Implementation  of  GUI 

The  design  of  the  GUI  was  performed  with  utilizing  concepts  of  human  interface  design. 
Large  buttons  were  included  that  could  be  utilized  via  the  touch  screen  and  would  not 
require  fine  motor  controls  such  as  drop  down  menus  do.  Since  the  GUI  would  be  located 
in  an  engineering  space,  fine  controls  would  make  the  interface  more  difficult  for  the  crew 
to  use. 

The  GUI  has  multiple  screens  varying  in  the  depth  of  the  information  presented.  The 
intitial  screen  contains  a  summary  of  the  system  operation  and  have  all  of  the  information 
needed  during  normal  operation  and  is  shown  in  Figure  2-4.  A  second  window  provides 
more  detailed  information  in  the  event  a  fault  is  detected.  The  summary  screen  contains  a 
simplified  diagram  of  the  system  which  is  updated  by  the  program  to  reflect  the  status  of 
the  system.  The  color  of  the  equipment  is  changed  as  the  status  changes  the  color  scheme 
is  intuitive  and  shown  in  Table  2.1. 


Table  2.1:  Color  Scheme  for  System  Diagram 


None 

Equipment  secured  (pumps  only) 

Green 

Normal  Operation 

Yellow 

Potential  Fault  Detected 

Red 

Fault  Detected 
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Figure  2-2:  Spiral  Development  Cycle  of  Software[4] 


When  the  system  is  operating  normally,  the  equipment  color  will  be  changed  to  green  to 
indicate  a  component  is  operating.  An  example  of  this  is  shown  in  Figure  2-5  which  indicates 
that  a  Vacuum  Pump  has  started  and  is  opertating  normally.  Yellow  is  displayed  when  a 
fault  is  detected  but  for  an  insificient  number  of  samples  to  flag  it  as  a  fault.  This  threshold 
setting  prevents  spurious  alarms  in  the  system  but  allows  the  abnormal  indications  to  be 
logged  and  categorized  for  further  analysis. 

When  a  fault  in  the  system  is  received,  an  example  of  which  is  shown  in  Figure  2-6 
as  a  Level  Probe  Failure,  the  component  color  is  changed  to  red  and  a  audible  warning 
is  given.  To  access  the  detailed  information  screen,  the  operator  clicks  on  the  component 
using  the  pointing  device  or  the  electronic  pen  included  with  the  tablet  PC.  An  example  of 
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NILM  GUI 


Figure  2-3:  Data  Flow  Diagram  of  CHT  GUI  Program. 


the  detailed  information  screen  is  shown  in  Figure  2-7. 

The  failures  indicated  by  the  system  were  those  identified  through  previous  research 
on  the  CHT  system  [24,  32]  and  are  listed  in  Table  2.2.  The  Vacuum  Pump  Seal  Water 
Line  Clog  failure  indication  is  provided  by  the  classification  program  [34].  The  remaining 
failures  are  detected  by  routines  within  the  GUI  program.  The  failure  detection  algorythms 
are  contained  within  separate  modules  of  the  program.  This  allows  modifications  to  be 
made  quickly  arid  easily  if  the  algorithms  are  improved  or  require  adaptation  if  the  system 
is  installed  on  a  different  class  of  ship.  As  an  indication  of  the  effectiveness  of  this  approach, 
the  leak  detection  algorithms  were  developed  by  previous  research  [24,  32]  and  are  utilized 
by  this  program  by  a  standard  data  exchange  interface. 
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Figure  2-4:  CUT  GUI  at  Startup.  Shown  in  the  diagram  are  the  Collection  Tank,  the 
Vacuum  Gauge,  one  Vacuum  Pump  representing  both  physical  pumps  and  one  Discharge 
Pump  representing  both  physical  pumps. 


Table  2.2:  Failures  in  the  CHT  System  Indicated  on  the  GUI  Interface 

Vacuum  Tank  Level  Probe  Failure 
Vacuum  Pump  Seal  Water  Line  Clog 
Vacuum  Line  Clog 
Small  System  Leak 
Large  System  Leak 
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Ia*t24hr  VP  Runs 


•  G  X 


Non-lntrusive  Load  Monitoring 

Probe  Failure  Detected 


VP  Avt  Run  Tim* 


DP  Aye  Run  Time 

lest  24 hr  DP  Runs 


Figure  2-6:  CHT  GUI  with  Indication  of  Level  Probe  Failure 


06  31:35  EC0| 


Quit  NIL M 
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•  -  * 

NILM  Alarm  Information 

Probe  Failure  Detected 

CC315- 
Ol.kMulWb  Util  * 

0?.juwxik  i*.  ri  st 


Figure  2-7:  CHT  GUI  Alarm  Information  Screen.  The  yellow  trace  represents  the  Real 
Power  in  KW  and  the  red  trace  the  Reactive  Power  in  KVAR.  This  screen  shows  the  power 
traces  from  a  Discharge  Pump  start  and  stop.  The  Discharge  Pump  normal  run  time  is 
approximately  45  to  60  seconds.  These  traces  show  a  run  time  of  approximately  1  second 
and  is  an  indication  of  a  Vacuum  Tank  Level  Probe  failure. 
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Chapter  3 


Land-Based  Engineering  Site 
(LBES) 


3.1  Background 

The  LBES  facility  provides  an  opportunity  to  test  the  functionality  and  operation  of  the 
NILM  equipment  on  Navy  systems  without  impacting  the  operation  of  Navy  ships.  The 
LBES  facility  is  located  on  the  Naval  Business  Center  in  Philadelphia,  Pennsylvania.  The 
LBES  was  built  to  replicate  the  equipment  and  systems  of  the  #2  Engine  Room  of  the 
Navy’s  Arleigh  Burke  (DDG-51)  class  of  destroyers.  The  facility  provides  for  testing  by 
the  U.S.  Navy  and  for  contractors  working  on  shipboard  systems.  It  is  also  utilized  by  the 
crews  of  newly  built  destroyers  as  a  means  of  training  prior  to  operating  their  engineering 
equipment  for  the  first  time. 

The  equipment  installed  at  LBES  are  the  major  systems  found  in  #2  Engine  Room: 
two  LM2500  gas-turbines  propulsion  main  engines  (GTMs)  along  with  the  reduction  gears, 
shafting  and  bearings,  three  gas-turbine  generators  (GTGs),  and  their  associated  auxiliary 
systems  such  as  Fuel  Oil  (FO),  Lubricating  Oil  (LO),  Low  Pressure  Air  (LP  Air)  and  cooling 
water. 

Since  it  is  a  land-based  test  facility,  LBES  is  configured  differently  than  the  #2  Engine 
Room  on  a  DDG-51.  For  instance,  LBES  has  a  complete  Zonal  Electrical  Distribution 
System  (ZEDS)  with  all  three  GTGs  instead  of  the  single  GTG  installed  in  #2  Engine 
Room  on  DDG-51  class  ships.  Since  the  LBES  is  a  land-based  facility,  a  large  waterbrake 
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is  installed  to  simulate  the  Controllable  Pitch  Propeller  (CPP)  that  would  normally  be  on 
the  end  of  the  propulsion  shaft.  The  shaft  length  at  LBES  is  shorter  than  it  is  on  a  DDG-51 
class  ship,  but  the  waterbrake  appears  the  same  as  a  CPP  to  the  propulsion  train  [23,  3]. 
A  perspective  view  of  the  LBES  facility  is  shown  in  Figure  3-1  and  a.  picture  is  shown  in 
Figure  3-2 


Waterbrake 


Gas  Turbine  Generators 


Gas  Turbine 
Propulsion  Engines 


I  T  PI 


Gas  Turbine  Exhaust 


Switchboards 


Figure  3-1:  Layout  of  LBES  Facility  in  Philadelphia,  PA  [18] 


3.2  Monitored  Systems 

The  following  is  a  description  of  the  NILM  systems  installed  to  monitor  the  equipment 
at  the  LBES  facility.  Each  section  has  a  short  summary  of  the  equipment  monitored  and 
the  NILM  installation  parameters  needed  to  process  the  data.  A  diagram  of  a  portion  of 
the  electrical  distribution  system  used  to  power  engineering  loads  is  shown  in  Figure  3-3. 
The  location  of  the  NILM  sensors  is  indicated  in  the  figure  by  the  red  rectangles.  The 
LBES  facility  has  two  electrically  isolated  distribution  systems.  The  engineering  space 
equipment  is  powered  from  commercial  sources.  The  output  of  the  electrical  generators  in 
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Figure  3-2:  Picture  of  LBES  Facility  in  Philadelphia,  PA  [23].  The  waterbrake  is  in  the 
foreground.  The  propeller  shaft  can  be  seen  leading  from  the  waterbrake  to  the  shaft 
bearing,  which  is  located  just  aft  of  the  LM2500  Gas  Turbines. 

the  engineering  spaces  is  connected  to  a  bank  of  load  cells.  Therefore  no  NILM  systems  are 
installed  on  that  portion  of  the  electrical  distribution  system. 

3.2.1  Low  Pressure  Air  Compressor  #2 

The  LPAC  #2  NILM  system  was  installed  by  LCDR  Greene  [18].  The  NILM  for  LPAC 
#2  is  colocated  with  the  NILM  for  LPAC  #3  because  the  LPACs  are  physically  located 
adjacent  to  each  other.  Each  LPAC  provides  pressurized  air  for  the  use  in  pneumatic  valve 
control  and  pneumatic  power  for  various  equipment.  Only  one  LPAC  is  required  to  supply 
all  loads  at  the  LBES  facility.  In  the  event  both  LPACs  are  not  operational,  the  LP  Air 
system  is  supplied  via  a  site  air  system.  It  should  be  noted  that  there  is  no  LPAC  #1  at 
LBES  because  the  on  a  DDG-51  LPAC  #1  is  located  in  #1  Engine  Room. 

The  site  LP  Air  system  consists  of  two  large,  high  capacity  air  compressors,  air  receivers 
and  associated  valves  and  piping  which  supply  compressed  air  to  the  building  in  which  the 
LBES  facility  is  located.  An  isolation  valve  (FA-V007)  and  check  valve  (APL-V003)  are 
located  in  the  line  where  the  site  LP  Air  system  connects  to  the  LBES  LP  Air  system.  The 
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1NILM  Sensor  Location 

Unknown  Locations: 
-GTG  til 
-LOSP  2A  (3SA) 


Figure  3-3:  LBES  Electrical  Distribution  Diagram  shown  the  relevant  loads  that  are  moni¬ 
tored  by  NILM  systems. 


check  valve  prevents  the  LBES  LP  Air  system  from  supplying  air  to  the  site  system. 

The  NILM  voltage  tap  and  CT  sensor  are  installed  inside  the  controller  for  the  LPAC. 
The  details  of  the  NILM  configuration  are  shown  in  Table  3.1.  The  physical  installation  is 
shown  in  Figure  3-5  and  Figure  3-4. 


Table  3.1:  LPAC  #2  NILM  Data 


NILM 

Channel 

Measurement 

Phase 

Measurement 

Resistors 

Tranducers 

Gain 

Code 

1 

Voltage 

A  and  B 

18012 

LEM  LV-25P 

0 

2 

LPAC  Cur¬ 
rent 

c 

3612 

LEM  LA-305S 

0 

3 

Pressure 

None 

None 

Omega  PX209 

0 

The  LPAC  #2  NILM  is  set  to  monitor  4  channels  of  d 

ata  in  raw  data  mode. 
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Figure  3-4:  NILM  Intallation  on  LPAC  #2.  The  CT  is  connected  on  phase  A  and  the 
voltage  tap  is  connected  across  phases  B  &  C. 


3.2.2  Low  Pressure  Air  Compressor  #3 

The  Low  Pressure  Air  Compressor  (LPAC)  #3  NILM  was  installed  as  part  of  the  research 
for  this  thesis.  The  NILM  voltage  tap  and  CT  sensor  are  installed  inside  the  controller  for 
the  LPAC.  The  details  of  the  NILM  configuration  are  shown  in  Table  3.2.  The  physical 
installation  is  shown  in  Figure  3-5  and  Figure  3-6. 


Table  3.2:  LPAC  #3  NILM  Data 


NILM 

Channel 

Measurement 

Phase 

Measurement 

Resistors 

Tranducers 

Gain 

Code 

1 

Voltage 

B  and  C 

1300 

LEM  LV-25P 

0 

2 

Current 

A 

13.5S2 

LEM  LA-305S 

0 

The  LPAC  #3  NILM  is  set  to  monitor  2  channels  of  cl 

ata  in  raw  data  mode. 

3.2.3  Gas  Turbine  Generator  Controller 

Gas  Turbine  Generator  (GTG)  #1  was  chosen  for  monitoring  because  of  it  is  equiped  with  a 
Redundant  Independent  Mechanical  Start  System  (RIMSS).  The  RIMSS  consists  of  a  small 
DC  motor,  a  small  gas  producing  turbine  and  a  small  RIMSS  power  turbine  connected  to 
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Figure  3-5:  LPACs  #2  and  #3  at  LBES.  Both  NILM  Monitoring  systems  are  installed  on 
the  cart  in  the  foreground  of  the  picture. 


Figure  3-6:  NILM  Intallation  on  LPAC  #3.  The  CT  is  connected  on  phase  A  and  the 
voltage  tap  is  connected  across  phases  B  &  G. 
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the  reduction  gear  of  the  GTG.  The  GTG  is  designated  AG9140RF,  which  indicates  that 
it  is  an  Allison  501-K34  gas  turbine  and  a  KATO  3000  kW  electrical  generator.  GTG  #1 
can  operate  independent  of  the  ship’s  electrical  distribution  system  to  allow  it  to  be  started 
in  completely  dark  conditions.  To  accomplish  this,  it  is  equiped  with  a  No  Break  Power 
Supply  (NBPS)  which  is  essentially  an  uninteruptable  power  supply  for  the  Full  Authority 
Digital  Control  (FADC)  Local  Operating  Panel  (LOCOP).  The  FADC  LOCOP  controls  and 
monitors  the  gas  turbine,  the  RIMSS,  the  electrical  generator  and  all  associated  auxiliary 
systems.  The  NILM  is  connected  to  the  FADC  LOCOP  to  monitor  all  control  signals  to 
the  GTG. 


Table  3.3:  Gas  Turbine  Generator  NILM  Data 


NILM 

Channel 

Measurement 

Phase 

Measurement 

Resistors 

Reference 

Resistors 

Tranducers 

Gain 

Code 

1 

Voltage 

A  and  B 

110ft 

51ft 

LEM  LV-25P 

0 

2 

Current 

c 

147ft 

51ft 

LEM  LA-305S 

IT 

The  GTG  NILM  is  set  to  monitor  2  channels  of  data  in  preprocess  mode. 

Figure  3-7:  NILM  Intallation  on  the  GTG 
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3.2.4  Gas  Turbine  Propulsion  Main  Engine  Controller 


The  Gas  Turbine  Propulsion  Main  Engine  (GTM)  Controller  NILM  monitors  Power  Panel 
1-282-1  and  three  specific  loads  on  the  panel:  the  Universal  Engine  Controller  (UEC)  Plus 
for  the  2A  engine,  the  UEC  Plus  for  the  2B  engine,  and  the  the  Shaft  Control  Unit  (SCU). 
Power  Panel  1-282-1  provides  120 VAC  power  to  these  components  and  to  the  GTM  igniters 
and  the  GTM  module  light.  The  NILM  system  in  this  location  monitors  all  four  CTs  and 
the  voltage  tap  on  the  power  panel  which  is  common  to  all  channels.  Figure  3-8  shows 
the  power  panel  and  Figure  3-9  shows  the  CTs  within  the  power  panel.  The  details  of 
the  NILM  configuration  are  shown  in  Table  3.4,  and  details  on  each  monitored  system  (i.e. 
UEC,  SCU)  are  presented  in  separate  subsections  below. 


Figure  3-8:  NILM  Intallation  on  Panel  1-282-1  [3].  Several  loads  on  the  panel  are  annotated. 
The  SCU,  UEC  2A  and  UEC  2B  are  separately  monitored. 
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Figure  3-9:  NILM  CT  Connections  on  Panel  1-282-1  [3],  The  CTs  for  the  individual  equip¬ 
ment  are  on  the  left  of  the  panel  and  the  aggregate  panel  CT  is  located  in  the  upper  right. 


Table  3.4:  Panel  Aggregate  Current  NILM  Data  |3] 


NILM 

Channel 

Measurement 

Phase 

Measurement 

Resistors 

Tr  and  uce  rs 

Gain 

Code 

1 

Voltage 

B  and  C 

110  12 

LEM  LV-2.5P 

0 

2 

UEC  2A  Cur¬ 
rent 

A 

13012 

LEM  LA-55S 

3 

3 

UEC  2B  Cur¬ 
rent 

A 

13012 

LEM  LA-55S 

3 

4 

SCU  Current 

A 

13012 

LEM  LA-55S 

2 

5 

Panel  Current 

A 

6712 

LEM  LA-205S 

1 

The  GTN 

[  NILM  is  set  to  monitor  6  c 

lannels  of  data  in  raw  data  mode. 

Universal  Engine  Controller  Plus  2A  &  2B 


The  Universal  Engine  Controller  (UEC)  Pluses  provide  all  control  and  monitoring  func¬ 
tions  for  the  LM2500  gas  turbines  [29J.  Each  GTM  has  a  separate  UEC  that  operates 
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independently  of  the  other.  The  control  signals  required  to  start  and  operate  the  GTM  can 
be  observed  by  monitoring  the  power  to  the  UEC.  The  necessary  monitoring  and  analysis 
process  is  described  in  [3|.  Detailed  information  on  the  UECs  can  be  found  in  [31). 


Shaft  Control  Unit 

The  Shaft  Control  Unit  (SCU)  provides  a  local  control  interface  to  both  GTMs.  The  control 
functions  include  starting  and  stopping  the  gas  turbine,  fuel  valve  control,  Lube  Oil  Service 
pump  control,  and  Fuel  Oil  Service  Pump  Control.  The  SCU  is  normally  unmanned  with 
control  of  the  GTMs  at  the  Propulsion  and  Auxiliary  Control  Console  (PACC)  located  in 
the  LBES  Central  Control  Station  (CCS).  The  PACC  is  connected  to  the  SCU  via  fiber 
optic  cabling  [3].  The  SCU  is  monitored  by  the  NILM  to  enable  the  observation  of  the 
control  signals  for  the  control  of  the  GTMs  and  their  associated  auxiliary  equipment.  The 
necessary  monitoring  and  analysis  process  is  described  in  [3].  Detailed  information  about 
the  SCU  can  be  found  in  [28]. 


3.2.5  Fuel  Oil  Service  Pump  #2A 

The  Fuel  Oil  Service  Pumps  (FOSP)  supply  fuel  to  the  GTMs  directly  and  to  the  gravity 
head  tank  used  by  the  GTGs.  The  FOSPs  are  two-speed  pumps.  Only  one  pump  is 
required  for  normal  operation  of  the  engineroom  equipment.  2A  FOSP  is  monitored  by  a 
NILM  system.  The  details  of  the  NILM  configuration  are  shown  in  Table  3.5,  and  a  photo 
of  the  NILM  installation  is  shown  in  Figure3-10  and  Figure  3-11. 


Table  3.5:  Fuel  Oil  #2 A  NILM  Data 


NILM 

Channel 

Measurement 

Phase 

Measurement 

Resistors 

Tranducers 

Gain 

Code 

1 

Voltage 

A  and  B 

17712 

LEM  LV-25P 

0 

2 

2A  FOSP 

Current 

c 

6112 

LEM  LA-55S 

0 

The  FO  NILM  is  set  to  monitor  2  channels  of  data  in  preprocess  mode. 
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Figure  3-10:  NILM  Intallation  on  Fuel  Oil  Service  Pump  #2A  [18J. 


Figure  3-11:  NILM  Intallation  on  the  2A  FO  Pump 


3.2.6  Lube  Oil  Service  Pump  #2A 

The  Lubricating  Oil  Service  Pump  (LOSP)  supplies  lubricating  oil  to  the  Main  Reduction 
Gear  (MRG)  and  the  propulsion  shaft  thrust  bearing.  The  LOSP  has  two  operating  speeds. 
2A  LOSP  is  monitored  by  a  separate  NILM.  The  details  of  the  NILM  configuration  are 
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shown  in  Table  3.6  and  a  photo  of  the  NILM  installation  is  show  in  Figure  3-12  and  Figure  3- 
13. 


Figure  3-12:  NILM  Intallation  on  Lube  Oil  Pump  #2A.  [18] 


Table  3.6:  Lube  Oil  Pump  #2A  NILM  Data 


NILM 

Channel 

Measurement 

Phase 

Measurement 

Resistors 

Tranducers 

Gain 

Code 

1 

Voltage 

A  and  B 

i  ion 

LEM  LV-25P 

0 

2 

LOSP  Cur¬ 
rent 

c 

20ft 

LEM  LA-305S 

0 

The  LO  NILM  is  set  to  monitor  2  channels  of  data  in  preprocess  mode. 

3.2.7  ISA  &  3SA  Switchboards 

The  ISA  and  3SA  switchboards  are  located  on  the  upper  level  of  the  LBES  facility.  They 
are  supplied  by  a  public  utility  through  the  circuit  breakers  '’MAIN  BKR.  1SA-M1”  and 
"MAIN  BKR.  3SA-M1".  There  are  two  NILMs  installed,  one  on  each  switchboard,  with 
their  sensors  connected  on  the  upstream  side  of  these  circuit  breakers.  The  voltage  taps  are 
on  phases  A  &  B  arid  the  CTs  are  on  phase  C.  The  connections  to  panel  3SA  are  shown  in 
Figure  3-14.  The  connections  to  ISA  are  similar. 
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Figure  3-13:  NILM  Intallation  on  the  2A  LO  Pump 


The  ISA  and  3SA  NILMs  each  record  four  channels:  two  for  voltage  and  two  for  current. 
The  CTs  used  in  these  installations  are  different  than  the  ones  normally  used  in  NILM 
installations.  The  primary  reason  for  this  change  is  that  the  test  director  requested  the  use 
of  a  split-core  CT  in  order  to  avoid  breaking  electrical  connections.  A  Fluke  i3000S  AC 
Clamps  were  chosen  for  use  as  the  CTs.  The  specification  sheets  for  the  Fluke  clamps  are 
presented  in  Appendix  A. 

Both  switchboard  NILMs  use  Fluke  iSOOO’s  CTs  set  to  the  3000A  range.  At  this  setting, 
they  are  not  damaged  by  the  typical  operating  currents,  but  they  do  produce  very  small 
output  voltage  signals.  An  amplifer  circuit  was  designed  by  Professor  Cox  to  raise  the  signal 
to  a  more  reasonable  level.  The  amplifer  card  has  two  amplifiers  installed,  and  the  single 
input  from  the  CT  is  split  and  amplified  at  two  levels  in  the  current  configuration,  one 
at  lOx  and  one  at  100x.  Further  testing  will  determine  the  necessary  amplification  levels 
needed.  The  levels  can  be  easily  adjusted  by  changing  resistors  in  the  circuit.  The  amplifier 
card  is  shown  in  Figure  3-15. 

The  details  of  the  NILM  configuration  are  shown  in  Table  3.7  and  Table  3.8. 
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Table  3.7:  ISA  NILM  Configuration 


NILM 

Channel 

Measurement 

Phase 

Measurement 

Resistors 

Tr  and  uce  rs 

Gain 

Code 

1 

Voltage 

A  and  B 

non 

LEM  LV-25P 

1 

2 

Current 

c 

None 

Fluke  i3000S 

1 

3 

Voltage 

A  and  B 

non 

LEM  LV-25P 

1 

4 

Current 

c 

None 

Fluke  i3000S 

4 

The  ISA 

NILM  is  set  to  monitor  4  c 

lannels  of  data  in  pre process  mode. 

Table  3.8:  3SA  NILM  Configuration 


NILM 

Channel 

Measurement 

Phase 

Measurement 

Resistors 

Tranducers 

Gain 

Code 

1 

Voltage 

A  and  B 

non 

LEM  LV-25P 

4 

2 

Current 

c 

None 

Fluke  i3000S 

1 

3 

Voltage 

A  and  B 

lion 

LEM  LV-25P 

4 

4 

Current 

c 

None 

Fluke  i3000S 

4 

The  3SA  NILM  is  set  to  monitor  4  c 

lannels  of  data  in  preprocess  mode. 

3.3  Data  Analysis 

During  the  course  of  this  thesis  project,  two  tests  were  performed  on  the  LBES  systems.  The 
first  test  was  to  simultaneously  record  power  signals  at  both  the  auxiliary  equipment  and 
at  the  switchboards.  The  second  was  to  determine  the  cause  of  abnormal  signals  detected 
in  the  LP  Air  system  by  previous  researchers  [3,  23]. 

3.3.1  Power  Signal  Test  of  Auxiliary  Equipment 

The  NILM  systems  from  all  previous  research  were  installed  locally  at  each  load  or  on  a 
power  panel  near  the  circuit  breaker  for  the  load.  To  test  the  effectiveness  of  installing 
NILMs  at  the  main  distribution  switchboards  each  load  was  started  and  operated  in  all 
available  modes  while  the  local  NILM  and  the  main  switchboard  NILMs  were  monitoring 
the  system.  Neither  the  GTMs  nor  the  GTGs  were  available  for  operation  during  the 
testing  period.  Therefore,  for  the  test  of  the  3SA  NILM,  LPAC  #2  was  started  first  in 
automatic  mode  to  provide  some  variation  in  the  signal  as  the  LPAC  cycled  during  it’s 
normal  operation.  The  other  loads  known  to  be  powered  from  this  switchboard  were  the 
2A  LOSP  and  2A  FOSP.  The  LOSP  is  the  larger  of  the  two  motors  and  was  chosen  to  be 
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Figure  3-14:  NILM  Intallation  on  the  3SA  Switchboard.  The  connections  for  the  ISA 
Swithboard  are  similar.  The  voltage  tap  wires  can  be  seen  on  Phases  A  &  B.  The  Fluke 
CT  can  be  seen  clamped  around  Phase  C. 


Figure  3-15:  Amplifier  card  used  to  amplify  the  signal  from  the  Fluke  CT  for  use  within 
the  NILM. 
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Figure  3-16:  Photograph  of  the  NILM  location  for  the  ISA  &  3SA  Switchboards. 


started  first.  It  was  started  in  slow  speed,  ran  for  approximately  one  minute  then  shifted 
to  fast  speed  without  stopping.  The  FOSP  was  then  started  in  slow  speed  and  shifted  to 
fast  speed  approximately  one  minute  later.  It  was  opperated  at  fast  speed  then  shifted 
back  to  slow  speed  prior  to  securing.  Once  the  FOSP  was  secured,  the  LOSP  was  secured. 
Figure  3-17  shows  the  results  of  this  test.  Table  3.9  gives  an  explanation  of  the  events 
annotated  in  the  figure.  The  results  of  the  test  provided  the  expected  indications  and  all 
signals  were  detectable  on  both  the  local  NILMs  and  on  the  3SA  NILM. 

The  second  portion  of  the  test  was  to  perform  a  similar  test  for  the  loads  powered  from 
the  ISA  switchboard.  LPAC  #3  was  started  first  and  placed  in  automatic  mode  to  provide 
a  cycling  load  while  the  loads  on  the  switchboard  were  operated.  The  loads  thought  to  be 
powered  from  ISA  were  the  B'  series  of  pumps,  2B  LOSP  and  2B  FOSP.  These  pumps  do 
not  have  local  NILMs  installed  to  monitor  them,  therefore,  only  the  ISA  NILM  was  expected 
to  record  their  operating  signals.  The  3SA  NILM  was  also  operating  during  this  test  but 
was  not  expected  to  gather  any  useful  data.  The  LOSP  and  the  FOSP  were  operated  in 
the  same  order  and  manner  as  during  the  first  test.  It  was  noted  while  processing  the  data 
that  the  ’B’  series  of  pumps  are  indeed  powered  off  of  the  3SA  switchboard.  The  only  load 
monitored  on  the  ISA  switchboard  during  the  test  was  LPAC  #3.  The  LOSP  and  the  FOSP 
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were  monitored  by  the  3SA  NILM  and  as  shown  in  Figure  3-18  were  operated  without  any 
cycling  load.  Table  3.9  gives  an  explanation  of  the  events  annotated  in  Figure  3-18. 

The  NILM  monitors  installed  on  the  ISA  and  3SA  switchboards  are  able  to  monitor 
other  equipment  that  has  not  been  monitored  previously  at  the  LBES  facility.  An  example 
of  the  power  trace  from  one  piece  of  equipment  is  the  FO  Purifier  shown  in  Figure  3-19. 
Table  3.11  gives  an  explanation  of  the  events  annotated  in  the  figure.  The  interesting  shape 
of  the  power  curve  is  due  to  the  eight-step  starting  sequence  of  the  FO  Purifier,  which  is 
necessary  due  to  the  rotating  mass  of  the  bowl.  The  bowl  is  the  internal  assembly  that 
rotates  at  a  high  rate  of  speed  and  is  used  to  remove  impurties  from  the  FO.  An  interesting 
event  occurs  in  the  shutdown  sequence  that  is  detectable  in  the  power  trace.  The  operator 
'’shoots  the  bowl”  with  water  to  provide  some  cleansing  of  the  bowl  during  shutdown1.  This 
step)  is  peformed  twice  per  procedure.  This  process  can  be  detected  by  the  NILM  and  is 
shown  in  Figure  3-19. 

Another  aspect  of  observing  the  aggregate  signal  is  that  the  signals  of  interest  will  be 
contained  in  background  noise  of  other  loads  at  LBES.  An  example  of  this  is  shown  in 
Figure  3-20  which  shows  the  effects  of  the  Lube  Oil  (LO)  Heaters  operating  and  LPAC  #2 
is  cycling  simultaneously.  Table  3.12  gives  an  explanation  of  the  events  annotated  in  the 
figure.  The  LO  Heaters  have  this  type  of  signal  because  of  the  power  modulation  of  the 
signal  to  give  a  relatively  high  average  power  output  of  the  heaters  without  damaging  the 
heaters  due  to  excessive  current.  The  cycling  frequency  of  the  signal  is  shown  in  Figure  3-21. 
The  annotation  in  the  figure  is  the  same  as  in  Figure  3-20.  The  cycling  of  LPAC  #2  can  be 
seen  in  the  overall  signal  at  the  lower  edges  of  the  power  levels  of  the  cycling  heaters.  This 
can  be  effective  when  observing  only  two  loads,  however,  when  multiple  loads  are  operating 
the  individual  signals  will  not  be  able  to  be  extracted  easily.  A  potential  solution  for  this 
problem  is  the  ginzu  program  detailed  in  [34]. 


lrThe  theory  of  operation  of  the  FO  Purifier  was  obtained  through  discussion  with  Frank  Facciolo  at  the 
LBES  facility. 
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Figure  3-17:  Plot  of  Power  on  3SA  Switchboard  during  Power  Test  Procedure.  The  events 
are  annotated  with  letters  and  are  explained  in  Table  3.9 


Table  3.9:  Explanation  of  events  in  Figure  3-17 


Event 

Description 

A 

LPAC  #2  Started  in  Automatic 

B 

LPAC  #2  cycling 

C 

2A  LOSP  Started  in  Slow  Speed 

D 

2A  LOSP  Shifted  to  Fast  Speed 

E 

2A  FOSP  Started  in  Slow  Speed 

F 

LPAC  #2  cycling 

G 

2A  FOSP  Shifted  to  Fast  Speed 

H 

2A  FOSP  Shifted  to  Slow  Speed 

I 

2A  FOSP  Secured  &  LPAC  #2  cycling 

J 

2 A  LOSP  Secured 

K 

LPAC  #2  cycling 
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Figure  3-18:  Plot  of  Power  on  3SA  Switchboard  during  Power  Test  Procedure  Showing  TT 
Pumps.  The  events  are  annotated  with  letters  and  are  explained  in  Table  3.10 


Table  3.10:  Explanation  of  events  in  Figure  3-18 


Event 

Description 

A 

2B  LOSP  Started  in  Slow  Speed 

B 

2B  LOSP  Shifted  to  Fast  Speed 

C 

2B  FOSP  Started  in  Slow  Speed 

D 

2B  FOSP  Shifted  to  Fast  Speed 

E 

2B  FOSP  Shifted  to  Slow  Speed 

F 

2B  FOSP  Secured 

G 

2B  LOSP  Secured 

H 

Waste  Pump  Started 

I 

!  Waste  Pump  Secured 
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Figure  3-19:  Plot  of  Power  for  the  Start  of  the  Fuel  Oil  Purifier.  The  events  are  annotated 
with  letters  and  are  explained  in  Table  3.11 


Table  3.11:  Events  for  the  Start  of  the  Fuel  Oil  Purifier  Corresponding  to  Figure  3-19 


Event 

Description 

A 

FO  Purifier  Start  Sequence  Commenced 

B 

FO  Start  Sequence  Complete,  FO  Purifier  FO  Pump  Started 

C 

FO  Purifier  Online 

D 

FO  Purifier  Offline 

E 

’’Shooting  the  Bowl” 

F 

FO  Purifier  Secured 
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Time  (Mm) 


Figure  3-20:  Plot  of  Power  for  Cycling  of  LPAC  #2  with  Lube  Oil  Heaters  Energized.  The 
events  are  annotated  with  letters  and  are  explained  in  Table  3.12 


Table  3.12:  Events  for  LP  Air  Test  Corresponding  to  Figure  3-20 


Event 

Description 

A 

Powrer  due  to  cycling  of  LO  Heaters 

B 

Effect  of  LPAC  #2  Cycling  visible. 
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Figure  3-21:  Close  up  of  Power  Trace  for  Cycling  of  Lube  Oil  Heaters.  'A'  in  Figure  is  the 
same  as  Figure  3-20 


3.3.2  LP  Air  Tests 

The  power  signals  on  LPAC  #2  were  noted  to  be  abnormal  by  LCDR  Bennett  during  his 
research  [3].  Figure  3-22  shows  the  LPAC  #2  power  trace  recorded  by  LCDR  Greene  in 
April  2005  and  it  is  evident  that  the  loading  and  unloading  times  of  the  compressor  occur 
with  a  regular  frequency  during  times  of  normal  load.  When  a  large  load  is  placed  on  the 
system,  such  as  when  a  GTM  starts,  LCDR  Greene  noted  that  the  power  trace  indicated  a 
variation  in  its  normal  pattern.  By  comparison  Figure  3-23  shows  the  power  trace  of  LPAC 
#2  in  April  2007.  The  unloaded  states  of  the  compressor  are  no  longer  regular  and  show  an 
increase  in  the  unloaded  times  at  irregular  intervals.  Figure  3-24  shows  an  enlarged  portion 
of  the  power  trace  for  LPAC  #2  in  April  2007  along  with  the  LP  Air  Receiver  pressure. 
This  plot  shows  that  the  compressor  unload  time  and  the  LP  Air  Receiver  pressure  are 
correlated,  the  plot  also  shows  unusual  behavior  during  certain  unload  intervals.  The  cause 
of  the  extended  unload  periods  was  unknown  to  LCDR  Bennett  but  they  were  not  found 
to  relate  to  any  plant  operations. 


50 


Figure  3-22:  Plot  of  Power  for  LPACs  #2  from  April  2005.  The  load  and  unload  times  of 
the  compressor  are  at  a  regular  frequency  during  times  of  normal  loading.  [3] 


Figure  3-23:  Plot  of  Power  for  LPACs  #2  from  April  2007.  The  unloaded  times  of  the 
compressor  are  not  regular  and  increase  at  irregular  intervals.  [3] 


An  LP  Air  test  was  conducted  by  operating  #2  LPAC  in  automatic  mode,  first  with  the 
site  air  supply  valve  open  then  with  it  shut.  The  operation  of  the  LPAC  was  observed  during 
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Figure  3-24:  Close-up  Plot  of  Power  for  LPACs  #2  from  April  2007.  The  unload  times  of 
the  compressor  are  correlated  to  the  LP  Air  System  pressure.  [3] 


the  test  and  the  NILM  data  was  analyzed  upon  completion.  LPAC  #2  NILM  monitors  a 
pressure  transducer  installed  on  the  gauge  line  to  the  LP  Air  Receiver.  The  relative  locations 
of  components  withing  the  LP  Air  system  are  shown  on  the  LP  Air  System  Diagram  for 
LBES  in  Appendix  B.  The  test  was  then  repeated  on  LPAC  #3  to  confirm  the  results  and 
to  observe  the  cycling  behavior  of  each  LPAC.  The  LP  Air  Test  procedure  is  included  in 
Appendix  C. 

Figure  3-25  shows  the  results  of  the  test  on  LPAC  #2.  The  details  of  the  events  an¬ 
notated  in  the  figure  are  shown  in  Table  3.13.  The  upper  trace  shows  the  power  signal  to 
LPAC  #2.  A  25  SCFM  leak  was  introduced  into  the  system  to  provide  an  air  load  because 
normal  plant  loads  were  not  operating  during  the  test.  Note  that  the  LPAC  exhibits  the 
behavior  with  the  excessive  unload  time  when  FA-V007  is  open.  The  bottom  trace  shows 
the  LP  Air  System  pressure  decreases  to  approximately  111  PSIG  around  13  minutes  into 
the  test  then  increases  to  approximately  122  PSIG  .  It  can  be  noted  that  the  increase  in 
pressure  is  not  caused  by  the  LPAC  #2  as  the  power  signal  did  not  change  during  this  time 
and  the  LPAC  remained  in  its  unloaded  state.  At  about  15  mintues  into  the  test  the  Site 
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Air  Isolation  Valve,  FA-V007,  was  shut.  The  LPAC  #2  power  trace  and  the  pressure  trace 
indicate  that  the  pressure  cycles  and  LPAC  load  and  unload  states  are  now  correllated.  The 
LP  Air  System  pressure  increases  during  the  LPAC  load  states  and  decreases  during  LPAC 
unload  states.  Approximately  28  minutes  into  the  test  FA-V007  was  reopened  and  the  LP 
Air  System  pressure  returned  to  previous  cycling  and  LPAC  #2  returned  to  very  short  load 
states  and  excessively  long  unload  states. 

A  second  set  of  tests  were  conducted  using  LPAC  #3.  These  tests  also  involved  actuation 
of  FA-V007.  As  shown  in  Figure  3-26  similar  results  were  obtained.  The  details  of  the 
events  annotated  in  the  figure  are  shown  in  3.14.  It  is  easy  to  see  that  LPAC  #3  has 
a  different  operating  power  profile  than  #2  however,  the  same  extended  unload  state  is 
evident.  Approximately  32  minutes  into  the  test  the  operating  pressure  of  LPAC  #3  was 
changed  to  110  to  120  PSIG  to  observe  how  the  signals  would  change  as  monitored  by  the 
NILM  system. 

Conversations  with  an  LBES  technician  indicated  that  the  site  air  compressors  had 
been  replaced  with  larger,  higher  capacity  compressors.  Upon  investigation,  the  cycling 
pressure  range  of  the  site  air  compressor  was  found  to  be  115  to  125  PSIG  at  the  discharge 
of  the  compressor.  The  expected  discharge  pressure  for  the  site  air  compressors  is  90  PSIG. 
This  would  provide  approximately  80  PSIG  at  the  point  where  the  site  LP  Air  system 
taps  into  the  LBES  LP  Air  system,  which  is  lower  than  the  minimum  cycling  setpoint  for 
both  LPACs.  With  the  site  air  compressors’  operating  band  set  to  115  to  125  PSIG,  the 
pressure  on  the  site  side  of  check  valve  APL-V003  will  exceed  the  pressure  on  the  LBES  side 
during  normal  cycling  of  the  LBES  LPACs.  Therefore,  the  site  air  will  supply  the  LBES  LP 
Air  system.  The  LBES  LPACs  will  sense  the  pressure  being  maintained  above  their  lower 
cycling  setpoint  and  as  a  result  will  have  long  unload  states.  The  LBES  LPACs  will  not 
cycle  on  LP  Air  load  as  designed.  There  is  no  negative  effect  on  the  operation  of  LBES  as 
a  test  facility  or  as  a  training  environment  for  the  ship’s  crews. 


Figure  3-25:  Plot  of  Power  for  LPACs  #2  and  #3  with  LP  Air  Header  Pressure.  Upper 
trace:  LPAC  #2  power.  Middle  trace:  LPAC  #3  power,  Bottom  Trace:  Receiver  pressure. 
Label  locations  correspond  to  changes  in  system  conditions.  The  events  are  annotated  with 
letters  and  are  explained  in  Table  3.13 


Table  3.13:  Events  for  LP  Air  Test  Corresponding  to  Figure  3-25 


Event 

Description 

A 

Site  LP  Air  Isolation  valve  FA-V007  shut 

B 

Site  LP  Air  Isolation  valve  FA-V007  open 
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Figure  3-26:  Plot  of  Power  for  LPACs  #2  and  #3  with  LP  Air  Header  Pressure.  Upper 
trace:  LPAC  #2  power.  Middle  trace:  LPAC  #3  power,  Bottom  Trace:  Receiver  pressure. 
Label  locations  correspond  to  changes  in  system  conditions.  The  events  are  annotated  with 
letters  and  are  explained  in  Table  3.14 


Table  3.14:  Events  for  LP  Air  Test  Corresponding  to  Figure  3-26 


Event 

Description 

A 

Site  LP  Air  Isolation  valve  FA-V007  open 

B 

Site  LP  Air  Isolation  valve  FA-V007  shut 

C 

Output  Pressure  of  #3  LPAC  changed  to  110-120  PSIG 

Figure  3-27:  Plot  of  Power  for  a  Decreased  Air  Load  on  LPAC  #3.  The  events  are  annotated 
with  letters  and  are  explained  in  Table  3.15 


Table  3.15:  Events  for  LP  Air  Test  Corresponding  to  Figure  3-27 


Event 

Description 

A 

LPAC  #3  on  in  automatic  at  115  to  125  PSIG 

B 

Decreased  air  load  on  LP  Air  System 

C 

LPAC  #3  Secured 

Chapter  4 


Multi-Function  Monitor  (MFM) 

4.1  Introduction 

Modern  surface  warships  employ  some  form  of  Zonal  Electrical  Distribution  System  (ZEDS) 
to  increase  the  reliability,  quality  of  service  [15]  and  survivability  or  the  ship  in  the  event  the 
system  is  damaged  [16,  33].  The  U.S.  Navy  has  plans  to  develop  and  deploy  an  Integrated 
Power  System  (IPS)  for  future  surface  warships.  IPS  systems  combine  the  electrical  power 
generation  and  propulsion  power  systems  [13].  The  MFM  is  a  central  component  of  these 
systems  to  provide  switching  and  fault  isolation  [6]. 

4.2  MFM  System  Description 

In  naval  shipboard  electrical  distribution  systems  prior  to  utilizing  a  ZEDS,  fault  isolation 
was  performed  by  selecting  the  overcurrent  settings  on  circuit  breakers  such  that  those 
closest  to  the  fault  will  trip  first.  This  process  is  called  selective  tripping  [27].  If  the  system 
works  as  designed,  the  fault  will  be  isolated  while  minimizing  the  impact  on  the  rest  of 
the  system.  Various  studies,  analyses  and  tests  have  shown  that  selective  tripping  may  not 
adequately  protect  the  electrical  distribution  system  during  battle  damage  if  the  damage 
is  not  a  single  point  failure  [33].  With  the  implementation  of  ZEDS,  the  situation  was 
further  complicated  because  the  direction  of  current  flow  through  the  components  in  the 
system  depends  on  the  system  configuration.  This  necessitated  a  protective  system  which 
could  sense  the  direction  of  current  flow  and  provide  intelligent  fault  detect  ion  and  isolation 
actions  [33].  The  ZEDS  system  for  a  DDG-51  Flight  IIA  Arleigh  Burke  class  destroyer  is 
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shown  in  Figure  4-1. 


- >  ^ - *  Charnel  1  PTs  - ^Chamd  1  CTs 

Positive  Direction  Shunt  Trip  Signal  *  Charnel  2  PTs  Charnel  2  CTs 

of  Current  Flow 


Figure  4-1:  DDG-51  Flight  IIA  ZEDS  [6]  This  diagram  shows  the  addressing,  locations 
and  signal  inputs  of  the  MFMs  located  in  the  system.  Also  shown  is  the  defined  positive 
direction  of  current  flow  for  each  MFM. 

The  MFM-I  was  designed  to  accomplish  this  task.  The  MFM- III  is  the  improved  version 
of  the  MFM-I  and  is  discussed  in  this  thesis.  Each  MFM-III  has  two  sets  of  sensors,  one 
located  on  each  side  of  its  associated  circuit  breaker.  The  exception  to  this  are  the  generator 
output  breaker  which  utilize  a  single  set  of  sensors.  There  are  two  line  voltage  measurements: 
vab  and  Vbc  from  two  potential  transformers  (440:110),  and  three  phase  currents:  ia,H  and 
ic  from  current  transformers  (CTs)  (6000:5).  The  voltage  and  current  inputs  are  used  to 
calculate  the  power  flowing  through  the  circuit  breaker.  The  CTs  provide  an  indication  of 
the  direction  of  current  flow.  The  trip  status  of  the  associated  circuit  breaker  is  also  an 
input  to  the  MFM  logic  circuit  [6].  A  functional  diagram  of  the  MFM  is  shown  in  Figure  4-2. 

There  are  four  types  of  MFMs  based  on  where  they  are  to  be  placed  in  the  distribution 
system.  In  order  for  the  software  to  be  able  to  determine  the  correct  system  status  and 
fault  determination,  each  MFM  is  assigned  a  unique  identifier  which  indicates  the  specific 
location  in  the  system  where  it  is  located.  To  provide  a  means  for  the  MFMs  to  pass  data 


to  the  other  MFMs  in  the  system,  each  MFM  is  equipped  with  three  ethernet  ports.  Port 
#1  is  used  to  transmit  and  receive  data  on  the  ring  connection  with  all  other  MFMs.  Data 
is  transmitted  through  port  #1  every  5.0  ms.  Ports  #2  &  #3  are  used  to  transmit  and 
receive  data  point-to-point  with  the  MFMs  on  either  side.  Data  is  transmitted  through 
these  ports  every  1.0  ms.  This  ability  of  the  MFMs  to  exchange  data  throughout  the  entire 
system  allows  for  a  coordinated  response  to  a  detected  fault.  The  system  is  designed  to 
initiate  a  trip  response  within  10  ms  of  detection  (6J. 


Figure  4-2:  Functional  Diagram  of  MFM  [6]  The  diagram  shows  the  sensor  layout  for  the 
MFM.  The  three  phase  current  and  two  voltage  measurements  per  channel,  the  three  eth¬ 
ernet  communications  channels  and  the  shunt  trip  output  to  the  associated  circuit  breaker. 
The  shunt  trip  status  input  of  the  circuit  breaker  is  not  shown  in  the  diagram. 


4.3  MFM  System  Theory  of  Operation 

There  are  two  main  software  routines  that  the  MFM-III  uses.  The  first  is  called  the  High 
Speed  Relay  (HSR)  and  was  developed  by  Barrons  Associates  Inc  (BAI).  The  HSR  routine 
is  used  to  determine  if  certain  faults  exist  and  to  calculate  power  and  current  levels.  These 
values  are  inputs  to  the  second  routine,  the  Integrated  Protective  Coordination  System 
(IPCS).  IPCS  was  developed  by  the  Naval  Surface  Warfare  Center,  Carderock  Division 
Code  982  with  assistance  from  Micheal  Sieleman  at  Bath  Iron  Works  (BIW)  [6].  The  inputs 
to  the  HSR  routine  are  similar  to  those  used  by  NILM.  Therefore,  the  HSR  routine  will  be 
looked  at  in  detail. 


4.3.1  The  HSR  Algorithm 

The  HSR  routine  contains  three  main  routines:  the  Park’s  transformation,  a  fault  detection 
algorithm,  and  fault  direction  determination.  There  are  five  outputs  of  the  HSR  Algorithm 
that  are  used  by  the  IPCS  Algorithm:  Angle  and  Magnitude  (from  Park's  Transformation), 
the  change  in  the  Angle,  Instantaneous  Power,  and  Fault  Direction.  The  details  of  these 
routines  and  outputs  are  discussed  below  in  this  section. 


Park’s  Transformation 

A  Park's  transformation  is  used  to  change  the  reference  frame  from  that  of  the  physical 
stator  to  that  of  the  theoretical  rotating  3-phase  electrical  signal  .  The  general  Park's 
Transformation  matrix  that  is  used  in  the  HSR  algorithm  is 
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where  /  can  represent  voltage  or  current  [20].  The  Park's  Transformation  is  normally  per¬ 
formed  using  all  three  values  for  /  but  the  MFM  only  measures  two  line  voltages,  therefore 
the  equations  must  be  modified  to  reflect  the  measurements  taken.  A  useful  relation  of  a 
3-phase  electrical  circuit  is  that  the  sum  of  the  line  voltages  must  be  equal  to  zero, 


Vab  +  Vbc  T  Vca  ~~  0* 


(4.2) 


Looking  at  the  first  row  of  the  Park's  Transformation  to  compute  Fci  using  the  line  voltages 
for  /  yields 

y,i  =  VabCOsO  +  Vbccos  (0  -  ^)  +  Vcacos  (0  +  ^)  .  (4.3) 

Using  the  relationship  of  Eq.  (4.2)  and  substituting  for  Vca  gives 

Vd  =  Vabcos0  +  Vh,  cos  (0  -  ^)  -  (Vab  +  Vbc)  cos  (0  +  ^)  .  (4.4) 


Multiplying  through  and  gathering  like  terms 


Vd  =  Vn 
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(4.5) 


60 


Using  the  trigonometric  identity 


co.su  —  cosv 


—2sin 


u  +  v 
2 


(4.6) 


and  substituting  into  Eq.  (4.5) 
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Simplifying  gives 


Vd  =  2Vab 


sin  (o  +  0 


+  214, 


sinOsin 


(4.8) 


Similarly,  for  the  second  row  of  the  Park's  Transformation  to  compute  Fq  using  the  line 
voltages  for  /  yields 


Vq  =■  VabsinO  + 


0  — 


2n 


4"  Vcasin  (  0  -|- 


Using  the  relationship  of  Eq.  4.2  and  substititing  for  Vca  gives 

Vi  =  VabsinO  +  Vbcsin  (o  -  -  (K.&  +  Vbc )  «in  ((>  +  y )  . 

Multiplying  through  and  gathering  like  terms 
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(4.10) 


(4.11) 


Using  the  trigonometric  identity 


sinu  —  shiv 


2  cos 


(4.12) 


and  substituting  into  Eq.  (4.11) 


V,,  =  V, 


ab 


2  COS 
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0-0 -?f 
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(4.13) 


Simplifying  gives 


Vi  =  -2Vab 


cos 


2Vfcc 


cosOsin 


(4.14) 


Noting  that  2 sin  (|r) 
follows 


=  2 sin  Eq.  (4.8)  and  Eq.  (4.14)  can  be  further  simplifed  as 


Vd  =  2  sin 


VbcsinO 


(4.15) 


and 


Vq  =  —2  sin 


COS 


+  Vbccos9 


(4.16) 


The  Fq  value  is  not  calculated  in  the  HSR  routine.  This  value  checks  the  three  voltage  inputs 
to  verify  that  they  sum  to  zero.  Since  only  two  line  voltages  are  measured,  the  system  is 
assumed  to  be  a  balanced  3- phase  system  and  it  is  assumed  that  the  sum  of  the  voltages 
equals  zero  [18].  Once  the  Fq  and  Fd  values  are  calculated,  they  are  used  to  compute  the 
magnitude  and  angle  of  the  phasor 


Magnitude  =  -z\J  F<f  F^ 

(4.17) 

/  F.  \ 

Angle  =  arctan  yjjr J  * 

(4.18) 

Fault  Detection 

The  Fault  Detection  routine  uses  two  different  tests  to  determine  if  a  fault  condition  exists, 
a  large  change  in  the  current  angle  computed  by  the  Park’s  Transformation  from  the  average 
angle  from  past  eight  samples  and  comparing  the  magnitude  from  the  Park’s  Transformation 
to  high  and  low  thresholds. 

Angle  Fault  The  A  Angle  calculated  as  the  change  in  the  angle  computed  in  the  Park’s 
Transformation  from  the  average  of  the  past  eight  valid  samples.  If  a  large  change  in  the 
angle  is  detected,  it  is  an  indication  of  a  severe  increase  in  power  and/or  a  reversal  of  the 
power  flow.  Both  of  these  conditions  are  indicators  that  a  fault  condition  may  exist.  If 
the  calculated  angle  and  the  average  angle  differ  in  sign,  the  angle  buffer  is  unwrapped  by 
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adding  7r  to  the  angles  that  are  in  the  buffer  so  that  all  of  the  angles  are  in  the  quadrants  I 
and  IV.  The  absolute  value  of  the  A  Angle  is  compared  to  a  threshold  value.  If  it  is  exceeded 
the  fault  flag  is  set  to  ’’Angle  Fault.” 

Magnitude  Fault  The  newly  calculated  magnitude  is  compared  to  low  and  high  thresh¬ 
olds  to  verify  if  it  is  within  the  acceptable  band  of  what  is  expected.  If  it  is  outside  this 
band,  then  it  is  an  indication  that  a  severe  increase  in  power  flow,  which  is  an  indication 
that  a  fault  condition  exists.  The  fault  flag  is  set  to  ’’Magnitude  Fault.” 

Fault  Direction 

If  the  Fault  Detection  algorithm  indicates  that  a  fault  exists,  the  power  is  calculated  through 
both  input  channels  to  determine  if  the  fault  is  upstream  or  downstream  of  the  MFM.  The 
CTs  are  oriented  specifically  during  installation  such  that  positive  power  values  indicate 
current  is  flowing  from  upstream  to  downstream  of  the  MFM.  Therefore  the  power  calcula¬ 
tion  will  provide  an  indication  of  the  direction  of  the  electrical  energy  flowing  through  the 
circuit  breaker. 


Calculate  Average  Power  By  definition  the  instantaneous  power  is 

P  =  Vaia  +  Vbh  +  Vcic.  (4.19) 

The  MFM  does  not  monitor  va ,  or  vc  but  rather  vab  and  Therefore  the  power  must 
be  computed  in  terms  of  the  actual  measured  voltages.  This  can  be  accomplished  by  noting 
that  the  relationship  between  phase  voltages  va  and  and  line  voltage  vab  is 


vab  =  Va  -  Vb. 


(4.20) 


Similarly  for  Vbc  and  vca 


Vbc  =  Vb  -  vc 


and 


Vca  ~  Vc  -  Va 


(4.21) 


(4.22) 
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Substituting  Eq.  (4/20)  and  Eq.  (4.21)  into  Eq.  (4.19)  yields 


P  -  («ab  +  Vb)  ia  f  (>-’bc  +  Vc)  ib  +  Vcic.  (4.23) 

Noting  that  the  current  for  a  3-phase  electrical  system  must  satisfy  the  relation 

ia  4-  ib  +  i-c  —  0  (4.24) 

and  solving  for  ib ,  this  can  be  substituted  into  Eq.  (4.23)  resulting  in 

p  =  (vab  +  vbc  +  Vc)  ia  +  (*' be  +  vc)  (~ia  ~  ic)  +  vcic .  (4.25) 

Gathering  like  terms 

P  —  Vahid  4-  Vbcia  +  vcia  ~  vbcia  -  vbcic  -  Vcia  -  Vcic  +  vcic.  (4.26) 

Simplifying  this  equation  gives 

P  =  vabia  -  Vbcic •  (4.27) 

This  gives  an  accurate  calculat  ion  of  power  in  the  circuit  assuming  that  the  loads  are  equally 
balanced.  If  there  is  a  fault  on  phase  B,  the  power  would  not  be  accurate  since  ib  is  not 
included  in  the  equation.  To  eliminate  uncertainties  and  provide  an  indication  of  a  single 
phase  fault  a  method  must  be  used  such  that  all  measurements  receive  equal  weight.  Noting 
that  Eq.  (4.27)  was  arrived  at  by  substituting  for  ib  in  Eq.  (4.23),  two  more  equations  can 
be  found  by  substituting  for  ia  and  ic.  Solving  for  ia  in  Eq  (4.24)  and  substituting  into 
Eq  (4.23)  results  in 


P  —  ( vab  +  vbc  T  vc)  (  ib  ~~  ic)  4-  {vbc  4~  vc)  ib  4~  vcic.  (4.28) 


Simplifying  the  equation  results  in 

P  =  -vabib  -  vabic  -  vbcic.  (4.29) 
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Gathering  like  terms  results  in 


P  =  -Vab  (4  +  ic)  ~  vhcic.  (4.30) 

Similarly,  solving  Eq  (4.24)  for  ic  and  substituting  in  to  Eq  (4.23)  results  in 

p  =  Vabia  +  Vbcia  +  vja  +  vbcib  +  vcib  -  vcia  -  vcib.  (4.31) 

Simplifying  the  equation  results  in 

p  =  Vab  ia  +  vbcia  +  Vbcib.  (4.32) 

Gathering  like  terms  results  in 

P  ~  vabia  +  vbc  (ia  +  ib)  •  (4.33) 

Each  of  these  three  equations  for  power  compute  the  power  for  the  entire  circuit.  To  arrive  at 
an  accurate  power  computation  where  all  three  phase  currents  are  accounted  for  equally,  the 
average  of  Equations  (4.27),  (4.30),  and  (4.33)  is  calculated.  Therefore,  the  instantaneous 
power  in  the  circuit  including  all  three  phase  currents  is 

Power  =  \(vahia  ~  vbcic)  +  (vbc  (ia  4-  ib)  +  vabia)  +  (~vb<4c  -  vab  (ib  +  ic))\  /3  (4.34) 

Determine  Fault  Direction  If  no  fault  is  detected  by  the  fault  detection  routine,  the 
power  value  is  added  to  a  buffer  of  the  last  16  power  values  with  no  faults  detected.  This 
buffer  will  be  used  to  calculate  the  average  power  which  will  be  used  for  fault  direction 
determination  when  a  fault  is  detected. 

If  a  fault  is  detected,  the  calculated  power  is  subtracted  from  the  average  of  the  buffer 
to  determine  a  A  power  as  shown  in  Equation  4.35.  The  equation  used  calculate  A  power  is 

&  Power  =  -  Power.  (4.35) 

If  the  A  power  is  positive,  indicating  an  decrease  or  reversal  in  electrical  energy  flowing 
through  the  circuit  breaker,  the  Fault  Direction  flag  is  set  to  Upline,  otherwise  it  is  set  to 
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Downline. 


4.3.2  Example  of  MFM  HSR  Algorithm 

To  illustrate  the  operation  of  the  HSR  algorithm,  the  following  examples  are  presented. 
The  values  for  current  are  assumed  values  and  are  given  only  to  provide  an  indication  of  the 
changes  in  magnitudes  of  the  current  at  various  points  in  the  circuit.  Figure  4-3  shows  the 
circuit  at  some  instant  of  time  with  two  ideal  power  sources  and  three  loads  at  a  steady- 
state  condition  with  no  faults  in  the  circuit.  The  three  loads  are  each  supplied  from  a 
separate  switchboard.  The  circuit  breakers  between  the  switchboards  are  each  monitored 
by  separate  MFMs,  each  with  two  channels  monitoring  either  side  of  the  circiut  breaker. 
The  table  shows  the  output  variables  and  how  they  change  between  no  fault  and  fault 
conditions.  Arrows  are  used  to  give  an  indication  of  the  relative  change  to  the  steady-state 
values.  With  no  faults  in  the  circuit,  the  fault  direction  determination  is  not  performed, 
therefore  ”N/A”  is  entered  into  the  table. 
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Figure  4-3:  MFM  HSR  Algoritm  Example:  Steady  State  Conditon.  Two  ideal  power  sources 
supply  the  loads  off  of  the  switchboards  in  a  steady-state  condition  with  no  faults  in  the 
circuit. 

Figure  4-4  illustrates  wrhat  outputs  would  be  from  the  MFM  HSR  algorithms  if  a  fault 
exists  on  Switchboard  #1.  In  this  simple  example,  the  current  is  assumed  to  flow  almost 
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entirely  to  the  fault,  therefore  the  load  currents  decrease  to  near  zero.  As  can  be  seen, 
the  current  direction  through  MFM  #1  is  reversed  from  steady-state  and  has  increased  in 
magnitude.  The  Park’s  Transformation  for  MFM  #1  will  result  in  an  increase  in  magnitude 
and  a  large  change  in  the  angle  calculation  (it  is  shown  negative  assuming  that  the  steady- 
state  vector  was  in  quadrant  I  and  the  reversal  of  the  current  cause  the  shift  of  the  vector 
towards  quadrant  IV).  The  large  increase  in  the  magnitude  and  the  large  change  in  the  angle 
will  cause  both  fault  flags  to  be  set,  therefore  the  direction  determination  will  be  performed. 
Since  the  current  flow  is  reversed  through  MFM  #1  with  no  change  in  the  voltage  phase, 
the  instaneous  power  calculated  will  also  be  negative.  Calculating  A pawer  using  Eq.  4.35 
results  in  a  positive  value  which  indicates  the  fault  is  ’’upline”  or  upstream  of  the  MFM  #1 
in  relation  to  the  defined  norm.al  power  flow. 

The  fault  causes  the  power  through  MFM  #2  to  increase.  If  the  fault  is  large  enough, 
the  maginitude  and  angle  flags  will  be  set  due  to  the  large  change  in  the  associated  values. 
Therefore  the  direction  determination  will  also  be  performed  in  MFM  #2.  Since  the  power 
only  increased  and  didn’t  reverse,  A  power  will  be  negative,  indicating  to  the  MFM  that  the 
fault  is  ’’downline”  or  downstream  from  MFM  #2  in  relation  to  the  defined  nonnal  power 
flow. 

Figure  4-5  illustrates  a  fault  on  Switchboard  #2.  The  outputs  for  this  scenario  for  each 
MFM  will  be  similar  to  those  for  MFM  #2  for  the  fault  on  Switchboard  #1.  For  MFM  #2, 
it  will  see  no  difference  between  the  two  fault  locations,  both  are  downstream.  For  MFM 
#1,  the  power  will  not  reverse,  therefore  it  will  also  indicated  that  the  fault  is  downstream. 
The  magnitude,  angle,  A  Angle,  and  power  will  increase  on  both  MFMs. 

The  last  scenario  presented  is  that  of  a  fault  on  MFM  #2  itself,  between  the  sensor 
locations.  Figure  4-6  shows  the  system  diagram  and  the  output  table  for  this  scenario. 
MFM  #1  will  see  the  same  indications  as  the  fault  on  Switchboard  #2  and  will  indicate  the 
fault  is  downstream.  MFM  #2  will  not  show  the  same  indications  on  both  of  its  sensor  taps. 
For  the  #1  tap,  the  power  will  increase  in  the  same  direction  as  steady-state,  therefore  tap 
#1  will  indicate  that  the  fault  is  downstream.  The  power  flow  through  tap  #2  will  reverse, 
and  similar  to  MFM  #1  in  the  scenario  with  the  fault  on  Switchboard  #1,  tap  #2  will 
indicate  that  the  fault  is  upstream. 
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Figure  4-4:  MFM  HSR  Algoritm  Example:  Fault  on  Switchboard  #1.  The  fault  causes  the 
power  flow  through  MFM  #1  to  reverse.  MFM  #1  indicates  that  the  fault  is  UPLINE. 
The  fault  causes  an  increase  in  the  power  flow  through  MFM  #2  in  the  same  direction  as 
steady-state.  MFM  #2  indicates  that  the  fault  is  DOWNLINE. 

4.4  Integration  of  NILM  with  the  MFM 

The  MFM  system  uses  similar  input  signals  to  NILM  and  both  systems  have  similar  basic 
functionality,  to  monitor  the  electrical  power  being  used  by  the  electrical  loads  in  the  system. 
The  main  focus  of  the  NILM  system  is  to  determine  detailed  operating  characteristics  and 
conditions  of  the  equipment.  The  main  focus  of  the  MFM  system  is  for  protection  of  the 
electrical  system  during  fault  conditions  and  to  isolate  faults  such  that  as  much  of  the 
distribution  system  is  maintained  in  an  operable  condition  as  possible. 

The  MFM  system  must  operate  very  rapidly  to  be  effective  in  fault  isolation.  One  area 
where  the  NILM  system  could  provide  a  useful  benefit  is  to  maintain  state  information 
on  major  loads  in  the  electrical  distribution  system.  This  state  information  will  provide 
two  benefits  to  the  operation  of  the  MFM.  First,  if  all  of  the  largest  loads  in  the  system 
are  running  and  a  large  increase  in  power  is  detect  by  the  MFM,  the  MFMs  could  take 
action  more  quickly  because  time  would  not  have  to  be  taken  to  determine  if  the  increase 
in  power  is  a  piece  of  equipment  or  a  fault.  Secondly,  the  actual  current  setpoints  could  be 
dynamically  set  based  on  the  operating  conditions  of  the  plant  equipment. 
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Figure  4-5:  MFM  HSR  Algoritm  Example:  Fault  on  Switchboard  #2.  The  fault  causes  an 
increase  in  the  power  flow  through  MFMs  #1  &  #2  in  the  same  direction  as  steady-state. 
MFMs  #1  &  #2  indicate  that  the  fault  is  DOWNLINE. 
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Figure  4-6:  MFM  HSR  Algoritm  Example:  Fault  on  MFM  #2.  The  fault  causes  an  increase 
in  the  power  flow  through  MFM  #1  in  the  same  direction  as  steady-state.  MFM  #1 
indicates  that  the  fault  is  DOWNLINE.  MFM  #2  tap  #1  senses  an  increase  in  power  flow  in 
the  same  direction  as  steady-state,  therfore  tap  #1  indicates  that  the  fault  is  DOWNLINE. 
MFM  #2  tap  #2  senses  increase  in  power  flow  opposite  in  direction  to  steady-state,  therfore 
tap  #2  indicates  that  the  fault  is  UPLINE 
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Chapter  5 


Conclusions  and  Future  Work 


5.1  Conclusions 

The  NILM  system  has  proven  to  be  able  to  provide  useful  data  in  the  monitoring  of  ship¬ 
board  systems.  Through  the  use  of  a  GUI  the  shipboard  crews  will  gain  the  ability  to 
operate  the  system  during  operation,  providing  them  with  diagnostic  indicators  and  system 
status  while  the  shipboard  systems  are  in  operation. 

The  ISA  and  3SA  NILM  systems  installed  at  LBES  have  shown  in  preliminary  testing 
to  be  able  to  distinguish  the  individual  componenet  signals  from  the  aggregate  signal  of 
all  loads.  Further  testing  is  required  when  the  LBES  facility  is  in  full  operation  to  assess 
the  ability  of  the  ISA  &  3SA  NILMs  to  distinguish  the  individual  equipment  operating 
characteristics  in  the  aggregate  signals. 

The  cause  of  the  abnormal  operation  of  the  LPACs  at  LBES  was  determined  to  be  the 
installation  of  new  site  air  compressors  which  are  set  to  operate  in  the  normal  pressure 
ranges  of  the  LBES  LPACs.  The  site  air  compressors  will  supply  the  LP  Air  loads  of  LBES 
whenever  their  pressure  exceeds  that  of  the  LBES  LPACs. 

5.2  Future  Work 
5.2.1  LBES 

There  is  a  need  to  further  explore  utilizing  the  aggregate  signals  at  the  LBES  facility  to 
characterize  the  operating  states  of  the  equipment.  This  analysis  can  be  compared  to  the 
states  determined  by  the  NILMs  installed  on  the  individual  equipment.  Automated  routines 
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could  be  developed  to  provide  realtime  indication  of  plant  conditions.  This  would  facilitate 
the  NILM  being  incorporated  into  the  Next  Generation  Integrated  Power  System  (NGIPS) 
on  future  Navy  ships. 

5.2.2  NILM  GUI 

Further  work  is  needed  to  implement  a  simple  GUI  for  the  NILM  systems  installed  aboard 
ships.  The  crews  will  need  to  have  something  they  can  use  without  requiring  extensive 
training.  The  information  presented  should  be  informative  and  the  operation  of  the  GUI 
should  be  intuitive. 


72 


Bibliography 


[1]  CDR  John  V.  Amy  Jr.  IJSN.  Considerations  in  the  Design  of  Naval  Electric  Power 
Systems.  In  Power  Engineering  Society  Summer  Meeting ,  2002  IEEE ,  volume  1,  pages 
331-335,  Chicago,  IL,  USA,  July  2002.  IEEE. 

[2]  CDR  John  V.  Amy  Jr.  USN  (Ret.)  Ph.D.  P.E.  Modern,  High-Converter-Populations 
Argue  for  Changing  How  To  Design  Naval  Electric  Power  Systems.  Technical  report, 
BMT  Syntek  Technologies,  USA,  2005. 

[3]  Patrick  Lawrence  Bennett.  Using  the  Non-Intrusive  Load  Monitor  for  Shipboard  Su¬ 
pervisory  Control.  Master’s  thesis,  Massachusetts  Institute  of  Technology,  2007. 

[4]  Barry  W.  Boehm.  A  Spiral  Model  of  Software  Development  and  Enhancement.  Com¬ 
puter,  21(5)  :6 1—72,  1988. 

[5]  R.  J.  Bost,  J.  G.  Mellis,  and  P.  A.  Dent.  Is  the  Navy  Serious  About  Reducing  Manning 
on  its  Ships?  Office  of  Naval  Research ,  1999. 

[6]  NSWCCD-SSES  Code98/9345.  A  Guide  To  MEM  III  Operation  (DRAFT).  Techni¬ 
cal  report,  NAVSEA  Naval  Surface  Warfare  Center  Carderock  Division  Naval  Ships 
Systems  Engineering  Station,  Philadelphia,  Pennsylvania,  2003. 

[7|  R.  W.  Cox.  Minimally  Intrusive  Strategies  for  Fault  Detection  and  Energy  Monitoring . 
PhD  thesis,  MIT,  Aug  2006. 

[8]  R.  W.  Cox,  G.  Mitchell,  P.  Bennett,  M.  Piber,  J.  Paris,  W.  Wichakool,  and  S.  B.  Leeb. 
Improving  Shipboard  Maintenance  Practices  Using  Non-Intrusive  Load  Monitoring.  In 
Proceedings  ASNE  Intelligent  Ships  VII,  Philadelphia,  PA. 


73 


[9]  R.  W.  Cox,  J.  P.  Mosman,  D.  McKay,  S.  B.  Leeb,  and  T.  J.  McCoy.  Diagnostic  Indi¬ 
cators  for  Shipboard  Cycliing  Systems  Using  Non-Intrusive  Load  Monitoring.  ASNE 
Day  2006,  2006. 

[10]  T.  DeNucci,  R.  W.  Cox,  S.  B.  Leeb,  J.  Paris,  T.  J.  McCoy,  C.  R.  Laughman,  and 
W.  C.  Greene.  Diagnostic  Indicators  for  Shipboard  Systems  Using  Non-Intrusive  Load 
Monitoring.  In  Proceedings  IEEE  Electric  Ship  Technologies  Symposium ,  Philadelphia, 
PA. 

[11]  CAPT  Norbert  Doerry.  Zonal  Ship  Design.  American  Society  of  Naval  Engineers  Naval 
Engineer's  Journal ,  118(1):39~53,  Winter  2006. 

[12]  CAPT  Norbert  Doerry.  Designing  Electrical  Power  Systems  for  Survivability  and 
Quality  of  Service.  American  Society  of  Naval  Engineer's  Naval  Engineer'' s  Journal , 
1 19(2): 25-34,  2007. 

[13]  CAPT  Norbert  Doerry.  Next  Generation  Integrated  Power  System  NGIPS  Technology 
Development  Roadmap.  Technical  report,  Naval  Sea  Systems  Command,  November 
2007. 

[14]  CAPT  Norbert  H.  Doerry.  Systems  Engineering  and  Zonal  Ship  Design.  Technical 
report,  Naval  Sea  Systems  Command,  May  2006.  Presented  at  the  ASNE  Engineering 
the  Total  Ship  Symposium,  ETS  2006,  May  1-2,  2006. 

[15]  CAPT  Norbert  H.  Doerry  and  David  H.  Clayton.  Shipboard  Electrical  Power  Quality 
of  Service.  Technical  report,  IEEE,  July  2005.  Presented  at  the  IEEE  Electric  Ship 
Technologies  Symposium,  Philadelphia,  PA,  July  25-27,  2005. 

[16]  CAPT  Norbert  H.  Doerry  and  Howard  Fireman.  Designing  All  Electric  Ships.  In 
Proceedings  fo  the  Ninth.  International  Marine  Design  Conference ,  May  2006. 

[17]  ENVIROVAC  INC.,  1260  Turret  Drive,  Rockford,  Illinois  61111.  Instruction  Manual 
for  the  Envirovac  Vacuum.  Collection  System,  United  States  Coast  Guard  270'  B  Class 
WMEC,  Technical  Publications  Number  2839 ,  April  1983. 

[18]  William  C.  Greene.  Evaluation  of  Non-Intrusive  Load  Monitoring  for  Condition  Based 
Maintenance  Applications  of  U.S.  Navy  Propulsion  Plants.  Master's  thesis,  MIT,  2005. 


74 


[19]  Henry  Hegner  and  Bipin  Desai.  Integrated  Fight  Through  Power.  Technical  report, 
IEEE,  2002. 

[20]  Paul  C.  Krause,  Oleg  Wasynczuk,  and  Scott  D.  Sudhoff.  A  nalysis  of  Electric  Machinery 
and  Drive  Systems.  Piscataway,  NJ  :  IEEE  Press  ;  New  York  :  Wiley  Interscience, 
c2002.,  2nd  edition,  2002. 

[21]  Mark  W.  Maier  and  Eberhardt  Rechtin.  The  Art  of  Systems  Architecting,  Second 
Edition .  CRC  Press,  2002. 

[22]  Timothy  J.  McCoy.  Manning  and  Automation.  M  IT  Power  Point  Presentation,  March 
2007. 

[23]  Thomas  Duncan  McKay.  Diagnostic  Indicators  for  Shipboard  Mechanical  Systems  Us¬ 
ing  Non-Intrusive  Load  Monitoring.  Master's  thesis,  Massachusetts  Institute  of  Tech¬ 
nology,  2006. 

[24]  J.  R  Mosman.  Evaluation  of  Non-Intrusive  Load  Monitoring  for  Shipboard  Cycling 
System  Diagnostics.  Master's  thesis,  Massachusetts  Institute  of  Technology,  June  2006. 

[25]  Naval  Sea  Systems  Command.  Technical  Manual  for  AGQIJ+ORF  Ship  Service  Gas  Tur¬ 
bine  Generator  (SSGTG),  Description ,  Operation,  and  Scheduled  Maintenance,  S931 1- 
XX-MMO-010 ,  vol  1  edition. 

[26]  Naval  Sea  Systems  Command.  Technical  Manual  for  Air  Compressor  Low  Pressure, 
Oil  Free,  S6226-DE-MMA-010/07524 ,  November  1992. 

[27]  Naval  Sea  Systems  Command.  Naval  Ships  ’  Technical  Manual  Chapter  320,  Electrical 
Power  Distribution  Systems ,  revision  2  edition,  April  1998. 

[28]  Naval  Sea  Systems  Command.  Technical  Manual  for  Machinery  Control  System,  Shaft 
Control  Unit,  P/N  63E915603G1,  S9202-A N-M M F- 020,  volume  2  edition,  Oct  1998. 

[29]  Naval  Sea  Systems  Command.  Technical  Manual  Organizational  Level  Maintenance 
LM 2500  Propulsion  Gas  Turbine  Module  Description  and  Operation,  S9234-AD-MMO- 
010/LM2500 ,  revision  3  edition,  May  2000. 

[30]  Naval  Sea  Systems  Command.  Compressor,  Air,  Low  Pressure,  Oil-Free,  Model  Star 
200,  S6220-EE-MMA-010 ,  change  f  edition,  October  2002. 


[31]  Naval  Sea  Systems  Command*  IJnverval  Engine  Controller  (UEC)  Operations  and 
Maintenance ,  S6263-B7-MMA-010 ,  vol  1  &  vol  2,  change  a  edition,  Aug  2003. 

[32]  Mark  A.  Fiber.  Improving  Shipboard  Maintenance  Practices  Using  Non-Intrusive  Load 
Monitoring.  Master’s  thesis,  Massachusetts  Institute  of  Technology,  2007. 

[33]  Lyle  Plesnick,  Tracy  Hannon,  and  Daniel  Devine.  An  Intelligent  Fault  Detection  De¬ 
vice  For  Shipboard  Power  Systems.  Technical  report,  Naval  Surface  Warfare  Center, 
Carderock  Division  /  U.S.  Navy. 

[34]  Ethan  Proper.  Automated  Classification  of  Power  Signals.  Master’s  thesis,  Mas¬ 
sachusetts  Institute  of  Technology,  June  2008. 

[35]  J.  S.  Ramsey,  S.  B.  Leeb,  T.  DeNucci,  J.  Paris,  M.  Obar,  R.  W.  Cox,  C.  R.  Lauhman, 
and  T.  J.  McCoy.  Shipboard  Applications  of  Non-Intrusive  Load  Monitoring.  In 
Proceedings  ASNE  Reconfiguration  and  Sui'vivability  Symposium ,  A  tlandic  Beach ,  FL , 
2005. 

[36]  Jack  S.  Ramsey.  Shipboard  Applications  of  Non-Intrusive  Load  Monitoring.  Master’s 
thesis,  Massachusetts  Institute  of  Technology,  2004. 

[37]  S.  R.  Shaw.  System  Identification  Techniques  and  Modeling  for  Non-Intrusive  Load 
Diagnostics .  PhD  thesis,  MIT,  February  2000. 

[38]  Systems  Engineering ,  Principles  and  Practice .  Wiley  Interscience,  2003. 


76 


Appendix  A 


NILM  Transducer  Data  Sheets 
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Voltage  Transducer  LV  25-P 

For  the  electronic  measurement  of  voltages  :  DC,  AC,  pulsed..*, 
with  a  galvanic  isolation  between  the  primary  circuit  (high  voltage) 
and  the  secondary  circuit  (electronic  circuit). 


Electrical  data 

Primary  nominal  r.m.s.  current 

10 

mA 

K 

Primary  current,  measuring  range 

0  ..  t 

14 

mA 

R« 

Measunng  resistance 

p 

w  *n*> 

p 

with  ±  1 2  V  @  ±  10  mA  M 

30 

190 

U 

30 

100 

u 

with  1 15  V  CajtlGmA^ 

100 

350 

u 

@±14mAra> 

100 

190 

Q 

l$N 

Secondary  nominal  r.m.s.  current 

25 

mA 

K» 

Conversion  ratio 

2500  ; 

1000 

V, 

Supply  voltage  (±  5  %) 

±  12  .. 

15 

V 

•c 

Current  consumption 

10(@±15V)+IS 

mA 

v. 

R.m.s.  voltage  for  AC  isolation  test t;.  50  Hz,  1  mn 

2.5 

kV 

c 

Accuracy  -  Dynamic  performance  data 

J 

X, 

Overall  Accuracy  @  lPN  ,  TA  =  25X  @±12..  15  V 

±0.9 

% 

@  ±  15  V  (±  6  %) 

±0.8 

% 

£. 

Linearity 

<0.2 

% 

Typ 

Max 

to 

Offset  current  @  ~  0,  TA  *  25X 

±0  15 

mA 

'o, 

Thermal  drift  of  l0  OX  +  25X 

±  0.06 

±0.25 

mA 

+  25X  ..  +  ro°c 

±0.10 

±0.35 

mA 

t. 

Response  time21  @  90  %  of  VPw 

40 

MS 

General  data 


T* 

Ambient  operating  temperature 

0  ..  ♦  70 

X 

T, 

Ambient  storage  temperature 

-  25  ..  +  85 

X 

Primary  coil  resistance  @  TA  -  70*X 

250 

a 

Rs 

Secondary  coil  resistance  @  TA  ~  7QX 

110 

a 

m 

Mass 

22 

g 

Standards 

EN  50178 

Notes  :  Between  primary  3nd  secondary 

*  R  ,  =  25  kli  {UR  constant,  produced  by  the  resistance  and  inductance 
of  the  pnmary  circuit) 

31  A  list  of  corresponding  tests  is  available 


lPN  =  10  mA 
Vw  =  10  ..  500  V 


Features 

•  Closed  loop  (compensated)  voyage 
transducer  using  the  Hall  effect 

•  Insulated  plastic  case  recognised 
according  to  UL  94-VO. 

Principle  of  use 

•  For  voltage  measurements,  a  current 
proportional  to  the  measured  'voltage 
must  be  passed  through  an  external 
resistor  R  .  which  is  selected  by  the 
user  and  installed  in  series  with  the 
primary  circuit  of  the  transducer. 

Advantages 

•  ExceKent  accuracy 

•  Very  good  linearity 

•  tow  thermal  drift 

•  Low  response  time 

•  High  bandwidth 

•  High  immunity  to  external 
interference 

•  Low  disturbance  in  common  mode. 

Applications 

•  AC  variable  speed  drives  and  servo 
motor  drives 

•  Static  converters  for  DC  motor  drives 

•  Battery  supplied  applications 

•  Uninterruptible  Power  Supplies 
(UPS) 

•  Power  supplies  for  welding 
applications. 


981  125/14 
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Dimensions  LV  25-P  (in  mm.  1  mm  =  0.0394  inch) 


Bottom  view 


\  •; 


l\»!  r\|< 

on! 

M  CM 


t  \ 


26 

.1524 


2*  00635mm 


3  <8 1mm 


Back  view 


Right  view 

4 «/  0.3  . 


Top  view 


LCM® swi” 

made 


C€ 

LV  25-P 
00-00  00 


M 


/  \  v 

Standard  00  Year  Ww* 

or  N«  SP 


Secondary  terminals 

Terminal  ♦  :  supply  voltage  +  12 ..  15  V 

Terminal  M  :  measure 

Terminal  -  :  supply  voltage  -  12 ..  15  V 

Connection 

-mT 


LV  ?5  I*  M ?  {,}  C 

,r,  J. 2: _ 


Mechanical  characteristics 


Remarks 


•  General  tolerance 

*  Fastening  &  connection  of  primary 

•  Fastening  &  connection  of  secondary 

*  Recommended  PCB  hole 


t  0.2  mm 

2  pins 

0.635  x  0.635  mm 

3  pins  0  1  mm 
1.2  mm 


•  ls  is  oosftive  when  Vf,  is  applied  on  terminal  +HT 

•  This  is  a  standard  model.  For  different  versions  (supply 
voltages,  turns  ratios,  unidirectional  measurements.,.), 
please  contact  us. 


Instructions  for  use  of  the  voltage  transducer  model  LV  25 -P 

Primary  resistor  R ,  :  the  transducers  optimum  accuracy  is  obtained  at  the  nominal  primary  current  As  far  as  possible,  R  ,  should  be 
calculated  so  that  the  nominal  voltage  to  be  measured  corresponds  to  a  primary  current  of  10  mA. 

Example:  Vcrftage  to  be  measured  Vw  *  250  V  a}  R ,  -  25  kii  /2.5  W.  \p  -  10  mA  Accuracy  -  ±  0.8  %  of  (@  TA  »  ♦  25*0) 

b)  R .  ~  50  kii/ 1.25  W.  Ip  -  5  mA  Accuracy  =  ±  1.6  %  of  TA  ®  +25*C) 

Operating  range  (recommended)  .  taking  mto  account  the  resistance  of  the  primary  windings  (which  must  remain  low  compared  to  R ,  in  order 
to  keep  thermal  deviation  as  low  as  possible)  and  the  isotedion  this  transducer  is  suitable  for  measuring  nominal  voltages  from  10  to  500  V, 


LEM  re  serves  the  right  to  carry  out  modifications  on  its  transducers,  in  order  to  improve  them,  without  previous  notice. 
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Current  Transducer  LA  55-P 

For  the  electronic  measurement  of  currents  :  DC,  AC,  pulsed  .... 
with  a  galvanic  isolation  between  the  primary  circuit  (high  power) 
and  the  secondary  circuit  (electronic  circuit). 


K* 

Primary  nominal  r.m.s.  current 

50 

A 

K 

Primary  current,  measuring  range 

0  . 

±70  A 

K 

Measuring  resistance  @ 

ta  =  w-c 

T»  =  85X 

with  ±  12  V  @±50A^ 

10  100  I 

60  95  a 

@  ±  70  A  ^ 

10  50 

60’>  60n  a 

with  ±  15  V  (a)  ±  50  A  ^ 

50  160 

135  155  a 

50  90 

135*  135=3  n 

Secondary  nominal  r.m.s.  current 

50 

rnA 

Ks 

Conversion  ratio 

1  : 

1000 

Vc 

Supply  voltage  (±  5  %) 

±12..  15  V 

'c 

Current  consumption 

10{@±15V)+ts  mA 

R.m.s.  voltage  for  AC  isolation  test.  50  Hz. 

1  mn  2.5 

kV 

C 

Accuracy  -  Dynamic  performance  data 

1 

..  j 

X 

Accuracy  @  .  TA  *  2S°C  @±16V(±5%) 

±0.65 

% 

@±12..  15  V  (±5%) 

±0.90 

% 

eL 

Linearity 

<  0.15 

% 

Typ 

Max 

«o 

Offset  current  @  lp  *  0,  25°C 

±0.2 

mA 

W. 

Residual  current  @  ip  =  0,  after  an  overload  of  3  x  lpw 

±0.3 

mA 

i0, 

Thermal  drift  of  l0  0°C  ..  ±  70X 

±0.1 

±0.5 

mA 

-  25°C  ..  ♦  85°C 

±0.1 

±0.6 

mA 

t„ 

Reaction  time  @  10  %  of  lPw4ft 

<  500 

ns 

t. 

Response  time  @  90  %  of ^ 

<  1 

ps 

dl/dt  di/dt  accurately  followed 

>  200 

A/ps 

f 

Frequency  bandwidth  (-  1  dB) 

DC  ..  200 

kHz 

General  data 


T* 

Ambient  operating  temperature 

-25..  +  85 

°C 

T, 

Ambient  storage  temperature 

-  40  ..  ♦  90 

JC 

Secondary  cod  resistance  @ 

T,  =  70‘C 

80 

a 

Ta  =  85-C 

85 

n 

m 

Mass 

18 

g 

Standards  4 

EN  50178 

Notes  :  Measuring  range  limited  to  ±  60  A  ^ 

2:  Measuring  range  limited  to  ±  55  A 
3:  Result  of  the  coercive  ftekJ  of  the  magnetic  circuit 
A  list  of  corresponding  tests  is  available 


LEM  Components 


lPN  =  50  A 


Features 

•  Closed  loop  (compensated)  current 
transducer  using  the  HaB  effect 

•  Printed  circuit  board  mounting 

•  Insulated  piastic  case  recognized 
according  to  UL  94- VO 

Advantages 

•  Excellent  accuracy 

•  Very  good  linearity 

•  Low  temperature  drift 

•  Optimized  response  time 

•  Wide  frequency  bandwidth 

•  No  insertion  losses 

•  High  immunity  to  external 
interference 

•  Current  overload  capability 

Applications 

•  AC  variable  speed  drives  and  servo 
motor  drives 

•  Static  converters  for  DC  motor  drives 

•  Battery  supplied  applications 

•  Uninterruptible  Power  Supplies 
(UPS) 

•  Switched  Mode  Power  Supplies 
(SMPS) 

•  Power  supplies  for  welding 
applications. 


980706/8 
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Dimensions  LA  55-P  (in  mm.  I  mm  -  0  0394  inch) 


Bottom  view 


. t. 


12.7 


LEM®  LA  55-P 


C€ 


—f — 

i 


SWISS 

made 


M  —  00-00  00  + 

xi  t~ 


\ 


4-  5^Q8_  .  2?  86 
36.5 


/ 

Standard  00 
or  N5  SP.. 

Front  view 


\  v 

Year  Week 


Left  view 


t 


1  1 

r 

o.f 

1 

H 

r 

j 

— 

>. . \. .  s 

4 

_ . 21 A . _ 

Secondary  terminals 


Terminal  + 
Terminal  - 
Terminal  M 


supply  voltage  ♦  12  ..  15  V 
supply  voltage  -  12  ..  15  V 
measure 


Connection 


.A  55-P 


JY  J- . \S)  -  -  CH  o  ov 


Mechanical  characteristics 

Remarks 

•  General  tolerance 

±02  mm 

*  L  is  positive  when  lp  flows  In  the  direction  of  the  arrow 

♦  Primary  througlvhole 

12.7  x  7  mm 

•  Temperature  of  the  primary  conductor  should  not  exceed 

•  Fastening  &  connection  of  secondary 

3  pins 

90‘C. 

0.63  x  0.56mm 

•  Dynamic  performances  (di/dt  and  response  time)  are  best 

Recommended  PCB  hole 

0.9  mm 

with  a  single  bar  completely  filling  the  primary  hole. 

•  In  order  to  achieve  the  best  magnetic  coupling,  the  primary 
wrndmgs  have  to  be  wound  over  the  top  edge  of  the  device. 

•  This  is  a  standard  model.  For  different  versions  (supply 
voltages,  turns  ratios,  unidirectional  measurements...), 
please  contact  us 

LEM  reserves  the  rignt  to  carry  out  modifications  on  its  transducers,  in  order  to  improve  them,  without  previous  notice. 
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Current  Transducer  LA  205-S 

For  the  electronic  measurement  of  currents  :  DC,  AC,  pulsed..., 
with  a  galvanic  isolation  between  the  primary  circuit  (high  power) 
and  the  secondary  circuit  (electronic  circuit). 


r 

Electrical  data 

i 

u 

Primary  nominal  r.m.s.  current 

200 

A 

•p 

Primary  current,  measuring  range 

0.. 

±300 

A 

Lm 

^  Measuring  overload  1! 

600 

A 

R„ 

Measuring  resistance  @ 

Ta  «  70sC  I 

|  Ta  =  85X 

M 

K 

with  ±  12  V  @  1 200  A  ^ 

0 

68 

0  86 

a 

@±300Am» 

0 

33 

0  30 

a 

with  ±  1 5  V  @±2WA^ 

5  ’ 

95 

5  93 

a 

@*.300  A  ^ 

5 

50 

5  49 

a 

>*, 

Secondary  nominal  r.m.s.  current 

100 

mA 

Conversion  ratio 

1  ; 

2000 

Vc 

Supply  voltage  (±  5  %) 

±  12  . 

15 

V 

*c 

Current  consumption 

20*®±15V)+I 

mA 

V, 

R.m.s  rated  voltage  **  safe  separation 

1625 

V 

basic  isolation 

3250 

V 

j  Accuracy  -  Dynamic  performance  data 

X5 

Overall  accuracy  @  lPH  ,  TA  «  25X 

±0,8 

% 

A 

Linearity 

<0.1 

% 

!  Typ 

Max 

1, 

Offset  current  @  C  ~  9,  TA  ~  25° C 

±0.15 

mA 

i0M 

Residual  current ‘*)@  l„  =  0.  after  an  overload  of  3  x 

±0,50 

mA 

U 

Thermal  drift  ofl0  -  10X  .. 

+  85X  j 

|±  0.15 

±0.30 

mA 

tfS 

Reaction  time  @  10  %  of  \p  ((tix 

<500 

ns 

tf 

Response  time  **  %  90  %  of  lp 

<  1 

ps 

di/dt  di/dt  accurately  followed 

>  100 

Ayps 

f 

Frequency  bandwidth  (-  3  d8) 

DC 

..  100 

kHz 

{ — 
I 

General  data 

j 

T* 

Ambient  operating  temperature 

-10  ..  +  85 

X 

Ts 

Ambient  storage  temperature 

-  40  ..  +  90 

X 

", 

Secondary  coil  resistance  @  TA 

=  70X 

35 

y 

ta 

»asx 

37 

a 

m 

Mass 

110 

s 

Standards  5) 

£N  50178 

Notes  :  3  mn/hour  @  Vc  ~  z  15  V,  5  Q 

Tt  Pollution  class  nr  2.  With  a  non  insulated  primary  bar  which  fills  the 


through-hole 

The  result  of  the  coercive  field  of  the  magnetic  circuit 
With  a  di/dl  of  100  A/ps 
S}  A  fist  of  corresponding  tests  is  available 


LEM  Components 


I PN  »  200  A 


Features 

«  Closed  loop  (compensated)  current 
transducer  using  the  Hall  effect 

*  Insulated  plastic  case  recognized 
according  to  UL  94- VO 

*  Patent  pending. 

Advantages 

•  Excellent  accuracy 

*  Very  good  linearity 

♦  Low  temperature  drift 

*  Optimized  response  time 

*  Wide  frequency  bandwidth 

♦  No  insertion  losses 

•  High  immunity  to  external 
interference 

♦  Current  overload  capability 

Applications 

•  AC  variable  speed  drives  and  servo 
motor  drives 

*  Static  converters  for  DC  motor  drives 

♦  Battery  supplied  applications 

*  Uninterruptible  Power  Supplies 
(UPS) 

•  Switched  Mode  Power  Supplies 
(SMPS) 

•  Power  supplies  for  welding 
applications. 


980716/5 
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Dimensions  LA  205-S  (in  mm.  1  mm  =  0.0394  inch) 


Mechanical  characteristics  Remarks 

♦  General  tolerance  ±05  mm  •  ls  is  positive  when  1^  flows  in  the  direction  of  the  arrow. 

♦  Fastening  2  holes  0  5.5  mm  •  Temperature  of  the  primary  conductor  should  not  exceed 

♦  Primary  through-hole  23  x  18  mm  10QX. 

♦  Connection  of  secondary  Mofex  5G46-04/AG  •  Dynamic  performances  (di/dt  and  response  time)  are  best 

with  a  single  bar  completely  filling  the  primary  hole. 

•  This  is  a  standard  model.  For  different  versions  (supply 
voltages,  turns  ratios,  unidirectional  measurements.  ), 
please  contact  us. 


LEM  reserves  the  right  to  carry  out  modifications  on  its  transducers,  in  order  to  improve  them,  without  previous  notice. 
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Current  Transducer  LA  305-S 

For  the  electronic  measurement  of  currents  :  DC,  AC,  pulsed. 
with  a  galvanic  isolation  between  the  primary  circuit  (high  power) 
and  the  secondary  circuit  (electronic  circuit). 


Electrical  data 

u 

Pnmary  nominal  r  m  s  current 

300 

A 

L 

Primary  current,  measuring  range 

0 

±500 

A 

R« 

Measuring  resistance  @ 

t,  =  rare 

T„  =  85'C 

Ru»„  r 

R„ 

p 

ftr-T-.  rr.ax 

with  ±  12  V  @  ±  300  A  ^ 

0 

52 

0  50 

Q 

@±500A(BM 

0 

17 

0  15 

n 

with  ±  15  V  @±300A^ 

0 

75 

5  73 

a 

@±500A™ 

0 

31 

5  29 

a 

Secondary  nominal  r.m.s.  current 

120 

mA 

Conversion  ratio 

1  : 

2500 

vc 

Supply  voltage  (±  5  %) 

±  12  .. 

15 

V 

■c 

Current  consumption 

20(@±15V)+I 

mA 

R.m  s.  rated  voltage  safe  separation 

1750 

V 

basic  isolation 

3500 

V 

Accuracy  -  Dynamic  performance  data 

Overall  accuracy  @  lPW  ,  TA  »  25*C 

±0.8 

% 

eL 

Linearity 

<0.1 

% 

Typ 

Max 

'c 

Offset  current  @  =  0.  T  =  25°C 

±  0.20 

mA 

■<* 

Residual  current*’  @  lp=  0,  after  an  overload  of  3  x 

±0  40 

mA 

lOT 

Thermal  drift  of  l0  -  10"C  ..  ±  85CC 

±0.12 

±0.30 

mA 

xn 

Reaction  time  @10%  of  lp  ^ 

<  500 

ns 

Response  time 3>  @  90  %  of  fp  ^ 

<  1 

ps 

dl/dt  di/dt  accurately  followed 

>  100 

A/ps 

f 

Frequency  bandwidth  {-  3  dB) 

DC 

.  100 

kHz 

General  data 


T, 

Ambient  operating  temperature 

- 10..  +  85 

°C 

Ambient  storage  temperature 

-40  ±90 

*C 

R, 

Secondary  coil  resistance  @ 

TA  -  70°C  35 

Q 

T.  =  85X  37 

a 

m 

Mass 

200 

9 

Standards  4) 

EN  50178 

■P  N 


300  A 


Features 

•  Closed  loop  (compensated)  current 
transducer  using  the  Hall  effect 

•  Insulated  plastic  case  recognizsd 
according  to  UL  94-VO 

•  Copyright  protected. 

Advantages 

•  Excellent  accuracy 

•  Very  good  linearity 

•  Low  temperature  drift 

•  Optimized  response  time 

•  Wide  frequency  bandwidth 

•  No  insertion  tosses 

•  High  immunity  to  external 
interference 

•  Current  overload  capability. 

Applications 

•  AC  variable  speed  drives  and  servo 
motor  drives 

•  Static  converters  for  DC  motor  drives 

•  Battery  supplied  applications 

•  Uninterruptible  Power  Supplies 
(UPS) 

•  Switched  Mode  Power  Supplies 
(SMPS) 

•  Power  supplies  for  welding 
applications. 


Notes  :  r'  Pollution  class  2.  With  a  non  insulated  primary  bar  which  fills  the 
through-hole 

2)  The  result  of  the  coercive  field  of  the  magnetic  circuit 
*  With  a  di/dt  of  100  A/ps 

4;  A  list  of  corresponding  tests  is  available  980707/6 


LEM  Components 


wwwJem.com 
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Dimensions  LA  305-S  (in  mm  1  mm  =  0  0394  inch) 


2**55 
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. . .  -Lj.4.^. 

Top  view 


Left  view 


Secondary  terminals 


Pm  1  :  supply  voltage  +  12  15  V 

Pin  2  :  measure 
Pin  3  :  supply  voltage  -  12  ..  15  V 
Pin  4  NC 


Connection 


Mechanical  characteristics 

Remarks 

*  Genera!  tolerance 

±0.5  mm 

•  lg  is  positive  when  lP  Nows  in  the  direction  of  the  arrow 

•  Fastening 

2  holes  0  5.5  mm 

•  Temperature  of  the  primary  conductor  should  not  exceed 

•  Primary  through-hole 

25  6  x  25.5  mm 

100X. 

•  Connection  of  secondary 

Motex  5046-04/  AG 

•  Dynamic  performances  (di/dt  and  response  time)  are  best 
with  a  single  bar  completely  filling  the  primary  hole 

•  This  is  a  standard  model.  For  different  versions  (supply 

voltages,  turns  ratios,  unidirectional  measurements  ..}, 
please  contact  us. 

LEM  reserves  the  right  to  carry  out  modifications  on  its  transducers,  in  order  to  improve  them,  without  previous  notice. 
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FLUKE 


i 3000s 

AC  Current  Clamp 

instruction  Sheet 


Introducing  the  13000s 

The  12€€Os  -■*  a  damp-on  AC  Current  Comp-  far  us*  with 
cacwascopes  and  nwww.i  vvlh  bar  and  yipu!  usir.j  y*  suppfed 
Dual  Sanam  to  BhiC  Acs* par.  The  ciamp  Is  optimized  foe 
measurement  on  power  tf=siritwi<>c  systems 

Unpacking 

The  fbHcwig  ten’*  shcad  &e  rrriuded  sn  yoyr  Cwent  Cramp  bo*;. 

*  Cumm?  Clanst 

*  Duo)  Saufii  to  SfC.  Adapter  {.model  PMfiQ&f  <’36  1 ; 

*  I  nstwcson  -Sheet  {this  pape? » 

Oiec*  the  contents  cf  toe  shipping  box  faf'corrpfesefifris.  ,f 
scmeiWrig  to  the  be*  has  beto  damaged  a  contact  your 

ijartejtor  cr  the  nearest  FLilKE  saes  or  sendee  office  imm «S*teiy, 

Safety  Information 


Read  First;  Safety  Information, 
safe  operation  and  service  of  the 


c  urreof&Jamp,  follow  these  tnatnwtions: 


Read  the  operating  e.s&ocrborra  before  use  and  fofow  afl 
safely  mirucbcm* 

Us*  ew  Current  Clamp  any  as  speofotc  n  the  operating 
rnirucftcnv  c«th*rwise  the  darpY  safety  features  may 
not  protect  you 

Adhere  tc  local  and  nadorsau  safety  cooes.  VnraHtftdtad 
protective  eoupmera  must  be  y«»S  to  prevent  shock  and 
arc  hast  «njt#y  when*  haaarctous  -a*  conductors  arv 
exposed 

Oo  not  bold  the  Currant  Camp  an/Nbere  be/onc  th* 
tactile  oa<r«r  se«  Figure  S 

Befcre  each  use  -nspect  the  Current  Clamp,  tot*  for 
crack.*  or  mas-ng  pcrtwita  of  the  clamp  hcytfcftg  >  output 
cable  mulattcn.  Aisc  ->c*  for  loose  or  wtaAenec 
ccnpocwnts.  Pay  particular  abenfecn  tc  to#  stsuiafion 
surrounding  th«ja«&. 

Check  fh*  magnetic  mafcng  surfaces  ctf  the  oamp  ja**s; 
these  should  he  free  c i  dust  dirt,  rtist  and  otlie?  fcrugn 
ih30er 

Never  use  the  damp  cn  a  ercuit  wish  voi'^es  higher  than 

m C  V  CAT  111. 


CAY  fi  i  equipment  is  designed  to  protect  agorst 
tr-ansients  in  equipment  an  §wxf  eq^nent 
inwateftofts.  such  as  cisiflbutGO  park's,  feeders  irq 
short  branch  cireatis.  and  ■i^dmg  sys terns  in  large 
buSdtogs 


use  trireme  caufeon  when  wcning  around  bore 
conductors  or  bus  bars  Contact  w?to  the  aondMctw  wryd 
neasuS  ;n  esectsc  shoe*. 

Use  caution  when  «crt«g  wtfh  sc  tages  above  t5G  V  cfc, 
3G  V  ac  me  of  *2  y‘  ac  pea*  Such  voltages  pose  a  sheet 
hazard 


Symbols 


m 

SAay  be  s.s*«  cn  HAZAPCCiUS  LfVE  c«toucm. 

Produc!  is  protected  by  dou&fe  msdabon. 

A 

ifesit  of  Danger,  itrpert^  rstfortnaben  Set 
•nstr’jcbon  Sheet. 

fesk  of  Eiectfi;  Shoe* 

I^I 

Cdfitorms  to  refevart  £yrop*an  startiards 

rarthgroisnd 

Specifications 

SAFETY 


A 


Input  j#iyj  a 
Output  float  .cog* 


Complies  with  European 
slandartfa  Bi-lEC  9*5  t(M 
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ELECTRICAL  SFEC^fCATfOWA 
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GENERAL 

Clamp.  Dipt«»wns  310  x  12B  v.4S  fr*m 

1 12.2  x  4.7  x  1  &ln> 

:.V«gnt  1200  g  ?42.3oz} 

tTafcl*  length  2 1  tv  <£2.9  .aj 

Mamsum  concwctor  »«*  ^54  ?25  in} 

Maximum  J»v  Opening  9C  mm  (3.5  *n> 

Teinperafe^e 

CtoeraS?ig  - 10  to  +53*C  {+ 14  to  1 122*F) 

Nomop«^3tiriB  HCio  »65*C iHOtoi-IS^Ti 

Reta-  ve  rfemoity 

Cpefateng  0  to  55*4,  ^  fo  -*.3J-4C  ( -*•  &5*F1 

0  to  75%.  up  to  »5C*C  ;+ 1  J2*FJ 
McrUDperatmg  0  to  S5%.  up  to  *  10*C  t + 14*F} 

0  to  75%.  up  to  +32*C ; » *-S4-',F5 

Mftude 

Operasng  to  2000  m  (6500  tl\ 

Non-opera* mg  to  ?2Q£30  m  t:40GKJ  ft; 


Ff^QUEHCY  RE4PONSE 
30A  Range 


Figure  1.  Fregusncry  ftpjpow*  $  tQA 

Vm  Range 


Fsjy#?  2.  Fnequefscy  Response  0  10A 

S»6A  Range 


Figure  3  FetpuertCY  Rtspome  ig  f  8GA 


PHASE  SHIFT 
MA  Range 


Fig $  Ph*s*  Shft  Q  f  OOA 


F* gw*  7  ftiavnMm  Cumatfar  *  FnoOwn  oY  to*  Fra 


instrument  Compatibility 

The  iSOCEs  is  -combats  «  wrtr  any  Fluk*  SospeMet*-  test  tcoi, 

?W*r  Harroooci  Analyzer  Osc^osccpe,  Mufrnufcer.  cr  c’5<*r 

vantage  measurement  dev**  toas  has  the  foUowtng  features. 

*  BNC  input  ccnnecaw.  The  Dti*  Banana  to  BNC  Adapter 
incused  :o  C'»  package  can  be  uses  to  ocroeci  »  standard 
mquts  co  ^ultvn*ren.  For  the  122  senes  SoopeMwers  use 
toe  BB120  Shiecec  Banana  tc  BMC  Adapter 

«  npu  accuracy  of  2\  t-eoer  » take  kill  jovanuge  ol  the 

accuracy  of  the  Current  Clamp 

*  ihpua  nypeoamc*  cr  greater  than  or  equal  to  t  MO.  awi  fcr  fuH 

tandwtotf:  and  accuracy  a  rcaisnwrt  capacity  <v  AT  pF, 

«  A  passband  <4  more  than  four  tines  toe  frequency  of  the 

vwnefcrm  x  &*  measured 


Using  the  Current  damp 

To  use  the  Current  Clarrp  toe-ve  trwsrucsicrtv 
t  Connect  the  SCCOs  Current  Damp  to  the  <sas#*s  input  cr  to# 
mwsunrw  inttramenL  When  ytu  a-*  uang  a  nuttmeter  ute 
the  C  .**  Banana  to  BMC  A«p*®r  iPMSCS  1CC1)  to  corner  the 
Cjrre.it  <Clarrp  » the  input.  sSee  rigur*  S.t 

2.  On  th*  C unent  C i amp.  select  tot  least  sersrbv*  range  '0.1 
mVf.Aj. 

3.  Seiecr  toe  appropruhe  clamp  senstartty  on  you-  SccpeWeter 
test  bo  or  csdlbscoce. 

4.  3osdon  toe  Current  Carp  perpes®cu;ar?y  and  centered 
around  toe  ccnduaor. 

Ma*e  sure  to a:  the  anwt  marked  on  toe  clamp  ;aw  points 
tomard  me  ;cn*ct  ortenswian  tor  correct  phase  drtp!*v  on  toe 
oscilloscope  {See  Figure  9 

Use  the  markings  on  toe  to  center  the  aoncbctor 

5.  Observe  the  cwrrern  vafcre  atso  waveform  on  the  rstrumerr's 
ckspsy. 

§  d  desired,  select  a  >:*«•  range  on  toe  Cjrer?  C  amp  and  set. 
to*  comispc^jmg  senstv-ry  tm^A  setting) on  toe 
Sc  ape  Meier  ter  too  or  cscillcscop* 


if  the  corresponding  sans.wty  is  rwf  avp.able  on  toe  SccceMew 
oroscsccscctpe  afreet  the  closest  seemg  ana  saieMate  the  actual 
currant  value  from  toe  cfcsp  rayed 

Examptel  Q  jnrt  Gamp  W.  a  O.f  mV>A,  3ctx«Me<er  :*v  too 
set  to  t  rr.VVA 

Tt>e  ScopeMetar  test  toe1  csplays  5QA 

The  real  current  vaUe  cor  be  ca<x» x.vc  wdh  toe  faVownng  fcrmu  a 


Aetna)  rvfrtnt  _ 


3w*'f  narue  • 


jerwa»fly  mmtren f 
ieayfiv+r>  CurenfCtoirp 


fntW  -  A 
'■  C.#mV :'  A 


-50CA 


ExantfOt  2.  Current  Cranp  set  p  10  mV; A.  Mjtfcmetev  enspnays 

o.2sev 

Azmat  curmTi . 

3<moiayy»lue  _  0  235V  _  2S5toV  ^ 
5<#ishv»fy  Currenr  Cramp  ttJmV  /  .4  tOmV  ,  A  ” 


A  WamhtQ 

If  the  sensitivity  setting  (rwViA}  of  the  ScopeMet er  test 
tool  or  osoiKostope  does  not  correspond  erth  the 
eetting  ot  the  Current  Clamp,  the  instrument  may 
display  a  much  fcn*«  current  than  th*  actual  value. 
This  may  cause  nr-tstnUrpteiatton  and  as  a 
consequence  inc-orrect  haodhng 


Fyurv  i  Onsncwoo  cf  the  Currertf  Clarrp 


Measurement  Considerations 

Observe  the  loltowmg  ga>u*Urtes  fcr  o^stionirg  toe  Cuneni  Oamr 
Jomvs: 

•  Center  the  .'onductor  «j«»e  toe  carrp  jaws 

•  Moke  sure  toe  camp  s  perperdxxilanc  toe  conductor. 

•  Mafc*  sure  that  the  arw  mahseo  on  tot  jaw  o<  the  ‘Current 
Ci«rp  ports  toward  the  correct  directw  . 

Observe  toe  fcttowmg  guidelines  when  malting  m*3Si***ments: 

•  If  p-c-ssib^,  aacic  measureme nts  cose  to  other  currerw-carytng 
conducion. 

•  Onto*  Cur  rtnlClarp  s*ecMh*Twstappr  oprat*  range  for 
th*  current  a*ng  measured  to  get  to*  best  acouracy 


Afa/n  fenance 

Befcc*  each  vs*,  assure  ccnt-^uec  wk«y  c>y  inspwctng  the  clamp 
took  fcv  o-atts  or  rrvssng  scnfcms  -a!  to*  cianp  ixiuvng  arc  output 
cat**  intUat^g  cow  and  fcr  loose  ey  weakened  ccmipcnenis.  Pay 
pait-rulaf  attention  to  -Jw  .ns^ator  sunc-undnsg  toe  3arc  pws.  I?  a 
6anp  las  ttos  impecsioa  4ap«  r  sht  tc  pr«ve">?  vrirrterped 
eperaton  A  damp  ,n>  warranty  will  be  prompey  repaired  or 
replaces  i'as  FkJk6'sd.>scf€C>on)  are  returned  at  nc' charge 

Claanino  and  Storaioc 

Periodically  rvp*  the  •:««  w*to  a  d#np  ;k«t  detergent  'io  net 
use  abrasives  *  sdventi  Open  toe  jaws  and  rw*  the  "nagorjc 
ede  c-ecei  «ih  a  lighiy  ci»d  ckKh.  Cc  not  »:o«r  rust  a  oorroswn 
to  fcrrri  on  toe  >n*gr*t'5c  cote  ends. 


If  your  Current  Clamp  does  not  work 

f  toe  Current  Cianp  dees  not  perfcrm  propeny.  us«*  to*  fo-ScwnQ 

slept  to  be<c  i*oa!*  toe  oewemr 

•  nspect  to*  |aw  r.urtrg  surface  ky  seankness.  i*  any  toreign 
rraR«mai  ,s  prevent  the  >am»  mi  not  yose  praptry  and  anors 
**&  resuii- 

»  Verify  to*  toe  furschor*  selection  an  toe  ScxpeAAflto  test  topi  or 
oscilloscope  is  ccrreot.  ;.e ,  the  drspiwy  vemca<  resolutcn  «  net 
too  low  or  too  high. 


LIMITED  WARRANTY  &  LIMITATION  OF  LIABILITY 
This  fTuk*  product  w#  be  free  toxn  de*eas  r.  mater.’al  ar<? 
wortransr-C)  for  cp*  year  ton  toe  dart*  cf  ptr-toase  Tnns  warranry 
doat  rvx  ccwr  fuses  dsspesabe  or  carnage  from 

aootJenl  negtect  riww  or  »ocnr*-  oofK#t<ons  of  operation  <x 
handing.  Reseilms  ane  not  authciiMd  to  *jt*no  a^y  oto*'  warranty 
or  Fluke  s  behalf  To  cbuwT  serve*  durog  to*  warranty  p*ro& 
venc  your  >*tectv*  crcduct  to  to*  ofeius:  Ruke  Authaize-c  Servics 
Ooaw  w*  a  ciescr  pboo  of  the  prociem 


TH5  WARRANTY  IS  VSUR  ONLY  REMEDY  MO  OTmER 
WARRAfTlfiS  SUCH  A3  FsTNcSS  FCR  A  ^ARTICULAR 
P  JRPOSE  ARE  EXPRESSED  OR  kfH-.EO  FLu+  S  IS  YCT 
L  ABLE  FCR  ANY  SPECIAL.  lNOiREDT.  NCiDEN'AL 
CCMSEwUENTIAL  DAMAGES  OR  LOSSES  ARISING  FROM  AMY 
CAUSE  CR  THEORY, 

■S<noe  some  states  a  ecunthes  dc  not  atio-w  ttie  exclusion  .y 
Imtat-en  of  an  rp'ied  w»ranry  or  erf  ncidentai  or  oensepwemiar 
damages,  this  lirmaticn  of  vab&ty  may  not  apply  to  you 
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FkAe  imfustra-  B.V 
BoifiC 

~<CC  A3  Aimelc 
The  hesh*hjrs.« 


SERVICE  CENTERS 

"c  roc  at*  an  autovued  service  oerm>-. 
visit  us  on.  the  Work!  Wet*  'A*b 


http.  If  www.  fluke,  ewn 

cr  c».  Fuse  usr»g  ary  c f  the  phone  ryjrrfaefs  isard  heiow 
♦  f  -EES-rJ93-6653  in  J.S  A  and  Canaea 
wT*40-M7-eaflD  in  Europe 
•*■  !-42544f-5CCC  from  other  court  net 


RUGGED  SOLID  STATE  TRANSDUCERS 

WITH  AMPLIFIED  OUTPUTS 

GAGE ,  ABSOLUTE,  AND  COMPOUND  PRESSURES 


tXClUSfVE! 


PX209/PX219  Series 

0-15  to  0-300  psi 
0-1  to  0-20  bor 


Starts  at 


v*  Stainless  Steel 
Fitting  and  Body 
v0  5-Point  NIST-Traceable 
Calibration  included 
Solid  State  Media 
Isolation  (Suitable  for 
Use  with  Most  Industrial 
Fluids  and  Oils) 
v"  Broad  Temperature- 
Compensated  Range 
of  -20  to  80  C 
(-4  to  1 76°F)  Yields  High 
Stability  with  Changing 
Temperatures 
v0  Electrical  Isolation  to 
100  Mq  Ensures  Long- 
Term  Reliability 
v0  Rugged  High  Shock  and 
Vibration  Design  for 
Tough  OEM  Applications 
V  100.000  Hr  MTBF  Typical 

Based  on  proprietary  sensor 
technology  developed  Dy  OMEGA 
to  meet  the  high  relfebtlity  and 
accuracy  demanded  by  military 
applications,  the  PX209^219  Series 
voltage  and  current  output  pressure 
transducer  offers  superior 
performance  in  non-corrosive 
applications,  including: 
engine/powertrafn  testing,  well 
monitoring,  jet  fuel  pressure 
metering,  and  ground  and  race 
water  monitoring.  The  transducer 
uses  a  4-actrve-arm  bridge  sensor 
with  a  micro- machined  diffused 
silicon  diaphragm  and  proprietary 
thin-film  media,  plus  dielectric 
isolation  barriers. 


PIN  WIRE 

1  *  E.XC  RED 

2  COM  BK 

3  .OUT  WT 


Dimensions:  mm  (in) 


4  to  20  mA  Output 
Excitation  24  Vdc  (7  to  35  Vdc) 
reverse  polarity  protected 

Output.  4  to  20  mA  (2-wine)  ±1%  FSO 
Zero  Balance  4  mA  ±2%  FSO 
Max  Loop  Resistance: 

50  x  (supply  voltage  - 10)  Q 

Common  Specifications 
Accuracy:  0.25%  FS  (including 
Hr^aitty,  hysteresis  and  repeatabiMy) 
Operating  Temperature: 

-54  to  12rC  (-65  to  250T) 
Compensated  Temperature: 

-20  to  80’C  /~4  to  176T) 

Thermal  Effect*: 

0.04%  fs/*c  (0.oe%  fs  pf) 

Proof  Pressure:  150% 

Burst  Pressure:  000%  range  max 
Response  Time:  2  ms  typical 


This  same  core  sensing  element 
technology,  which  includes  multiple 
types  of  signal  conditioning  and  the 
ability  to  survive  extremes  of  shock 
and  vibration,  provides  a  modular 
building  block  for  OMEGA’S 
revolutionary  family  of 
pressure-sensing  Instruments. 


SPECIFICATIONS 
Voltage  Output 
Excitation.  24  Vdc  ©  1 5  mA 
5  Vdc  Output:  7  to  35  Vdc 
10  Vdc  Output:  12  to  35  Vdc 
Output  0  to  5  Vdc  or  0  to  10  Vdc 
±  1 .5%  FSO.  3-wtre 

Zero  Balance  0  Vdc  ±2%  FSO 
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VOLTAGE  OUTPUT 
PRESSURE  T  RANSDUCERS 


Appendix  B 


LPAC  System  Diagram 


.P  SERVICE  AIR  SYSTEM 


Appendix  C 


LPAC  Test  Procedure 
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LBES  Non-Intrusive  Load  Monitor  LPAC  Testing 


C.l  Objective 

Observe  the  operation  of  the  Low  Presure  Air  Compressor  (LPAC)  using  the  Non-Intrusive 
Load  Monitor  (NILM)  to  attempt  to  determine  the  cause  of  abnormal  operation  detected 
during  previous  operations. 

C.2  Procedure 

1.  Verify  the  NILM  for  LPAC  #1  is  operational. 

2.  Verify  aggragate  sensor  NILMs  are  operational. 

3.  Verify  the  Low  Pressure  x\ir  System  is  lined  up  to  support  operation  of  the  LPAC  #1. 

4.  Start  LPAC  #1  in  automatic  mode. 

5.  Insert  a  12.5  SCFM  into  the  LP  Air  System  using  the  flow  meter  installed  on  ALP- 
V033. 

6.  Observe  operation  of  the  LPAC  through  a  number  of  load  and  unload  cycles. 

7.  Isolate  shop  air  by  shutting  FA-V007. 

8.  Observe  operation  of  the  LPAC  through  a  number  of  load  and  unload  cycles. 

9.  Open  FA-VO07. 

10.  Secure  LPAC  #1 

11.  Verify  the  NILM  for  LPAC  #2  is  operational. 

12.  Verify  the  Low  Pressure  Air  System  is  lined  up  to  support  operation  of  the  LPAC  #2. 

13.  Start  LPAC  #2  in  automatic  mode. 

14.  Observe  operation  of  the  LPAC  through  a  number  of  load  and  unload  cycles. 

15.  Isolate  shop  air  by  shutting  FA-V007. 
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16.  Observe  operation  of  the  LPAC  through  a  number  of  load  and  unload  cycles. 


17.  Open  FA-V007. 

18.  Secure  LPAC  #2. 
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Appendix  D 


Amplifier  for  ISA  &  3SA  NILM 


This  amplifier  was  designed  by  Professor  Cox  to  amplify  the  voltage  signals  from  the 
Fluke  i3000S.  The  schematic  diagram  for  the  circuits  are  shown  in  Figure  D-l.  The  design 
allows  the  scale  resistors  to  be  changed  out  easily,  allowing  the  amplication  factors  to  be 
changed  as  required.  The  OpAmps  can  also  be  changed  easily  to  any  8-pin  OpAmp  desired. 
The  power  terminals  of  the  amplifier  are  connected  the  corresponding  terminals  on  the  NILM 
power  connection  block.  The  output  terminals  of  the  amplifier  are  connected  to  the  current 
positions  of  SCSI  terminal  board  inside  of  the  NILM  enclosure.  Unlike  a  normal  NILM 
setup,  reference  resistors  are  not  required  for  the  current  measurements  because  the  Fluke 
i3000S  CTs  output  a  voltage  directly.  The  voltage  terminals  of  the  SCSI  terminal  board 
are  connected  to  the  output  terminals  of  the  voltage  sensing  board.  Figure  D-‘2  shows  the 
design  of  the  amplifier  board  in  ExpressPCB  1  M software. 


Screu  Terminal  Block 


Figure  D-l:  Schematic  Diagram  of  the  Amplifier  for  the  ISA  &  3SA  NILM  Systems  usinig 
Fluke  iSOOOS  CTs. 
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Figure  D-‘2:  ExpressPCB™ Image  of  the  Amplifier  for  the  ISA  &  3SA  NILM  Systems  usinig 
Fluke  i3000S  CTs. 
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Appendix  E 


Processing  Scripts  for  NILM  Data 
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These  scripts  were  written  by  previous  researchers  including  LCDR  Bennett  [3],  Profes¬ 
sor  Rob  Cox  [7]  and  Jim  Paris  [?].  They  have  been  modified  as  required  during  this  research 
to  account  for  new  data  sources  and  local  data  storage  and  processing. 


E.l  Reading  Data  DVDs 
E.l.l  ISA 

#! /bin/bash 
for  i  in  *.gz 
do 

echo  $i ; 

gunzip  -c  $i  >  tempfile; 

grep  -v  ~#  tempfile  |  . /chansif t2_12bit .pi  1  2  I  prep  >  lSA_fullprep${i/ . gz/ . txt} 
cat  lSA_fullprep${i/ .gz/ . txt}  I  . /chansif t2_12bit .pi  2  >  lSA${i/ . gz/ . txt} 
rm  -f  $i 
done 

rm  -f  tempfile 


E.l. 2  3SA 

# ! /bin/bash 
for  i  in  *.gz 
do 

echo  $i; 

gunzip  -c  $i  >  tempfile; 

grep  -v  ~#  tempfile  I  . /chansif t2_12bit .pi  1  2  |  prep  >  3SA_fullprep${i/ . gz/ . txt} 
cat  3SA_fullprep${i/ .gz/ . txt}  I  . /chansif t2_12bit .pi  2  >  3SA${i/ . gz/ . txt} 
rm  -f  $i 


done 

rm  -f  tempfile 
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E.1.3  LPAC  #2 


The  script  for  LPAC  #2  requies  different  processing  because  the  data  for  the  pressure 
transducer  is  included  in  the  data  file. 

# ! /bin/bash 

#mount  /cdrom 
#cp  /cdrom/*  . 

for  i  in  * . 1 .gz 
do 

echo  $i; 

gunzip  -c  $i  >  tempfile; 

grep  -v  ~#  tempfile  |  . /chansif t2_12bit .pi  1  2  I  prep  >  LPAC2_f ullprep${i/ . gz/ . txt} 

cat  LPAC2_fullprep${i/.gz/.txt>  I  . /chansif t2_12bit .pi  1  2  >  LPAC2${i/ . gz/ . txt} 
rm  -f  $i 
done 

for  i  in  *.2.gz 
do 

echo  $i; 

gunzip  -c  $i  >  tempfile; 

grep  -v  ~#  tempfile  I  . /chansif t2_12bit .pi  1  >  LPAC2_pressure_${i/ .gz/ . txt> 

rm  -f  $i 

done 

rm  -f  tempfile 
#rm  -f  LPAC* 

#umount  /cdrom 
#eject 

The  \textit{Prune  Script}  reduces  the  size  of  the  pressure  files  by  a  factor  specified  by  t 
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#! /bin/bash 


for  i  in  *.2. txt 
do 

echo  $i; 

cat  $i  |  ./prune. pi  $1  >  ${i/ . 2 . txt/_prune$l . txt} 
done 

#rm  -f  LPAC* 


E.1.4  LPAC  #3 

#! /bin/bash 

#mount  /cdrom 
#cp  /cdrom/*  . 

for  i  in  * . 1 . gz 
do 

echo  $i; 

gunzip  -c  $i  >  tempfile; 

grep  -v  ~#  tempfile  I  . /chansif t2_12bit . pi  1  2  |  prep  >  LPAC3_fullprep${i/ . gz/ . txt} 
cat  LPAC3_fullprep${i/ .gz/ . txt}  |  . /chansif t2_12bit . pi  1  2  >  LPAC3${i/ . gz/ . txt} 
rm  -f  $i 
done 

rm  -f  tempfile 
#rm  -f  LPAC* 

#umount  /cdrom 
#eject 


102 


E.1.5  FOSP 


# ! /bin/bash 

for  i  in  *.gz 
do 

echo  $i; 

gunzip  -c  $i  >  tempfile; 

grep  -v  ~#  tempfile  I  . /chansif t2_12bit .pi  1  2  I  prep  >  F0_f ullprep${i/ .gz/ . txt} 
cat  FO_fullprep${i/.gz/.txt}  |  . /chansif t2_12bit .pi  2  >  F0_Temp${i/ . gz/ . txt} 
grep  -v  ~k  FO_Temp${i/.gz/.txt}  >  F0${i/ . gz/ . txt} 
rm  -f  $i 

rm  -f  FQ_Temp${i/.gz/.txt} 
done 

rm  -f  tempfile 


E.1.6  LOSP 

# ! /bin/bash 

for  i  in  *.gz 
do 

echo  $i; 

gunzip  -c  $i  >  tempfile; 

grep  -v  ~#  tempfile  I  ./chansif t2_12bit .pi  1  2  >  LO_fullprep${i/ .gz/ . txt} 
cat  LO_fullprep${i/.gz/.txt}  |  . /chansif t2_12bit . pi  2  >  L0_Temp${i/ . gz/ . txt} 
grep  -v  ~k  L0_Temp${i/ .gz/ . txt}  >  LO${i/ . gz/ . txt} 
rm  -f  $i 

rm  -f  LO_Temp${i/.gz/.txt} 
done 

#rm  -f  tempfile; 


E.1.7  Panel  Aggregate 


# ! /bin/bash 

#mount  /cdrom 

#cp  /cdrom/*. gz  . 
for  i  in  * . 1 .gz 
do 

if  [  -e  ${i/. l.gz/.2.gz>  ]&&[-©  ${i/ . 1 . gz/ . 3 . gz}  ] 
then 

echo  $i ; 

gunzip  -c  $i  >  tempfilel; 
echo  ${i/. l.gz/.2.gz>; 

gunzip  -c  ${i/. l.gz/.2.gz}  >  tempfile2; 
echo  ${i/.l.gz/.3.gz>; 
gunzip  -c  ${i/. l.gz/.3.gz}  >  tempfile3; 
grep  -v  ~#  tempfilel  >  tempfilela; 
rm  -f  tempfilel; 

grep  -v  ~#  tempfile2  >  tempfile2a; 
rm  -f  tempfile2; 

grep  -v  ~#  tempfile3  >  tempfile3a; 
rm  -f  tempfile3; 

paste  tempfilela  tempfile2a  tempfile3a  >  tempfile; 

cat  tempfile  |  . /chansif t2_12bit . pi  1  2  I  prep  >  UEC2A_f ullprep${i/ . gz/ . txt} 

cat  tempfile  I  . /chansif t2_12bit . pi  1  3  |  prep  >  UEC2B_fullprep${i/ . gz/ . txt} 

cat  tempfile  |  . /chansif t2_12bit .pi  1  4  |  prep  >  SCU_fullprep${i/ .gz/ . txt} 

cat  tempfile  I  . /chansif t2_12bit . pi  1  5  I  prep  >  Panel_fullprep${i/ .gz/ . txt} 

cat  UEC2A_fullprep${i/.gz/.txt}  I  . /chansif t2_12bit . pi  2  >  UEC2A${i/ . gz/ . txt} 
cat  UEC2B_fullprep${i/ .gz/ . txt}  |  . /chansif t2_12bit . pi  2  >  UEC2B${i/ .gz/ . txt} 
cat  SCU_fullprep${i/ .gz/ .txt}  I  . /chansif t2_12bit . pi  2  >  SCU${i/ . gz/ . txt} 
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cat  Panel_fullprep${i/ .gz/ . txt>  I  . / chansif t2_12bit . pi  2  >  Panel_${i/ . gz/ . txt} 

rm  -f  tempfilela; 

rm  -f  tempfile2a; 

rm  -f  tempfile3a; 

rm  -f  tempfile; 

rm  -f  $i 

rm  -£  ${i/.l.gz/.2.gz> 
rm  -f  ${i/.l.gz/.3.gz> 

done 

#umount  /cdrom 
#eject 

E.1.8  GTG 

# ! /bin/bash 

for  i  in  *.gz 
do 

echo  $i; 

gunzip  -c  $i  >  tempfile; 

grep  -v  ~#  tempfile  I  . /chansif t2_12bit .pi  1  2  |  prep  >  GTG_fullprep${i/ . gz/ . txt} 
cat  GTG_fullprep${i/ .gz/ . txt}  I  . /chansif t2_12bit . pi  2  >  GTG${i/ . gz/ . txt} 
rm  -f  $i 
done 

rm  -f  tempfile 

E.1.9  ch  ansi  ft:2  \  2  bi  t  .pi 

This  script  vselects  columns  of  data  in  the  file  as  specified  by  the  user. 

# ! / usr/bin/per 1  -w 

while  (<STDIN>)  { 
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Qcols  =  split; 
foreach  $k  « 3ARGV)  { 

#$var  =  int( ($cols  [$k-l] )) ; 

$var  =  ($cols  [$k-l] ) ; 
print  "$var  \t"; 

> 

print  "\n" ; 

> 

E.2  Matlab  Scripts  to  Plot  Data 

E.2.1  Example  of  Plotting  Data  for  all  Files  for  ISA  NILM 

7.  A  test  script  to  run  through  an  entire  directory 


D=dir ( * * . txt 7 ) ; 

Ptot  =  []  ; 

for  i=l : 1 : length(D) 


V  =  440; 

Rm  =  110; 

K  =  1000; 
g  =  5; 

P  =  [V*prep_to_current (g,  Rm,  K,  load(D(i) . name) )]  ; 
t  =  [0:1: length (P)-l]*( 1/120) *(1/60) ; 

plot (t ,P) ; 

saveas(gcf ,  regexprep(D(i) .name, ' txt ’ , }fi g;)) ; 

Ptot=  [Ptot ;P] ; 

end 

Ttot  =  [0:1: length (Ptot ) -1] *(1/120) * ( 1/60) ; 
plot (Ttot , Ptot) ; 
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saveas (gcf ,  regexprep ( ’ ISATotalPower . f ig J ) ) ; 


E.2.2  Example  of  Time  Aligned  Plot  of  Multiple  Source  Files 

l . 

%  Plot  3SA  Data 

7. . 

x  -  load  (' 3SAsnapshot-20080505- 140001 . 1 . txt ? ) ; 

V  =  440; 

Rm  =  130; 

K  =  1000; 
g  =  1-25; 

P  =  [V*prep_to_current (g,  Rm,  K,  x)] ; 
t  =  [0:1: length (P)~l] *(1/ 120) * (1/60) ; 

h(l)  =  subplot (5, 1 , 1) ; 
plot ( t , ~P ) ; 
title (,UEC2A>); 
xlabeK'Time  (min)  y )  ; 
ylabel ( ; Power  (W) ’ ) ; 


7o  . . 

7,  Plot  L0SP  2A  Data 

V~~~ . ~~~~ 

x  =  load  (’ LDsnapshot-20080505- 140001. l.txtO ; 
V  =  440; 

Rm  =  20; 

K  =  2500; 
g  -  10; 

P  =  [V*prep_to_current (g,  Rm,  K,  x)] ; 
t  =  [0:1: length (P) -!]*(!/ 120) * ( 1/60) ; 
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h(2)  =  subplot (5, 1,2) ; 
plot (t ,-P) ; 
title ('LOSP  2A')  ; 
xlabel ( 'Time  (min) ' ) ; 
ylabel( 'Power  (W) ' ) ; 


% . ~ . 

7o  Plot  FOSP  2A  Data 

°/o . . 

x  =  load  ( 'FOsnapshot-20080505- 140001. l.txt') ; 
V  =  440; 

Rm  =  61; 

K  =  1000; 
g  -  10; 


P  =  [V*prep_to_current (g,  Rm,  K,  -x)]  ; 
t  =  [0:1: length (P) -1] * (1/ 120) * ( 1/60) ; 

h(3)  =  subplot (5 , 1 , 3) ; 
plot (t,P) ; 
title ('FOSP  2A' ) ; 
xlabel ('Time  (min) ') ; 
ylabel ( ' Power  (W) ' ) ; 


% . . . 

I  Plot  #2  LPAC  Data 

l . 

x  =  load  (' LPAC2snapshot-20080505- 14000 1. l.txt') ; 
V  =  440; 

Rm  =  36; 

K  =  2500; 
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g  =  10; 


P  =  [V*prep_to_current (g,  Rm,  K,  -x)] ; 
t  «  [0:1: length (P) -1] * (1/ 120) * (1/60) ; 

h(4)  =  subplot (5, 1,4) ; 
plot (t ,P) ; 
title  ( 'LPAC  #20; 
xlabelOTime  (min) * ) ; 
ylabel ( J  Power  (W) ? ) ; 


°/o  . . 

l  Plot  #3  LPAC  Data 

% . . 

x  =  load  ( ,LPAC3snapshot-20080505- 140001 . 1 . txt O ; 
V  =  440; 

Rm  =  13.5; 

K  =  2500; 
g  -  10; 

P  =  [V*prep_to_current (g,  Rm,  K,  -x)]  ; 
t  -  [0: 1 : length (P)-l]*( 1/120) *(1/60) ; 

h(5)  =  subplot (5, 1,5) ; 
plot (t , P) ; 
title ( 'LPAC  #30; 
xlabelOTime  (min)  O  ; 
ylabel ('Power  (W) O ; 

saveas(gcf  ,  '3SATimeAlignPlot '  ,  }f  igO  ; 
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F.l  Main  NILM  GUI  Program 


function  varargout  =  NILM.GUI (varargin) 


y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /e  /e  /«  /c  ft  ft/tftftft/t  ft  /«  /©  /o  /o  /e  /o  /e  ft  /o  /©  /©  /o  ft  ft  ft  ft  /o  /o  /o  /o  /o  /©  /e  /o  ft  /«  /o  /o  /o  ft  /o  /o  /#  /o  ft  /o  /o  /o  /o  ft  /o  /©  /o  /©  ft  /o  /o  ft  ft  ft  ft  /o  /o  /o  /o  /o  /o  /o  /e  /o  /a  ft 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


X  X 

X  INITIALIZATION  TASKS  OF  THE  GUI  X 


X 


V  7 7  7 7 7 7 7 7  7  7 7  7 7 7  7 7 7 7  7 7 7 7  7 7 7 7 7  7  7  7  7  7  7  7 7  7 7 7 7  7  7  7  7 7  7  7  7 7 7  7  7  7 7 7  7  7 7  7 7 7  7  7  7  7  7 7 7  7  7  7  7 7  7  7 

/©  /o  /o  /o  /o  /o  /o  /«  /o  /o  /c  /o  /O  ft  ft  /©  /©  /o  /o  /o  /©  /o  /o  /o  /o  /o  /o  /©  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  It  It  ft  ft  It  ft  It  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  it  ft  ft 


i:i:itUm:itUmtm:i:^^^^ 


global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS.LIST  EVENT.STATS  PICTURES 
GUI .LAYOUT 


7© 


7©  Create  the  figure,  the  indications  and  the  system  diagram 

7o . . . ~ . . 

DrawSystemDiagramGUI () 

InitializeVariables () 


7o . . . . 

7©  Find  and  stop  all  running  timer  functions  to  prevent  multiples  running 
7©  on  the  system  should  it  be  shutdown  incorrectly 

7© . . . . . . . 

t imer s=t imerf ind ; 
delete (timerfind) ; 


7© . 

7©  Initialize  the  timer  function 

7©  ~  the  timer  will  execute  functions  automatically  to  keep  the  data 
7©  up  to  date 

7© - - - - - - 

t.main  =  timer ( J period , ,2.0, 'ExecutionMode ? , ’fixedrate' , 'StartDelay7 ,0.5) ; 
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t_main.TimerFcn  =  ({OMyTimerFunction}) ; 
start (t_main) ; 


7o  . ~ . . . ~ . . 

%  Initialize  the  alarm  timer  function 

°/9  -  the  timer  will  execute  functions  automatically  to  keep  the  data 
7o  up  to  date 

V~~~ . . . . . ~ . 

t_alarm  ~  timer ( > period , ,4.0, 'ExecutionMode ? , 'fixedrate' , ,StartDelay> ,0.5) ; 
t_alarm.TimerFcn  =  (OAlarmSound) ; 
start (t_alarm) ; 


mmmmmmmnmmm immmmmmmmmmmmiimm 

7o  7. 

%  UTILITY  FUNCTIONS  OF  THE  GUI  l 

%  7. 


0  /  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of  Of 

to  /o  /o  to  to  /o  /o  /o  lo  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to 


7,  Main  GUI  Function 


function  MainGUlQ 


y. . . . 

'/,  Update  the  GUI 

V. . 

UpdateGUIDataO 


% . . . . 

7  Check  for  New  *.evt  Files 

X . . 
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CheckForNewF i les ( ) 


7«  End  Main  GUI  Function 


% . . . . . . . . . 

*/,  Play  the  Alarm  Function 

*/. . 

function  AlarmSound(obj ,  event) 

global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS. LIST  EVENT. STATS  PICTURES  . . . 
GUI .LAYOUT 

if  (ALARMS. SHOW.WARNING) 

sound (PARAMETERS . ALARM.SOUND) ; 

end 

% . . . . — . . . . . . . . 

7,  End  Play  the  Alarm  Function 

7. - - - . . . . . . . . . 


7. . . . . . 

7,  Update  GUI  Data  Function 

7. - - - - - - - . 

function  UpdateGUIDataO 

global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS.LIST  EVENT. STATS  PICTURES  . . . 
GUI .LAYOUT 


7. . 

7.  Show  correct  pump  states  for  running  conditions 

7. . . . 

if  (STATE. VP. ON) 

set (PICTURES. RUNNING. VP. PIC, ’Visible’ ,  ’on’); 
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else 


set (PICTURES . RUNNING . VP . PIC , ' Visible ’ ,  ' of f ' ) ; 

end 


if  (STATE. DP. ON) 

set (PICTURES. RUNNING. DP. PIC, 'Visible’ ,  'on') ; 

else 

set (PICTURES . RUNNING . DP . PIC , ' Visible ’ ,  ' of f ’ ) ; 

end 


% . . . 

%  Show  correct  states  for  alarm  conditions 


if  (ALARMS. PF. ACTIVE) 

set (PICTURES. ALARMS. PF. PIC,  'Visible' ,  'on'); 

set (GUI.LAYOUT. Warning,  'String',  AL ARMS. PF. WARNING, 

' BackgroundColor ’ ,  ' red ’ ) ; 
set(GUI_LAYOUT. ResetButton,  'Visible',  'on'); 

ALARMS . SHOW_WARNING  =  true; 

end 

if  (ALARMS. CSW. ACTIVE) 

set (PICTURES. ALARMS. CSW. PIC,  'Visible',  'on'); 
set (GUI .LAYOUT. Warning,  'String',  ALARMS. CSW. WARNING, 
' BackgroundColor ’ ,  ’ red ' ) ; 
set (GUI.LAYOUT. ResetButton,  'Visible',  'on'); 

ALARMS . SHOW.WARNING  =  true; 

end 

if  (ALARMS. CGL. ACTIVE) 

set (PICTURES. ALARMS. CGL. PIC,  'Visible',  'on'); 
set (GUI .LAYOUT. Warning,  'String',  ALARMS. CGL. WARNING, 
’ BackgroundColor ’ ,  ' red ’ ) ; 
set (GUI.LAYOUT. ResetButton,  'Visible',  'on'); 
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ALARMS . SHOW_WARNING  =  true; 


end 

if  ( ALARMS. FL. ACTIVE) 

set (PICTURES. ALARMS. FL. PIC,  ’Visible’,  ’on’); 

set (GUI.LAYOUT. Warning,  ’String’,  ALARMS. FL. WARNING,  ... 

’ BackgroundColor ’ ,  ’ red ’ ) ; 
set(GUI_LAYOUT.ResetButton,  ’Visible’,  ’on’); 

ALARMS. SHOW. WARNING  =  true; 

end 

if  (ALARMS . SL . ACTIVE) 

set (PICTURES. ALARMS. SL. PIC,  ’Visible’,  ’on’); 

set (GUI .LAYOUT. Warning,  ’String’,  ALARMS. SL. WARNING,  ... 

’ BackgroundColor ’ ,  ’ red ’ ) ; 
set(GUI_LAYOUT.ResetButton,  ’Visible’,  ’on’); 

ALARMS . SHOW.WARNING  =  true; 

end 

if  ("ALARMS. SHOW.WARNING) 

set(GUI_LAYOUT. ResetButton,  ’Visible’,  ’off’); 

end 


% . 

%  Display  the  information  for  the  current  event 

% . . . . 

[m,n]  =  size (EVENTS.LIST) ; 

set (GUI.LAYOUT . DateValue , ’String’ ,  datestr (EVENTS.LIST (n) . DATE.TIME, 1) ) ; 
set (GUI.LAYOUT . TimeValue , ’String’ ,  datestr (EVENTS.LIST (n) . DATE.TIME, 13) ) ; 
set (GUI.LAYOUT. Event Value, ’String’ ,  EVENTS.LIST (n) . EVENT.TEXT) ; 


% . . . . . 

%  Update  the  statistical  information  for  the  pumps 

% . 
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if  (isf ield(EVENT_STATS. VP,  ’ RUN.TIME ’ ) ) 

ART, VP  =  num2str (round (mean ( [EVENT_STATS . VP. RUN.TIME] ) ) ) ; 
set (GUI_LAYOUT . VPARTValue ,  ’String',  strcat (ART, VP,  ’  sec’)); 

end 

if  (isf ield(EVENT_STATS . DP,  ’RUN_TIME’ ) ) 

ART_DP  =  num2str (round (mean ( [EVENT.STATS . DP .RUN_TIME] ) ) ) ; 
set (GUI_LAYOUT . DPARTValue ,  ’String’,  strcat (ART_DP,  ’  sec’)); 

end 


Time24hrs  =  EVENTS, L 1ST (n) . DATE_TIME  -  1; 

% . . . 

*/,  Delete  all  times  in  the  TRACKING .VP . RUNS24  structure  which  are 
*/.  >  24hrs  old 

•/. . . . - . 

if  (max (TRACKING. VP. RUNS24)  >  0) 

OldTimeMax  =  1 ; 

X . . . . - . 

X  Find  the  highest  index  of  the  time  >  24hrs  ago 

while  (TRACKING. VP. RUNS24 (OldTimeMax) <  Time24hrs) 

TRACKING. VP. RUNS24 (OldTimeMax)  <  Time24hrs; 

OldTimeMax  =  OldTimeMax  +  1 ; 

end 


X . . . . . 

X  Delete  the  times  >  24hrs  ago 

x — . — . — 

if  (OldTimeMax  >  1) 

TRACKING . VP . RUNS24 ( 1 : OldTimeMax) = [] ; 

end 
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[m,n]  =  size (TRACKING. VP. RUNS24) ; 

set (GUI_LAYOUT . VP24HValue ,  ’String’,  n) ; 

end 


7. . 

7,  Delete  all  times  in  the  TRACKING .DP . RUNS24  structure  which  are 
%  >  24hrs  old 

7. . . . . . — . 

if  (max (TRACKING. DP. RUNS24)  >  0) 

for  i=l: length (TRACKING. DP. RUNS24) 

if  (EVENTS_LIST(n) . DATE.TIME(i)  <  Time24hrs) 

TRACKING. DP. RUNS24(i)  =  [] ; 

end 

end 

[m,n]  =  size (TRACKING. DP. RUNS24) ; 
set (GUI_LAY0UT . DP24HValue ,  ’ String ’ ,  n) ; 

end 

drawnow 

7,  End  Update  GUI  Data  Function 

7. . — . . . . . . . . . 


7. 

7. 

7. 

7. 

7. 


7. . . . . . 

7.  My  Timer  Function 

function  MyTimerFunction(obj ,  event) 

global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS_LIST  EVENT, STATS  PICTURES  . . . 
GUI .LAYOUT 
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MainGUI () 


v~~ . . . 

70  End  My  Timer  Function 
7o . 


% . ~ . ~ . ~ . . . ~~~~ 

°/e  Function  To  Check  For  New  Data  Files 

v~~ . ~ . 

function  CheckForNewFilesO 

global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS.LIST  EVENT.STATS  PICTURES  . . 
GUI.LAYOUT 


7o . . . ~~~~ . ~ . ~ . 

7  Check  for  new  files  in  directory 

l . ~ . . . ~ . 

FileList  =  dire  (str  cat  (PARAMETERS  .  NEWDATA,  ?/*.evt?)  ,n,,d,>; 

FileListSize  =  size (FileList) ; 

while  (FileListSize(l)  >  0) 

70  If  there  are  data  files  in  the  directory,  read  the  information 

%— . ~~ . . . . . . 

LongFileName  =  s treat (PARAMETERS . NEWDATA,  char (FileList (1 , 1) )) ; 
[DateTimeStamp,  EventCode,  event^class,  window_size,  P,  Q,  local^index] 

Read_File_Inf o (LongFileName) ; 

DateNum  =  datenum(DateTimeStamp,  ^yyymmdd-HHiMM: SS J ) ; 

DateVector  =  datevec (DateNum) ; 
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X  Append  data  to  Event  Data  Matrices 

X . . 

[m,n]  =  size(EVENTS_LIST) ; 

EVENTS_LIST(n+l) .DATE_TIME  =  DateNum; 
EVENTS_LIST(n+l) .TIMESTAMP  =  DateTimeStamp ; 
EVENTS_LIST(n+l) .EVENTCODE  =  str2num(EventCode) ; 
EVENTS_LIST(n+l) .FILENAME  =  char (FileList (1 , 1) ) ; 
EVENTS, LI ST (n+1) .EVENT.TEXT  =  event.class; 


7. . 

7,  Classify  the  Event 

X . . . 

EventClassif ication(n+l) 


7. . . - . . . . 

X  Call  the  Diagnostic  Function 

r~~ . — . 

Diagnostics (n+1) 


7. . . 

7.  Update  the  GUI 

X . 

UpdateGUIDataQ 


X . 

X  move  files  to  Old  Data  directory 

7. . . . 

FN.from  =  strcat (PARAMETERS. NEWDATA,  char (FileList (1, 1))) ; 
FN_to  =  strcat (PARAMETERS . OLDDATA,  char(FileList(l, 1))) ; 
movefile (FN_from,  FN_to) ; 
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•/. 


'/»  Delete  First  Row  of  FileList  Array 

7. . 


FileList (1 ,  :  )=  []  ; 


7.  Recheck  Directory  for  New  Files 

7.-~~ . 

FileList  =  dirc(’NewData/*.evt’ , ’ ’ , ’  d’) ; 
FileListSize  =  size (FileList) ; 


end 


7,  End  of  Check  For  New  Files  Function 


7. 


mmmmmmmxmmmmmmmra ;rammmmrammn 


CLOSING  FUNCTIONS  OF  THE  GUI 


7. 


•/ 7 7  7  7  7  7  7  7  7  7 7  7 7  7  7  7  7  7  7  7  7  7  7 7 7  7  7  7  7 7  7  7 7 7  7  7  7  7 7  7 7  7  7  7 7 7 7  7 7 7  7  7 7 7  7  7 7 7 7 7 7  7 7 7 7 7 7  7  7 7 7 7 7 7 

/o  /o  /e  ft  ft  /o  /o  /o  ft  /o  ft  ft  It  /o  /o  /o  /o  /o  /o  /o  ft/t  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  /o  /o  /o  /c  ft  ft  ft  ft  /«  /#  /o  /«  /o  /o  /o  /e  /o  /e  /o  /o  /o  / t  ft  /o  /o  ft  ft  ft  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  ft  /o  /o  /o 


o  /  y  y  o/  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 
/o  ft  /o  /o  /o  /«  ft  /o  /o  /o  /e  ft  /o  /o  /o  /o  ft  /o  /o  /o  /o  ft  ft  /e  /o  /o  ft  /o  /o  /«  /o  ft  /o  /o  /o  /o  ft  /o  /o  /o  /o  /o  /o  /o  /o  /o  h  / o  /o  /o  /o  /e  /o  /o  /e  /o  /o  /o  /o  /o  h  h  /o  /o  /e  /o  /o  /#  /c  /o  / o  /o  /o  /o  /o 


y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 
/o  /o  /o  /o  h  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  It  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  It  ft  ft  ft  It  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft 


°/e  End  My  Functions 


y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 
ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  It  It  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  It  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft 
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7*  Non- Intrusive  Load  Monitor  (NILM) 

7o  Initialize  Variables  Function 
7,  Written  by:  Richard  Jones 


function  InitializeVariables () 


V  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/a  /o  /a  /o  /a  /o  /o  /o  /o  /o  /o  /o  /o  /a  /o  /o  /o  /o  /a  /a  /o  /a  /o  /a  /a  /o  /a  /o  /o  / o  /o  /o  /o  /a  /o  /o  /a  /a  /o  /o  /o  /o  /o  /o  /a  /o  /o  /a  /o  /a  /o  /c  /e  /a  /o  /o  /e  /o  /o  /o  /o  /c  /o  /o  /a  /o  /o  /a  /a  /«  /a  /a  /o  /o  /a 


7o  Global  Variables 


y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/a  /a  /a  /o  /o  /o  /o  /o  /e  /o  /o  /e  /a  /a  /a  /a  /a  /a  /«  /e  /o  /o  /o  /o  /o /a  /a  /a  /o  /a  /e  /o  /a  /o  /o  /e  /o  /o  /o  /o  /o  /o  /a  /o  /a  /a  /a  /a  /o  to  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a 


global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS.LIST  EVENT_STATS  PICTURES  . . . 
TEMPRUN 


1 . . . 

7,  PARAMETERS 
°/t . . . 

PARAMETERS . OLDDATA  =  ’ OldData/ ’ ; 

PARAMETERS . NEWDATA  =  'NewData/’; 

PARAMETERS.  PICTURES  =  'Pictures/'’; 

PARAMETERS. BACKUP  =  'DataBackup/' ; 

PARAMETERS. SOUNDS  =  ’Sounds/'; 

PARAMETERS. SYSTEM_PICTURE  =  'CHT_SYSTEM_DP_FAULT.jpg' ; 
PARAMETERS .ALARM_SOUND  =  wavread( 'Sounds/DangerNew . wav’ ) ; 
PARAMETERS. INITIALIZE  =  false; 


% . 

%  STATE 

7. . 

STATE. VP. ON  =  0; 
STATE. DP. ON  =  0; 


STATE . VP . START  =  false; 
STATE. VP. STOP  =  false; 
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STATE. DP. START  =  false; 
STATE. DP. STOP  =  false; 
STATE . VP . ERROR  =  false; 
STATE. DP. ERROR  =  false; 


7o . 

t  ALARMS 

°/o . 

ALARMS. CGL. ACTIVE  =  false; 

ALARMS. CGL. WARNING  =  'Clogged  Gauge  Line  Detected'; 

ALARMS. CGL. VISIBLE  =  false; 

AL ARMS. PF. ACTIVE  =  false; 

ALARMS. PF. WARNING  =  'Probe  Failure  Detected'; 

ALARMS. PF. VISIBLE  -  false; 

ALARMS. PF. LIST  =  0; 

ALARMS. SL. ACTIVE  =  false; 

ALARMS. SL. WARNING  =  'Slow  Leak  Detected  in  the  System'; 

ALARMS. SL. VISIBLE  =  false; 

ALARMS. FL. ACTIVE  =  false; 

ALARMS. FL. WARNING  =  'Fast  Leak  Detected  in  the  System'; 

ALARMS. FL. VISIBLE  =  false; 

ALARMS. CSW. ACTIVE  =  false; 

ALARMS. CSW. WARNING  =  'Clogged  Seal  Water  Line  to  the  Vacuum  Pumps  Detected' 
ALARMS. CSW. VISIBLE  =  false; 

ALARMS . SHOW. WARNING  =  false; 


7, . 

7,  TRACKING 

X . 

TRACKING. VP. RUNS24  =  0; 
TRACKING. DP. RUNS24  =  0; 
TRACKING . VP . AVERUNTIME  =  0; 
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TRACKING. DP. AVERUNTIME  -  0; 


r~~ . 

%  EVENTS  LIST 

7e . 

[m,n]  =  size (EVENTS. LI ST) ; 
if  (n  <  2) 

EVENTS.LIST . DATE.TIME  -  0; 
EVENTS.LIST. TIMESTAMP  = 
EVENTS.LIST . EVENTCODE  =  0; 
EVENTS.LIST. FILENAME  =  > , ; 
EVENTS_LIST.EVENT.TEXT  * 

end 


. 

I  EVENT  STATS 

r~~~ . ~~ 

if  ( i sf i e Id (EVENT.STATS ,  ’VP’)) 

[m,n]  =  size (EVENT.STATS . VP) ; 

else 

n  *  0; 

end 

if  (n  <  2) 

EVENT.STATS . VP . START  «  0; 

EVENT. STATS. VP. STOP  =  0; 

EVENT.STATS . VP . START.FILENAME  «  ’  ’  ; 
EVENT.STATS. VP. STOP.FILENAME  =  J , ; 
EVENT.STATS . VP . RUN.TIME  =  0; 
EVENT.STATS . VP . OFF.TIME  =  0; 

end 
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if  (isfield(  EVENT.  ST  ATS ,  'DP')) 

[m,n]  «  size (EVENT_STATS . DP) ; 

else 

n  «  0; 

end 

if  (n  <  2) 

EVENT.STATS . DP . START  =  0; 

EVENT. ST ATS. DP. STOP  =  0; 

EVENT.STATS . DP . START.FILENAME  -  ”  ; 
EVENT.  STATS.  DP.  STOP.FILENAME  =  > 
EVENT.STATS . DP . RUN.TIME  =  0; 

EVENT. STATS. DP. OFF.TIME  =  0; 

end 


7o - 

°/o  PICTURES 

7o . - 

P I CTURES. ALARMS. CGL. SHOW  «  false; 
PICTURES. ALARMS. PF. SHOW  =  false; 

PI CTURES. ALARMS. CSW. SHOW  =  false; 
PICTURES. RUNNING. VP. SHOW  =  false; 
PICTURES. RUNNING. DP. SHOW  =  false; 
for  i=l : 3 

P I CTURES. AL ARMS. SL(i) .SHOW  =  false; 

end 

for  i=l:3 

PICTURES. ALARMS. FL(i) .SHOW  =  false; 

end 


7o . . 

%  TEMPRUN 
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Vo . 

if  ("isfield (TEMPRUN,  >VP>)) 
TEMPRUN . VP . TIME  -  0; 

TEMP RUN . VP . FILE  »  ’ , ; 

end 

if  ("isfield (TEMPRUN,  'DP’)) 
TEMPRUN. DP. TIME  =  0; 
TEMPRUN. DP. FILE  = 

end 


0/ « ,  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 
/e  to  /o  /o  /a  /«  to  /o  /o  /o  /o  to  /o  /o  /o  /o  to  /o  /o  /o  /o  /«  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /a  /a  /o  to  to  to  to  to  to  t a  /a  /o  to  to  to  to  to  to/ofo/o/o/o  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to 

70  End  of  Function 

y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 
to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to  to 
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function  DrawSystemDiagramGUI () 

global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS.LIST  EVENT.STATS  PICTURES  . . . 

VP.ON  VPON  GUI.LAYOUT 

scrsz  =  get (0, ’ScreenSize’ ) ; 

GUIWidth  =  750; 

GUIHeight  -  600; 

GUI.LAYOUT. MainFigure  =  figure(...  °/0  The  main  GUI  figure 

’MenuBar’ , ’none’ ,  ... 

’Toolbar’ , ’none’ ,  ... 

’HandleVisibility’ , ’callback’ ,  ... 

’Color’ ,  get(0, ’ def aultuicontrolbackgroundcolor ’ ) ,  ... 

’Position’,  [ (scrsz (3) -GUIWidth) /2  (scrsz (4) -GUIHeight) /2  GUIWidth 
’Name’ , ’Non- Intrusive  Load  Monitor  (NILM) ’ ,  ... 

’NumberTitle ’ ,  ’off’,  ... 

’HandleVisibility’ , ’on’) ; 

GUI. LAYOUT. Plot Axes  =  axes(...  %  Axes  for  plotting  the  system  diagram 

’Parent’ ,  GUI.LAYOUT . MainFigure ,  . . . 

’Units’,  ’normalized’,  ... 

’HandleVisibility’ , ’callback’ ,  ... 

’Position’ , [0.0  0.01  0.85  0.90],  ... 

’DataAspectRatioMode ’ ,  ’manual’ ,  ... 

’PlotBoxAspectRatio ’ ,  [111],  ... 

’PlotBoxAspectRatioMode ’ ,  ’manual’ ,  ... 

’XTick’  ,  [],  ’  YTick’  ,  []  ,  ... 

’NextPlot’,  ’add’,  ... 

’HandleVisibility’ , ’on’) ; 

GUI.LAYOUT. ResetButton  =  uicontrol (GUI.LAYOUT. MainFigure, ’Style’ ,  ’pushbutton’,  ... 

’String’,  ’Reset’,  ... 

’Position’ ,  [660  480  90  30] ,  ... 

’Callback’ ,  {OAlarmReset. Callback} ,  . . . 

’FontSize’,  12,  ’FontWeight ’ ,  ’bold’,  ... 
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'Visible ' ,  'off'); 


GUI.LAYOUT.CloseButton  =  uicontrol (GUI_LAYOUT. MainFigure , 'Style' ,  'pushbutton', 

'String',  'Quit  NILM',  . . . 

'Position' ,  [660  30  90  30] ,  ... 

'Callback' ,  {@ExitButton_Callback} ,  . , . 

'FontSize',  12,  'FontWeight ' ,  'bold',  ... 

'Value',  1); 

GUI_LAY0UT. VPARTTitle  =  uicontrol(  ...  °/0  Text  Title  for  VP  Ave  Run  Time 

GUI_LAYQUT . MainFigure ,  ... 

'Style ' ,  'text ' ,  ... 

'String',  'VP  Ave  Run  Time',  ... 

'Units',  'Normalized',  ... 

'Position',  [0.80,  0.65,  0.2,  0.04],  ... 
'BackgroundColor ' ,  'blue',  ... 

'ForegroundColor ' ,  'yellow' ,  ... 

'FontSize',  12,  'FontWeight',  'bold'); 

GUI_LAY0UT .  VPARTValue  =  uicontrol(  ...  °/0  Text  Title  for  VP  Ave  Run  Time 

GUI_LAY0UT . MainFigure ,  ... 

'Style' ,  'text' ,  ... 

'String' ,  'VP  ART' ,  ... 

'Units',  'Normalized',  ... 

'Position',  [0.80,  0.61,  0.2,  0.04],  ... 
'BackgroundColor',  'white',  ... 

'FontSize',  12); 

GUI.LAY0UT. VP24HTitle  =  uicontroK  ...  70  Text  Title  for  VP  Ave  Run  Time 

GUI  ..LAYOUT .  MainFigure ,  ... 

'Style' ,  'text' ,  ... 

'String',  'last  24hr  VP  Runs',  ... 

'Units',  'Normalized',  ... 

'Position',  [0.80,  0.55,  0.2,  0.04],  ... 
'BackgroundColor' ,  'blue' ,  ... 

'ForegroundColor' ,  'yellow' ,  ... 
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’FontSize’,  12,  ’FontWeight ’ ,  'bold'); 


GUI.LAYOUT . VP24HValue  =  uicontrol(  ...  %  Text  Title  for  VP  Ave  Run  Time 

GUI^LAYOUT . MainFigure ,  ... 

’Style’ ,  ’text’ ,  ... 

’String’ ,  ’VP  24h’ ,  . . . 

’Units’ ,  ’Normalized’ ,  ... 

’Position’,  [0.80,  0.51,  0.2,  0.04],  ... 
’BackgroundColor ’ ,  ’white’,  ... 

’FontSize’,  12); 

GUI_LAYQUT.DPARTTitle  =  uicontrol(  ...  °/0  Text  Title  for  VP  Ave  Run  Time 

GUI.LAYOUT . MainFigure ,  . . . 

’Style’ ,  ’text’ ,  ... 

’String’,  ’DP  Ave  Run  Time’,  ... 

’Units’,  ’Normalized’,  ... 

’Position’,  [0.80,  0.35,  0.2,  0.04],  ... 
’BackgroundColor’ ,  ’blue’ ,  ... 

’ForegroundColor ’ ,  ’yellow’ ,  ... 

’FontSize’,  12,  ’FontWeight ’ ,  ’bold’); 
GUI_LAY0UT.DPARTValue  =  uicontrol(  ...  %  Text  Title  for  VP  Ave  Run  Time 

GUI.LAYOUT . MainFigure ,  ... 

’Style’ ,  ’text’ ,  ... 

’String’ ,  ’DP  ART’ ,  . . . 

’Units’,  ’Normalized’,  ... 

’BackgroundColor’ ,  ’white’ ,  ... 

’Position’,  [0.80,  0.31,  0.2,  0.04],  ... 
’FontSize’,  12); 

GUI_LAY0UT . DP24HTitle  =  uicontrol(  ...  %  Text  Title  for  VP  Ave  Run  Time 

GUI.LAYOUT. MainFigure,  ... 

’Style’ ,  ’text’ ,  ... 

’String’,  ’last  24hr  DP  Runs’,  ... 

’Units’,  ’Normalized’,  ... 

’ Position ’ ,  [0.80,  0.25,  0.2,  0 . 04] ,  ... 


'BackgroundColor 7  ,  'blue7,  ... 

'ForegroundColor 7 ,  'yellow7 ,  ... 

'FontSize7,  12,  'FontWeight  7 ,  'bold'); 
GUI_LAY0UT.DP24HValue  =  uicontrol(  ...  %  Text  Title  for  VP  Ave  Run  Time 

GUI_LAYOUT . MainFigure ,  ... 

'Style7 ,  7 text 7 ,  ... 

'String7 ,  7 DP  24h7 ,  ... 

'Units7,  'Normalized7,  ... 

'Position7,  [0.80,  0.21,  0.2,  0.04],  ... 
'BackgroundColor 7 ,  'white7 ,  ... 

'FontSize7,  12); 

GUI_LAY0UT.DateTitle  =  uicontrol(  ...  °/0  Text  Title  for  VP  Ave  Run  Time 

GUI_LAY0UT. MainFigure ,  . . . 

'Style 7 ,  7  text 7 ,  ... 

'String7 ,  'Date : 7 ,  ... 

'Units',  'Normalized7,  ... 

'Position7,  [0.00,  0.0,  0.10,  0.04],  ... 
'BackgroundColor 7 ,  'blue7 ,  ... 

'ForegroundColor 7 ,  'yellow7 ,  ... 

'FontSize7,  12,  'FontWeight 7 ,  'bold'); 
GUI_LAY0UT.DateValue  =  uicontrol(  ...  %  Text  Title  for  VP  Ave  Run  Time 

GUI_LAY0UT. MainFigure,  ... 

'Style 7 ,  7  text 7 ,  ... 

'String7,  'Date  Here7,  ... 

'Units',  'Normalized7,  ... 

'Position7,  [0.10,  0.0,  0.15,  0.04],  ... 
'BackgroundColor 7 ,  'white7 ,  ... 

'FontSize7,  12); 

GUI_LAY0UT.TimeTitle  =  uicontrol(  ...  70  Text  Title  for  VP  Ave  Run  Time 

GUI_LAYQUT . MainFigure ,  . . . 

'Style7 ,  'text7 ,  ... 

'String7 ,  'Time: 7 ,  ... 
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'Units’,  'Normalized',  ... 

'Position',  [0.25,  0.0,  0.10,  0.04],  ... 
'BackgroundColor ' ,  'blue' ,  ... 

'ForegroundColor ' ,  'yellow' ,  ... 

'FontSize',  12,  'FontWeight ' ,  'bold'); 
GUI.LAYOUT.TimeValue  =  uicontrol(  ...  %  Text  Title  for  VP  Ave  Run  Time 

GUI.LAY0UT . MainFigure ,  . . . 

'Style' ,  'text' ,  ... 

'String',  'Time  Here',  ... 

'Units',  'Normalized',  ... 

'Position',  [0.35,  0.0,  0.15,  0.04],  ... 
'BackgroundColor',  'white',  ... 

'FontSize',  12); 

GUI_LAY0UT.EventTitle  =  uicontrol(  ...  °/0  Text  Title  for  VP  Ave  Run  Time 

GUI.LAYOUT . MainFigure ,  ... 

'Style ' ,  'text ' ,  ... 

'String',  'Event:',  ... 

'Units',  'Normalized',  ... 

'Position',  [0.50,  0.0,  0.10,  0.04],  ... 
'BackgroundColor',  'blue',  ... 

'ForegroundColor ' ,  'yellow' ,  ... 

'FontSize',  12,  'FontWeight',  'bold'); 
GUI.LAYOUT.EventValue  -  uicontrol(  ...  °/0  Text  Title  for  VP  Ave  Run  Time 

GUI.LAYOUT . MainFigure ,  . . . 

'Style' ,  'text' ,  ... 

'String',  'Event  Here',  ... 

'Units',  'Normalized',  ... 

'Position',  [0.60,  0.0,  0.45,  0.04],  ... 
'BackgroundColor' ,  'white' ,  ... 

'FontSize',  12); 

GUI_LAY0UT. Warning  =  ui control (  ...  °/0  Warning  text 

GUI .LAYOUT . MainFigure ,  . . . 
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'Style 7 ,  7 text 7 ,  ... 

'String7 ,  77 ,  ... 

'Units7,  'Normalized7,  ... 

'Position7,  [0.0,  0.85,  1.0,  0.05],  ... 
'BackgroundColor 7 ,  'blue7,  ... 
'ForegroundColor 7 ,  'white7,  ... 

'FontSize7,  16,  'FontWeight 7 ,  'bold'); 
GUI_LAY0UT . NILMTitle  =  uicontrol(  ...  %  NILM  Title  text 

GUI.LAYOUT . MainFigure ,  . . . 

'Style7 ,  'text 7 ,  ... 

'String7,  7Non-Intrusive  Load  Monitoring7, 
'Units',  'Normalized7,  ... 

'Position7,  [0.0,  0.90,  1.0,  0.10],  ... 
'BackgroundColor 7 ,  'blue7 ,  ... 
'ForegroundColor 7 ,  'white7 ,  ... 

'FontSize7,  18,  'FontWeight7,  'bold'); 


1 . . . . . . . 

%  Create  object  groups  for  controlling  group  settings  in  the  GUI 

7o . . . . . 

GUI.LAYOUT. AlarmGroup  =  hggroup; 

GUI.LAYOUT. RunningGroup  -  hggroup; 

hold (GUI.LAYOUT . Plot Axes ,  7  on 7 ) ; 

axis (GUI. LAYOUT. Plot Axes,  [0,  800,  0,  600]); 

axis (GUI.LAYOUT . Plot Axes ,  7  of f  7 )  ; 


7o 
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7.  Update  the  System  Diagram 
7. . . . — 


7. . 

X  Draw  Border 

X . 

‘/.rectangle (’ Position ’,  [0,0,  800,  600]); 


X . - . 

X  Draw  all  Circular  objects 

X . 

X  Vacuum  Pump 
AA(1) .x  =  635; 

AA(1) .y  =  400; 

AA(1) . r  =  50; 

AA(1). color  =  ’red’; 

7.  Discharge  Pump 
AA(2) .x  =  635; 

AA(2) .y  =  200; 

AA(2) .r  =  50; 

AA (2). color  =  ’green’; 

7.  Gauge  Body 
AA(3) .x  =  50; 

AA(3) .y  =  400; 

AA(3) .r  =  30; 

AA(3). color  =  ’green’; 

[m,n]  =  size(AA); 

for  i  =  l:n 
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OpenCircleGUI (AA(i) . x,  AA(i).y,  AA(i).r); 
drawnow 


end 


X . — 

X  Alarm  Conditions 

7. . 

%  Vacuum  Pump 
AA(1) .x  =  635; 

AA(1) .y  =  400; 

AA(1) .r  =  50; 

AA(1). color  =  'red’; 

PICTURES. ALARMS. CSW. PIC  =  FillCircleGUKAA(l) .x,  AA(l).y,  AA(l).r,  AA(1) . color) ; 
set (PICTURES. ALARMS. CSW. PIC,  ’Visible’,  ’off’,  ’Parent’,  GUI_LAYOUT. AlarmGroup) ; 
drawnow 

X  Discharge  Pump 
AA(2) .x  =  635; 

AA(2) . y  =  200; 

AA(2) .r  =  50; 

AA(2). color  =  ’yellow’; 

XPICTURES. ALARMS. CSW. PIC  =  FillCircleGUI (AA(2) .x,  AA(2) .y,  AA(2).r,  AA(2) . color) ; 
Xset (PICTURES. ALARMS. CSW. PIC,  ’Visible’,  ’off’,  ’Parent’,  GUI_LAYOUT. AlarmGroup) ; 
drawnow 

X  Gauge  Body 
AA(3) .x  =  50; 

AA (3) .y  =  400; 

AA(3) . r  =  30; 

AA(3). color  =  ’red’; 

PICTURES. ALARMS. CGL. PIC  =  FillCircleGUI (AA(3) .x,  AA(3).y,  AA(3).r,  AA(3) . color) ; 
set (PICTURES. ALARMS. CGL. PIC,  ’Visible’,  ’off’,  ’Parent’ ,  GUI.LAYOUT. AlarmGroup) ; 


drawnow 


7. . 

%  Running  Conditions 

7. . 

7.  Vacuum  Pump 
AA(1) .x  =  635; 

AA(1) .y  =  400; 

AA(1) .r  =  50; 

AA(1). color  =  'green'; 

PICTURES. RUNNING. VP.PIC  =  FillCircleGUKAA(l) .x,  AA(l).y,  AA(l).r,  AA(1) . color) ; 
set (PICTURES. RUNNING. VP. PIC,  'Visible',  'off',  ... 

'Parent',  GUI_LAY0UT. RunningGroup) ; 

drawnow 


7,  Discharge  Pump 
AA(2) .x  =  635; 

AA(2) . y  =  200; 

AA(2) . r  =  50; 

AA(2). color  =  'green'; 

PICTURES . RUNNING . DP . PIC  =  FillCircleGUI (AA (2) . x ,  AA(2).y,  AA(2) .r,  AA(2) . color) ; 
set (PICTURES. RUNNING. DP. PIC,  'Visible',  'off',  ... 

'Parent ' ,  GUI_LAYOUT . RunningGroup) ; 


drawnow 


7. . 

7.  Vacuum  Receiver 
7. . 

7.  Normal  Vacuum  Receiver 

rectangle ( 'Position' , [100, 100,300,400] , 'Curvature ' , [.3] ,  'LineWidth' ,  4) ; 
drawnow 


%  Alarm  Condition  for  Vacuum  Receiver 

PICTURES. ALARMS. PF. PIC  =  rectangle ( ’Position’ ,  [100, 100 , 300,400] ,’ Curvature ’ ,  . 

[.3],  ’LineWidth’,  4,  ’FaceColor’ , ’r’) ; 

set (PICTURES. ALARMS. PF. PIC,  ’Visible’,  ’off’,  ’Parent’,  GUI_LAY0UT . AlarmGroup , 

’ButtonDownFcn’ ,  <2PF_ButtonDownFcn) ; 


drawnow 


7. . 

7.  Gauge 

7. . 

7.  Gauge  Line 
X  =  [80,  100] ; 

Y  =  [400,  400] ; 

line (X ,  Y,  ’Color’,  ’black’,  ’LineWidth’,  8); 


7.  Gauge  Needle 
X  =  [35 ,  65]  ; 

Y  =  [415,  385] ; 

line(X,  Y,  ’Color’,  ’black’,  ’LineWidth’,  2); 

X  =  [35,  35] ; 

Y  =  [405,  415] ; 

line(X,  Y,  ’Color’,  ’black’,  'LineWidth’,  2); 

X  =  [35,  45] ; 

Y  =  [415,  415] ; 

line (X,  Y,  ’Color’,  ’black’,  ’LineWidth’,  2); 


7.  . 

7.  Level  Probes 

7. . 

X  =  [275,  275]  ; 
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Y  =  [150,  500] ; 

line(X,  Y,  ’Color’,  ’blue’,  ’LineWidth’,  3); 


X  =  [250 ,  250] ; 

Y  =  [250,  500] ; 

line(X,  Y,  ’Color’,  ’blue’,  ’LineWidth’,  3); 

X  =  [225,  225]; 

Y  =  [350,  500] ; 

line(X,  Y,  ’Color’,  ’blue’,  ’LineWidth’,  3); 


l . 

’/,  Fast  Leak 

X . 

X_FL(l).xl  =  100; 

X_FL(1) . x2  =  75; 

Y_FL(1) . yl  =  550; 

Y_FL(1) . y2  =  575; 

X_FL(2) . x 1  =  100; 

X_FL(2) . x2  =  100; 

Y_FL (2) .yl  =  550; 

Y_FL (2) . y2  =  585; 

X_FL(3) . x 1  =  100; 

X_FL(3) . x2  =  125; 

Y_FL(3) . yl  =  550; 

Y_FL(3) . y2  =  575; 

for  i=l:3 

PICTURES. ALARMS. FL(i) .PIC  =  line(  [X_FL(i) .xl,  X_FL(i).x2], 

[Y_FL(i).yl,  Y_FL(i).y2], 
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'Color’,  'red',  'LineWidth',  4); 
set(PICTURES. ALARMS. FL(i) .PIC,  'Visible',  'off',  ... 

’ Parent ’ ,  GUI_LAYOUT . AlarmGroup) ; 

drawnow 

end 

7. . 

7.  Slow  Leak 

7, . 

X_SL(1) . x 1  =  150; 

X_SL(1) . x2  =  135; 

Y_SL(1) .yl  =  550; 

Y_SL(1) . y2  =  565; 

X_SL(2) . xl  =  150; 

X_SL(2) . x2  =  150; 

Y_SL(2) .yl  =  550; 

Y_SL(2) . y2  =  575; 

X_SL(3) .xl  =  150; 

X_SL(3) . x2  =  165; 

Y_SL(3) .yl  =  550; 

Y_SL(3).y2  =  565; 

for  i=l:3 

PICTURES. ALARMS. SL(i). PIC  =  line([X_SL(i).xl,  X_SL(i).x2],  ... 

[Y_SL(i).yl,  Y_SL(i).y2],  ... 
'Color',  'red',  'LineWidth',  4); 
set (PICTURES. ALARMS. SL(i) .PIC,  'Visible',  'off',  ... 

’ Parent ’ ,  GUI_LAY0UT . AlarmGroup) ; 

end 
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'/. . 

*/,  Draw  in  Piping  System 

% . 

7,  Vacuum  Piping 
X  =  [50,  550] ; 

Y  =  [550,  550] ; 

line(X,  Y,  'Color’,  ’black’ 

X  =  [250 ,  250] ; 

Y  =  [500,  550] ; 

line(X,  Y,  ’Color’,  ’black’ 

X  =  [550,  550] ; 

Y  =  [550,  400] ; 

line(X,  Y,  ’Color’,  ’black’ 

X  =  [550 ,  635] ; 

Y  =  [400,  400] ; 

line(X,  Y,  ’Color’,  ’black’, 

X  =  [630,  700] ; 

Y  =  [450,  450] ; 

line(X,  Y,  ’Color’,  ’black’, 

7,  Discharge  Piping 
X  =  [250,  550] ; 

Y  =  [50,  50] ; 

line(X,  Y,  ’Color’,  ’black’, 

X  =  [250,  250]; 

Y  =  [50,  100]  ; 

line(X,  Y,  ’Color’,  ’black’. 


’LineWidth’,  4) 


’LineWidth’,  4) 


’LineWidth’,  4) 


’LineWidth’,  4) 


’LineWidth’,  4) 


’LineWidth’,  4) 


’LineWidth’,  4) 
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X  =  [550,  550] ; 

Y  =  [50,  200] ; 

line (X ,  Y,  'Color’,  'black',  'LineWidth' ,  4); 

X  =  [550,  635] ; 

Y  =  [200,  200] ; 

line(X,  Y,  'Color',  'black',  'LineWidth',  4); 

X  =  [630,  700] ; 

Y  =  [250,  250] ; 

line(X,  Y,  'Color',  'black',  'LineWidth',  4); 

function  h  =  FillCircleGUI (x ,y , r ,  color) 
global  GUI.LAYOUT 
hold  on 

th  =  0 : pi/50 : 2*pi ; 
xunit  =  r  *  cos(th)  +  x; 
yunit  =  r  *  sin(th)  +  y; 
h  =  fill (xunit,  yunit,  color,  . . . 
'LineWidth',  4,  ... 

'Parent',  GUI.LAYOUT . PlotAxes) ; 
hold  off 


function  h  =  OpenCircleGUI (x,y,r) 
global  GUI.LAYOUT 
hold  on 

th  =  0 :pi/50 : 2*pi ; 
xunit  =  r  *  cos(th)  +  x; 
yunit  =  r  *  sin(th)  +  y; 

h  =  plot (xunit,  yunit , 'Color ' ,  'black',  'LineWidth',  4,  'Parent',  GUI_LAY0UT. Plot Axes) 
hold  off 
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t  Alarm  Reset  Function 

function  AlarmReset_Callback(hObject ,  eventdata) 

global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS.LIST  EVENT. STATS  PICTURES 
VP_ON  VPON  GUI.LAYOUT  ALARM.NOW 


%  Clear  All  Active  Alarms 

ALARMS. PF. ACTIVE  =  false; 
ALARMS . CSW . ACTIVE  =  false; 
ALARMS. CGL. ACTIVE  -  false; 
ALARMS. FL. ACTIVE  =  false; 
ALARMS. SL. ACTIVE  =  false; 


70  Set  All  Alarm  Pictures  to  'off' 

set (PICTURES. ALARMS. PF. PIC,  'Visible',  'off'); 
set (PICTURES. ALARMS. CSW. PIC,  'Visible',  'off'); 
set (PICTURES. ALARMS. CGL. PIC,  'Visible',  'off'); 

for  i=l:3 

set (PICTURES. ALARMS. FL(i) .PIC,  'Visible',  'off'); 
set (PICTURES. ALARMS. SL(i) .PIC,  'Visible',  'off'); 

end 


%  Reset  the  Alarm  Banner  back  to  original  blue  w/o  text 


r~~ 
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set (GUI_LAYOUT . Warning ,  } String J , 


?  > 

>  ... 

J  BackgroundColor  * ,  J  blue } ) ; 


°/0  Turn  of  alarm  sound 

ALARMS. SHOW„WARNING  «  false; 

°/e  Hide  Reset  Button 

set (GUI.LAYOUT . ResetButton ,  > Visible ’ ,  > of f  ’ ) ; 


%  Exit  Function 

function  ExitButton_Callback(hObject ,  eventdata) 

%  hObject  handle  to  ExitButton  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 
%  handles  structure  with  handles  and  user  data  (see  GUIDATA) 


%  Backup  all  data  prior  to  closing 
t  BackupData 

%  Find  and  stop  all  running  timer  functions 

t imers=timerf ind ; 
delete (timerfind) ; 
close  all; 
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7.  Probe  Failure  Button  Down  Function 

7, . 

function  PF_ButtonDownFcn(hObject,  eventdata) 

global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS.LIST  EVENT. STATS  PICTURES  . . 

VP_ON  VPON  GUI.LAYOUT  ALARM.NOW 

% . 

7«  Clear  PF  Active  Alarms 

% . 

ALARMS . PF . ACTIVE  =  false; 


% . . . . 

7,  Set  PF  Alarm  Pictures  to  'off' 

7. . 

7.  set  (PICTURES.  ALARMS.  PF.  PIC,  'Visible’,  'off'); 


7. . . - - - - 

7.  Reset  the  Alarm  Banner  back  to  original  blue  w/o  text 

7. . . . . . . . 

set (GUI.LAYOUT. Warning,  'String',  ”,  ... 

’ BackgroundColor ' ,  ’ blue ’ ) ; 


7, . 

7.  Turn  of  alarm  sound 

7. . — . 

ALARMS. SHOW.WARNING  =  false; 


7. . . . 

7.  Open  the  Alarm  Information  GUI 

7, . 

ALARMS . ALARMCODE  =  'PF'; 


DrawAlarmGUI () 


7. . 

7,  Clogged  Gauge  Line  Button  Down  Function 

7. . 

function  CGL_ButtonDownFcn(hObject,  eventdata) 

global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS.LIST  EVENT.STATS  PICTURES  . . . 
VP_0N  VPON  GUI.LAYOUT  ALARM.NOW 

7. . 

7.  Clear  CGL  Active  Alarms 

7. . 

ALARMS. CGL. ACTIVE  =  false; 


7. . 

7.  Set  PF  Alarm  Pictures  to  ’off’ 

7. - - - - - 

7.  set  (PICTURES.  ALARMS.  CGL.  PIC,  ’Visible',  'off'); 


7. . . . 

7t  Reset  the  Alarm  Banner  back  to  original  blue  w/o  text 

7. . . . . . 

set (GUI_LAYOUT. Warning,  'String',  ”,  ... 

’ BackgroundColor ’ ,  ’ blue ' ) ; 


7.  Turn  of  alarm  sound 

7. . 

ALARMS. SHOW.WARNING  =  false; 


r~ . - . - 

%  Open  the  Alarm  Information  GUI 

X . 
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ALARMS . ALARMCODE  =  ' CGL ' ; 

DrawAlarmGUK) 

70  Clogged  Seal  Water  Button  Down  Function 

X . . . . . . 

function  CSW_ButtonDownFcn(hObject ,  eventdata) 

global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS, LIST  EVENT, STATS  PICTURES  . 
VP.ON  VPON  GUI.LAYOUT  ALARM, NOW 

°/o . ~ . ~~~~ 

°/e  Clear  CSW  Active  Alarms 

ALARMS. CSW. ACTIVE  =  false; 

70  Set  CSW  Alarm  Pictures  to  'off' 

°/o  set  (PICTURES.  ALARMS.  CSW.  PIC,  'Visible',  'off'); 

% . . . . * . 

°/0  Reset  the  Alarm  Banner  back  to  original  blue  w/o  text 

set (GUI.LAYOUT. Warning,  'String',  ",  ... 

' BackgroundColor ' ,  ' blue ' ) ; 

X . 

70  Turn  of  alarm  sound 

x - 

ALARMS. SHOW.WARNING  =  false; 


X . . . 

7o  Open  the  Alarm  Information  GUI 
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r~ . — . 

ALARMS . ALARMCODE  =  ’CSW’; 
DrawAlarmGUI () 


7. . . . - . . . . . . . 

7.  Fast  Leak  Button  Down  Function 

7. . . . 

function  FL_ButtonDownFcn(hObject,  eventdata) 

global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS_LIST  EVENT_ STATS  PICTURES  . . . 
VP_0N  VPON  GUI.LAYOUT  ALARM_NOW 

7. . . 

7.  Clear  FL  Active  Alarms 

7. . . . . 

ALARMS. FL. ACTIVE  =  false; 


7. . . . - 

70  Set  FL  Alarm  Pictures  to  ’off’ 

7. . — . 

7.  set  (PICTURES.  ALARMS.  FL.  PIC,  ’Visible’,  ’off’); 


7. . . . . . . 

7.  Reset  the  Alarm  Banner  back  to  original  blue  w/o  text 

7. - - - - - - - 

set (GUI.LAYOUT . Warning ,  ’ String ’ ,  ”,  ... 

’ BackgroundColor ’ ,  ’ blue ’ ) ; 


7,  Turn  of  alarm  sound 

7. . . . 

ALARMS. SHOW.WARNING  =  false; 
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’/o  Open  the  Alarm  Information  GUI 

7. . 

ALARMS. ALARMCODE  =  ’FL’; 
DrawAlarmGUI () 


7. . . . - . . . 

7»  Slow  Leak  Button  Down  Function 

7. . . . — . . . . . . 

function  SL_ButtonDownFcn(hObject ,  eventdata) 

global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS.LIST  EVENT, STATS  PICTURES  . . . 
VP_ON  VPON  GUI, LAYOUT  ALARM, NOW 

7. . . 

7.  Clear  SL  Active  Alarms 

7. . . - 

ALARMS . SL . ACTIVE  =  false; 


7. . . . 

7o  Set  SL  Alarm  Pictures  to  ’off’ 

% . 

7.  set  (PICTURES.  ALARMS.  SL.  PIC,  'Visible’,  ’off’); 


7. . 

7,  Reset  the  Alarm  Banner  back  to  original  blue  w/o  text 

7, . . . - 

set (GUI_LAYOUT . Warning ,  ’ String ’ ,  ”,  ... 

’BackgroundColor ’ ,  ’blue’); 


7. . . . 

7  Turn  of  alarm  sound 

7. . 

ALARMS. SHOW.WARNING  =  false; 
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%  Open  the  Alarm  Information  GUI 

% . 

ALARMS . ALARMCODE  =  ’SL’ ; 

Dr awAlarmGUI ( ) 
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°/e  Non-Intrusive  Load  Monitor  (NILM) 

%  Event  Classification  Function 
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70  Written  by:  Richard  Jones 


7  7  7  7  7 7  7  7  7  7  7 7  7  7  7  7  7 7 7  7  7  7  7 7  7 7  7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7 

/o  /o  /o  /a  /a  /a  /o  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /a  /a  /a 


function  EventClassif ication(EventNum) 


‘/a  Global  Variables 


7  7  7  7  7  7  7  7  7  7  7 7  7  7 7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7 

/a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /c  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a 


global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS _ LIST  EVENT.STATS  PICTURES 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%% 


7o  Local  Variables 


7  7  7 7  7  7  7  7  7  7  7 7  7 7 7  7  7  7  7  7  7  7  7  7  7 7  7 7  7  7  7  7  7  7  7  7 7 7  7 7 7  7  7  7  7 7  7  7  7 7  7  7 7 7  7  7  7 7 7 7 7  7  7 7 7 7 7  7  7  7 7  7 7 7 7 

/o  /o  /a  /a  /a  /a  /a  /a  /a  /o  /o  /a  /a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /o  /o  /a  /a  /a  /o  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /a  /a  /o  /«  /a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  / o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /a  /a 


STATE. VP. ERROR  =  false; 


STATE. DP. ERROR  -  false; 


°/a . ~ . — 

7e  Event  Code  Functions 

switch  (EVENTS.LIST (EventNum) .EVENTCODE) 
case  0  l  Unclassif iable 

case  1  l  VP  On 

STATE. VP. ON  =  STATE. VP. ON  +  1; 
STATE . VP . START  =  true; 

PICTURES. RUNNING. VP. SHOW  =  true; 
case  2  l  VP  Off 

STATE. VP. STOP  =  true; 

PICTURES. RUNNING. VP. SHOW  =  false; 
case  3  %  DP  On 


149 


STATE. DP. ON  =  ST ATE. DP. ON  +  1; 

STATE . DP . START  =  true; 

PICTURES. RUNNING. DP. SHOW  =  true; 
case  4  7,  DP  Off 

STATE . DP . STOP  =  true; 

PICTURES. RUNNING. DP. SHOW  =  false; 
case  5  70  VP  On  (Clogged) 

STATE. VP. ON  =  STATE. VP. ON  +  1; 

ST ATE. VP. START  -  true; 

PICTURES. RUNNING. VP. SHOW  =  true; 

STATE. VP, ERROR  =  true; 
case  7  70  Unclassif iable  On  Event 

case  9  7«  All  Off  Event 

ST ATE. VP. STOP  -  true; 

STATE. DP. STOP  -  true; 

PICTURES. RUNNING. VP. SHOW  =  false; 

PICTURES. RUNNING. DP. SHOW  »  false; 
case  11  70  2  VPs  On 

STATE. VP. ON  -  STATE. VP. ON  +  2; 

STATE . VP . START  =  true; 

PICTURES. RUNNING. VP. SHOW  =  true; 
case  13  7  VP  on  and  DP  on 

STATE. VP. ON  =  STATE. VP. ON  +  1; 

STATE . VP . START  «  true; 

STATE. DP. ON  =  STATE. DP. ON  +  1; 

STATE . DP . START  =  true; 

PICTURES. RUNNING. VP. SHOW  =  true; 

PICTURES. RUNNING. DP. SHOW  -  true; 

case  21  %  VP  On  (Normal)  followed  by  a  it  immediately  turning  off 

STATE. VP. ON  =  STATE. VP. ON  +  1; 

STATE . VP . START  =  true; 

STATE . VP . STOP  =  true; 
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STATE. VP. ERROR  -  true; 


case  43  %  DP  On  followed  by  it  immediately  turning  off 

STATE. DP. ON  -  STATE. DP. ON  +  1; 

STATE . DP . START  =  true; 

STATE . DP . STOP  =  true; 

STATE . DP . ERROR  =  true; 

case  65  70  VP  On  (Clogged)  followed  by  it  turning  immediately  off 

STATE. VP. ON  =  ST ATE. VP. ON  +  1; 

STATE . VP . START  =  true; 

STATE . VP . STOP  =  true; 

STATE . VP . ERROR  =  true; 

case  87  %  Unclassif iable  On  Event  followed  by  an  Off  Event 

otherwise  °/0  The  Catch-All  in  the  case  of  unknown  code 

end 


% . ~ . . . . . 

70  Add  the  Date/Time  stamp  to  the  RUNS24  structure  if  a  pump  started 

% . . . . . 

[m,n]  =  size (TRACKING. VP. RUNS24) ; 
if  (STATE. VP. START) 

if  (n  ==  1)  &&  (TRACKING. VP. RUNS24(n)  ==  0) 

TRACKING. VP. RUNS24(n)  -  EVENTS.LIST (EventNum) . DATE.TIME ; 

else 

TRACKING. VP. RUNS24(n+l)  =  EVENTS. LIST (EventNum) .DATE_TIME; 

end 

end 


[m,n]  =  size (TRACKING. DP. RUNS24) ; 
if  (STATE. DP. START) 

if  (n  “  1)  &&  (TRACKING. DP. RUNS24(n)  —  0) 

TRACKING . DP . RUNS24 (n)  =  EVENTS.LIST (EventNum) . DATE.TIME; 

else 
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TRACKING . DP . RUNS24(n+l)  =  EVENTS_LIST(EventNum) . DATE.TIME; 


end 


end 


7.  End  of  Function 
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70  Non- Intrusive  Load  Monitor  (NILM) 


7o  Function  To  Perform  Leak  Analysis  of  the  CHT  System 
7o 


7o  Written  by:  Richard  Jones 


°/  v  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/a  /a  /a  /e  /a  /o  /a  /a  /a  /a  /o  /o  /o  /o  /a  /a  /o  /a  /o  /«  /o  /o  /o  /«  /a  /a  /a  /o  /a  /o  /o  /e  /o  /o  /  o  /a  /a  /o  /o  /o  /a  /o  /a  /o  /o  /a  /a  /a  /a  /a  /o  /a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a 


function  Diagnostics (n) 


y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/a  /a  /a  /a  /a  /a  /o  /a  /o  /a  /a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /a  /o  /a  /a  /o  /o  /a  /o  /a  /o  /a  /o  /o  /o  /a  /o  /o  /o  /c  /a  /a  /a  /a  /a  /a  /o  /a  /a  /o  /o  /a  /a  /o  /o  /o  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /a  /o 


7a  Global  Variables 


mmxx%y.y/ay.mxmy;mmnmxmmxmmxxmm%mmxxmxx%xmx 


global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS.LIST  EVENT. STATS  PICTURES  . . . 
TEMPRUN 


7a . . . . 

70  Initialize  process  state  variables 

7a . ~ . ~ . 

ContinueProcessDP  =  true; 
ContinueProcessVP  =  true; 
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7a  Diagnostic  Logic 

7a  If  the  CHT  system  is  determined  to  not  be  operating  correctly,  the  leak 
7a  detection  should  not  be  performed. 


7777 777777 777777777777 77777777 77 77777 7777 7 7 777 7777777777 7 77 77 77777777 77 777 7 


if  (STATE. VP. ERROR) 


ContinueProcessVP  =  false; 


end 


if  (STATE. DP. ERROR) 

ContinueProcessDP  =  false; 


end 
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%  Update  the  tracking  arrays  as  applicable 

7o . . 

if  (STATE. VP. START)  kk  (ContinueProcessVP) 
TempRun(n) ; 

end 

if  (STATE. DP. START)  kk  (ContinueProcessDP) 
TempRun(n) ; 

end 

if  (STATE. VP. STOP)  kk  (ContinueProcessVP) 
ConstructVPArray ; 

end 

if  (STATE. DP. STOP)  kk  (ContinueProcessDP) 
ConstructDPArray ; 

end 


7o7o7o7o7o7e7o7o7o7o7o7o7//o7o7o7o707//o7e7e7o767e7o7o7o7e7o7o7o7o7o7o7o7o7o7o7o7Jo7o7o7o7«7X/o7o7o7o7o7o7o7o7e7o7o7o7o7o7o7o7o7o7o707o7e7o7o7o7o7o 


%  Probe  Failure  Detection 

7o  If  there  is  a  failure  of  the  level  sensing  probes  the  in  the  Vacuum  tank, 
7o  the  Discharge  Pumps  will  cycle  on  and  off  frequently,  only  running  for  a 
70  short  period  of  time  (spikes  in  power  data)  . 


mmmmmm immmmmmmmmmmmmmmmmmm 


if  (EVENTS_LIST(n) .EVENTCODE  ==  43) 


ALARMS. PF. ACTIVE  =  true; 

[M_num,  num]  =  size (ALARMS . PF . LIST) ; 
if  (num  ==  1)  kk  (ALARMS . PF . LIST (num)  ==  0) 
ALARMS. PF. LIST (num)  =  n; 


else 


ALARMS. PF. LIST (num  +  1)  =  n; 

end 


end 


7o  Clogged  Gauge  Line  Detection 

°/e  If  the  gauge  line  becomes  clogged,  the  Vacuum  Pumps  will  run  for  an 
%  extended  time  to  lower  the  tank  vacuum  to  a  level  which  will  reach  the 
70  setpoint  of  the  control  circuit  for  Vacuum  Pump  operation. 
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7a . . . . . — 

7a  Update  the  Gauge  Clog  Array  if  VP  RunTime  >  2min  (120  Sec) 
7a . . . . . . 


EVENT. STATS . VP (n) . RUN.TIME ; 
if  (EVENT.STATS . VP (n) . RUN.TIME  >  120) 
if  (isf ield(GAUGE_CLOG,  ?ArrayO) 
[m,n]  =  size (GAUGE.CLOG. Array) ; 


if  (n  ==  5) 

for  i  =  1:4 

GAUGE.CLOG . Array (n)  -  GAUGE.CLOG . Array (n+1) ; 

end 

end 

if  n  ==  0 


n  =  1; 

end 

else 


n  =  1 ; 

end 


GAUGE.CLOG(l) .Array (n)  =  E VENT. ST ATS (n) .VP. Stop; 
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end 


y, . 

7e  Check  for  5  flags  set  in  the  Array, 
%  and  check  time  of  oldest  entry  if 
#/e  5  flags  are  present 

l . ~ . . . . 

if  (isf ield(GAUGE_CLOG,  'Array’)) 
if  (size (GAUGE_CLOG . Array)  ==  5) 


disp( ' I  '  'm  here  :)'); 


if  (GAUGE_CLOG . Array ( 1)  <  3600) 

GAUGE.CLOG . Warning ( 1) . Show  «  1; 

GAUGE.CLOG. Warning (1) .Text  =  'Clogged  Gauge  Line  Detected'; 
ProcessLeak  =  0; 

else 


GAUGE_CLOG . Warning ( 1) . Show  =  0; 
GAUGE.CLOG. Warning (1) .Text  =  '' 

end 


end 

end 
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#/0  Leak  Detection 


mmmmmnmnmmmmmmmmmmmmmmmmmn 


’/cProcessLeak  =  0; 


#/0if  (isf ield(PUMP_RUNS ,  ,VP^OFF?)) 
%  if  (PUMP_RUNS . VP_OFF  ==  1) 

7.  if  (ProcessLeak) 


LeakAnalysis(n) ; 


end 
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'/,  end 
'/.end 


*/.  End  of  Function 
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A  A  A  A  A  A  /o  A  A  /o  A  /*  /o  /o  /o  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  AAA  A  A  A  A  A  A  AAA  A  A  A  A  A  A  A  A  A  A  A  A  A  A 

X  Non- Intrusive  Load  Monitor  (NILM) 

°/0  Function  to  construct  the  array  for  the  DischargePumps  (DP)  .  This  is 
70  then  used  for  diagnostic  purposes. 

°/o 

e/0  Written  by:  Richard  Jones 

7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7 7  7 7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7  7 7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7 7 7  7  7  7 
A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  AAA  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A 

function  ConstructDPArray () 
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A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A 


7.  Global  Variables 


V  7 7  7 7 7 7 7 7 7 7 7  7  7 7  7  7 7 7 7 7 7 7  7  7 7 7  7 7 7  7 7 7 7 7  7  7  7 7 7 7 7  7  7  7 7  7 7  7  7 7 7 7 7 7 7 7  7  7  7  7  7 7  7  7  7 7 7 7  7 7 7  7  7 7 

A  A  AAA  A  A  A  A  A  A  A  A  AAA  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  AAA  A  A  A  A  A  A  A  A  A  A  A  AAA  A  A  A  A 


global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS.LIST  EVENT, STATS  PICTURES  . . . 
TEMPRUN 


7. . . . . . 

7,  Only  Perform  the  calulations  if  there  is  a  VP  running, 
7,  otherwise,  ignore  the  VP  shutting  off 

7, . . . 

if  (STATE. DP. ON  >  0) 

%  Find  size  of  arrays 

7. . 

[m,n]  =  size(EVENTS_LIST) ; 

[m,NumEvents]  =  size(EVENT_STATS.DP) ; 


7. . . . . . . 

7,  Update  the  DP  Statistics  When  the  DP  is  Secured 

r~ . — . . 

EVENT.STATS. DP (NumEvents+1) .START  =  TEMPRUN (1) .DP. TIME; 
EVENT_STATS.DP(NumEvents+l) .STOP  =  EVENTS_LIST(n) . DATE.TIME; 
EVENT.STATS. DP (NumEvents+1) . START.FILENAME  =  TEMPRUN . DP (1) .FILE; 
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EVENT.STATS. DP (NumEvents+1) . STOP.FILENAME  =  EVENTS.LIST(n) .FILENAME; 


7, . 

7«  Compute  the  Pump  Runtime 

% . 

Delta  =  datevec (EVENT_STATS . DP (NumEvents+1) . STOP  -  ... 

EVENT. STATS. DP (NumEvents+1) .START) ; 

EVENT. STATS . DP ( NumEvent s  + 1 ) . RUN.TI ME  =  Delta(5)  *  60  +  Delta(6) ; 

7. . 

7o  Compute  the  time  between  DP  runs 

7. . . . 

if  (NumEvents  >  1) 

Delta  =  datevec ( EVENT. STATS . DP (NumEvents+1) .START  -  ... 

EVENT.STATS . DP (NumEvents) . STOP) ; 

EVENT.STATS. DP (NumEvents+1) .OFF. TIME  =  Delta(4)/60  +  Delta(5)  +  Delta(6)*60; 

else 

EVENT.STATS. DP (NumEvents+1) .OFF. TIME  =  0; 

end 

7,  Reset  the  tracking  variables 

7, . — . . . — 

STATE. DP. ON  =  0; 

STATE . DP . STOP  =  false; 

else 

STATE . DP . STOP  =  false; 

end 
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7©  End  of  Function 
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%  Non- Intrusive  Load  Monitor  (NILM) 


7,  Function  to  construct  the  array  for  the  Vacuum  Pumps  (VP)  .  This  is 
°/0  then  used  for  diagnostic  purposes. 

7a 


l  Written  by:  Richard  Jones 


function  ConstructVPArrayO 


7  7  7  7 7 7  7  7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7 7  7  7  7  7 7  7  7 7  7  7  7  7  7  7  7  7 7  7  7 7  7 7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7 

/a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /a  / a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  h  h  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a 


#/o  Global  Variables 


7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7 7  7  7  7  7  7 7 7  7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7 7  7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7 

/a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  / a  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a 


global  PARAMETERS  STATE  ALARMS  TRACKING  EVENTS_LIST  EVENT, STATS  PICTURES  . . . 
TEMPRUN 


7. . . . 

7.  Only  Perform  the  calulations  if  there  is  a  VP  running, 
7.  otherwise,  ignore  the  VP  shutting  off 

7." . . . 

if  (STATE. VP. ON  >  0) 

7. . 

7.  Find  size  of  arrays 

7. . . . . 

[m,n]  =  size (EVENTS.LIST) ; 

[m.NumEvents]  =  size (EVENT.STATS . VP) ; 


7. . . . . . .  . 

7i  Update  the  VP  Statistics  When  the  VP  is  Secured 

7. . . . . . . . . . — 

7. . 

7.  Pump  #1 
% - 
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EVENT. STATS . VP (NumEvents+1) .START  =  TEMPRUN . VP( 1) . TIME; 

EVENT.STATS. VP (NumEvents+1) .STOP  =  EVENTS.LIST(n) . DATE.TIME; 
EVENT.STATS . VP (NumEvents+1) . START.FILENAME  =  TEMPRUN. VP (1) .FILE; 
EVENT.STATS. VP (NumEvents+1) . STOP.FILENAME  =  EVENTS_LIST(n) .FILENAME; 


%  Compute  the  Pump  Runtime 

% . . . 

Delta  =  datevec (EVENT.STATS . VP (NumEvents+1) . STOP  -  ... 

EVENT.STATS . VP (NumEvents+1) .START) ; 

EVENT. STATS . VP (NumEvents+1) . RUN.TIME  =  Delta(5)  *  60  +  Delta(6) ; 


% . . . . . . 

%  Compute  the  time  between  VP  runs 

% . 

if  (NumEvents  >  1) 

Delta  =  datevec (EVENT.STATS. VP (NumEvents+1) .START  -  ... 

EVENT. STATS . VP (NumEvents) . STOP) ; 

EVENT.STATS. VP (NumEvents+ 1) .OFF.TIME  =  Delta(4)/60  +  Delta(5)  +  Delta(6)*60; 

else 

EVENT.STATS. VP (NumEvents+1) .OFF.TIME  =  0; 

end 


•/. . 

'/.  Pump  #2 

7. . 

if  (STATE. VP. ON  ==  2) 

EVENT.STATS . VP (NumEvents+2) . START  =  TEMPRUN . VP (2) . TIME ; 

EVENT.STATS . VP (NumEvents+2) . STOP  =  EVENTS. LI ST (n) . DATE.TIME ; 
EVENT.STATS. VP (NumEvents+2) . START.FILENAME  =  TEMPRUN. VP (2) .FILE; 
EVENT.STATS. VP (NumEvents+2) . STOP.FILENAME  =  EVENTS.LIST(n) .FILENAME; 
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°/o  Compute  the  Pump  Runtime 

Delta  =  datevec (EVENT.STATS . VP (NumEvents+2) . STOP  -  ... 

EVENT.STATS. VP (NumEvents+1) .START) ; 

EVENT.STATS. VP (NumEvents+2) .RUN.TIME  =  Delta(5)  *  60  +  Delta(6) ; 

°/o  Compute  the  time  between  VP  runs 

Delta  =  datevec (EVENT.STATS . VP (NumEvents+2) . START  -  ... 

EVENT.STATS. VP (NumEvents) .STOP) ; 

EVENT.ST ATS. VP (NumEvents+2) .OFF.TIME  =  Delta(4)/60  +  Delta(5)  +  Delta(6)*60; 

end 


t 


/e  Reset  the  tracking  variables 


7. 


STATE. VP. ON  *  0; 

STATE. VP. STOP  =  false; 


else 


STATE . VP . STOP  -  false; 


end 


7  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/«  /o  /e  /e  /e  /o  /«  /o  /o  /o  /e  /o  /«  /g  /o  /o  /o  /o  /o  /e  /o  /a  /o  /o/e/o  /o  /o  /o  /9  /o  /o  /o  /e  /o  /«  /o  /e  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /e  /o  /o  /o  /o  /o  /«  /«  /«  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


'/,  End  of  Function 


mmitmmmmmmimmmmmmmmxmmmmmxiimra 


%  Non-Intrusive  Load  Monitor  (NILM) 

70  Function  To  Backup  All  Data 


7o  Written  by:  Richard  Jones 


V  7  V  7  7  7 7 7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7 7 7  7 7  7  7  7  7  7  7  7  7  7  7 7 7  7  7  7  7  7  7  7  7 7  7 7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7 

/a  /o  /a  /o  /a  /a  /a  /o  /a  /o  /o  /o  /a  /o  /o  /c  /e  /a  /a  /o  /#  /o  /a  /o  /o  /o  /a  /a  /a  /a  /e  /o  /c  /o  /a  /a  /o  /o  /a  /a  /o  /o  /o  /a  /e  /e  /o  h  h  / o  /a  /a  /a  /a  /o  /o  /o  /a  /o  /a  /a  /o  /o  /e  /o  /a  /a  /o  /a  /o  /a  /o  /o  / o  /a 


function  BackupDataO 


°/e  Global  Variables 


7  7  7  7  7 7  7  7  7  7  7 7  7 7 7  7  7 7  7  7  7  7  7  7  7 7 7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7 

/a  /o  /a  /o  /a  /a  /o  /o  /a  /a  /a  /a  /a  /a  /o  /a  /o  /a  /a  /o  /a  /o  /a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a 


global  EVENTS. LI ST  EVENT.STATS  EVENTS.SIZE  EVENTLIST  EVENTTIMES  VPRUN  . . . 
DPRUN  PUMP.RUNS  SHOW.WARNING  TOTAL.PUMP.RUNS  SLOW.LEAK 


%%%%%%%%%%%%%%%%%  %%%%%%%%%%%%%%% /a 


°/e  Save  EVENT.STATS  data  structure 


7  7  7  7 7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7 7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7 7  7  7  7  7  7 7  7  7  7  7  7  7  7  7 

/a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  / a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /o  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a 


save  DataBackup/BackupEventStats .mat  EVENT.STATS; 
save  DataBackup/TotalPumpRuns . mat  TOTAL.PUMP.RUNS ; 
save  DataBackup/SlowLeak .mat  SLOW.LEAK; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


7e  End  of  Function 


7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7  7  7  7  7  7  7 7  7  7  7  7  7  7  7 7  7  7  7  7  7 7  7  7  7  7  7  7 7  7  7  7  7  7 7  7  7  7  7  7  7  7  7 

/o  /a  /o  /o  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /«  /a  /a  /a  /o  /a  /a  /a  /a  /a  /o  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /a  /o 
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